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FOREWORD 


This  report  was  prepared  in  the  Components  3ranch  (AFTC),  Turbine  Engine  Division.  Air 
Force  Aero  Propulsion  Laboratory.  Wright- Patterson  Air  Force  Base.  Ohio,  under  Project 
3066,  “Gas  Turbine  Technology.”  Task  306803.  “Advanced  Engine  Studies.”  with  Charles  E. 
Bentz  as  Project  Engineer. 

This  report  covers  work  conducted  within  the  Components  Branch  in  the  time  period  between 
July  1965  and  June  1967  and  was  submitted  by  the  author  31  August  1967. 

This  technical  report  has  been  reviewed  and  is  approved. 


ERNEST  C.  SIMPSON 
Chief,  Turbine  Engine  Division 
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ABSTRACT 


This  report  describes  a  digital  computer  program  titled  SMOTE  (Simulcticn  of  Turbofan 
Engine).  SMOTE  is  a  computer  program  for  balancing-cycle  turbofan  engines  capable  of  running 
both  design  and  off-design  points.  The  program  is  written  in  Fortran  IY  language  and  was 
designed  for  use  on  an  IBM  7050  Digital  Computer,  although  it  has  also  bean  run  on  an  IBM 
System  3eo.  Performance  maps  (Bloch  Data  format)  of  V  6  major  engine  components  are 
required.  Information  for  setting  up  the  Block  Data  and  ir-cut  data  is  given  in  the  report.  Also 
included  in  the  report  is  a  complete  program  listing  w  description  of  each  subroutine  and 
a  sample  data  pack. 


(Distribution  of  the  abstract  is  unlimited.) 
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SYMBOLS 


STATION  NUMBERS 


1 

2 

21 

3 

4 

5 
55 

6 

7 

8 
9 

23 

24 
22 
23 
29 


ambient 
fan  entrance 

fan  exit/ compressor  and  duet  entrance 

cam  pressor  exit/burncr  entrance 

burner  exit/hi  pressure  turbine  entrance 

hi  press,  turbine  exit/lo  press,  turbine  entrance 

lo  proii.  ^rblne  exit 

-fterbemer  entrance 

afterburner  exit 

main  rc*xla  threat 

main  north*  exit 

duct  burner  entrance 

duct  burner  exit 

duct  exit  (if  mbaed-flew  engine) 

duct  nozzle  threat 

chct  nozzle  exit 


THERMCaJYNAlnC  PROPERTIES 

T  total  temperature  (degrees  R) 

H  total  sijthzlpy 

P  total  pressure  (atmospheres) 

S  total  entropy 

TS  eistic  temperature  (degrees  R) 

PS  static  pressure  (atmospheres) 

FAR  fuel-air  retie 

AM  Msch  number 

V  Telocity  (rset/secood) 


COMPONENT  SYMBOLS 


F 

C 

B 

COM 

T3? 

TLP 

A 

AFT 

D 

BUC 

U 

NOZ 

OB 

T 


fan 

compressor 

burner 

burner 

hi  pressure  turbine 

In  pressure  turbine 

afterburner 

afterburner 

duct 

duct 

mats  nozzle 
nozzle 
overboard 
total 


ENGINE  SYMBOLS 


ETAB 

CN 

PCN 

PR 

Z 

ETA 

TFF 

BKTC 

DHT 

VfA 

WF 

M3 

PCBL 

BL 

HFEXT 


ram  recovery 
corrected  speed  ratio 
percent  speed 
pressure  ratio 
pressure- rsiio  ratio 
efficiency 

turbine  flow  funcUco 

tnrfcias  delta  enthalpy  (temperature  corrected) 

turbine  delta  enthalpy 

air  flow  (potmds/sec) 

fa el  flow  (pocnds/eec) 

gas  flow  (pouDds/sec) 

percent  blood 

bleed  (pcuads/eec) 

horsepower  extracted 
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SYMBOLS  (CONT) 


DP 

pressure  drop 

DT 

temperature  increase 

MISCELLANEOUS 

A 

area 

ALT? 

altitude 

AM 

Mach  number  of  aircraft 

BYPASS 

bypass  ratio 

CF 

correction  factor 

CS 

ambient  speed  of  sound 

CV 

nozzle  velocity  coefficient 

DEL 

delta  degradation  coefficient 

DS 

design  value 

DUM 

dummy  (not  used) 

FG 

gross  thrust 

PGM 

momentum  ihrust 

PGP 

pressure  thrust 

FN 

net  thrust 

FRD 

ram  drag 

GU 

initial  or  guess  values 

ITRYS 

number  of  loops  thru  engine  before  quitting 

SFC 

specific  fuel  consumption 

TOLALL 

tolerance 

VA 

velocity  of  aircraft 

VJ 

jet  velocity 

NOTE:  Soma  symbols  may  be  truncated  when  combined  with  other  symbols  due  to  six 

character  limit  imposed  by  Fortran  Computer  Language. 
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SECTION  I 
INTRODUCTION 
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Part  n 

SECTION  n 
BLOCK  DATA 

To  provide  the  basic  component  performance  requires  the  maps  for  the  fan,  compressor, 
Combustor  end  both  turbines  to  be  converted  to  Block  Data  subroutines.  The  Block  Data  pre¬ 
sented  in  this  report  is  very  general  aod  does  not  represent  any  specific  engine.  For  the  follow¬ 
ing  discussion  on  setting  up  Block  Data,  refer  to  the  program  listing.  Section  V.  The  format 
for  all  the  maps  is  very  similar,  with  a  maximum  of  15  curves  and  a  maximum  of  15  pointe  for 
each  curve.  “N"  refers  to  the  number  of  curves  and  “NP”  is  an  array  of  the  number  of  points 
on  each  curve. 

For  the  feu  and  compressor  maps,  the  corrected  speed  (CN)  is  input  as  a  ratio  of  the  actual 
corrected  speed  curve  to  the  design  corrected  speed.  The  pressure  ratio  (PR),  corrected  air¬ 
flow  (WAC).  and  efficiency  (ETA)  are  input  in  groups,  with  a  group  for  each  corrected  speed 
curve.  The  corrected  speeds  and  pressure  ratios  must  be  loaded  low  to  high.  Note  that  there 
are  two  points  per  card. 

For  the  eominctor  map,  curves  of  pressure  entering  the  combustor  (P3)  in  PSI  are  loaded 
low  to  high,  as  well  as  delta-T  (T4  -  T3)  points  for  each  curve.  The  efficiencies  (ETA)  cor¬ 
responding  to  each  DELT  are  loaded  in  one  array.  Note  that  in  this  array,  each  column  is 
obtained  from  one  pressure  curve. 

For  the  turbine  maps,  turbine  flow  function  (TFF)  curves  are  input  low  to  nigh,  and  the 
corrected  speed  (CK),  delta-3  (DH),  and  efficiency  (ETA)  are  input  in  groups,  with  a  group  for 
each  TFF  line.  The  corrected  speed  must  be  loaded  low  to  high.  TFF,  CN,  and  DH  are  defined  as 


TFF  ■ 


CN  = 


WC-  yr5T 
?in 

PCNC 


D3  = 


«/tIN 

hIN  "  HOUT 


■IN 


Formats  other  than  those  presented  here  may  be  used  for  the  component  performance  maps. 
However,  some  change  In  logic  would  then  be  required. 
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SECTION  in 
INPUT  DATA 

,v  -  ,  )•  £  ■- ,  ■  -.  '  ■ 

The  input  data  Is  divided  into  two  sections;  data  cards  for  the  controlled  output,  and  data 
cards  in  Namelist  format  for  running  each  point.  For  the  following  discussion  on  setting  up 
the  input  data,  refer  to  the  listing  of  sample  data  Immediately  following  the  program  listing. 

1.  CONTROLLED  OUTPUT 

The  variables  that  are  to  be  output  are  selected  by  the  flr3t  section  of  data  cards.  Any  vari¬ 
able  that  is  in  one  of  the  main  commons  (DESIGN,  FRONT,  SIDE,  or  BACK)  may  be  selected 
for  output  by  punching  the  name  of  the  variable  as  it  appears  in  the  common  (with  trailing 
blanks,  if  necessary)  in  Columns  1  through  6.  Up  to  102  variables  (seventeen  lines  of  six  vari¬ 
ables)  may  be  chosen  for  a  particular  run.  During  the  output  phase,  the  name  of  the  variable  is 
printed  out,  with  its  value  printed  immediately  below  the  name. 

Another  feature  of  the  controlled  output  is  the  ability  to  change  the  name  of  a  variable  to  be 
output;  for  example,  it  may  be  desired  to  change  a  station  designation  to  one  more  common  to 
a  particular  programmer.  In  this  case,  the  variable  name  would  be  punched  in  Columns  1 
through  6  as  described  above,  but  in  addition,  the  desired  name  would  be  punched  in  Columns  13 
through  18.  Special  symbols,  such  as  /,  may  be  used  in  the  new  name.  The  last  card  of  the 
controlled  output  must  be  a  card  with  “THEEND”  punched  in  Columns  1  through  6. 

2.  NAMELIST  INPUT 

The  normal  data  for  running  the  dcrired  points  follows  the  controlled  output  data  and  is  in  a 
Namelist  format,  where  the  name  of  the  Namelist  is  DATAIN.  Usually  the  first  set  of  data  is 
the  design  point,  as  shown  in  the  sample  Input  data.  When  the  design  point  is  run  (IDES-1),  all 
map  scaling  or  correction  factors  are  printed  out, as  well  as  being  retained  in  common.  There¬ 
fore,  it  is  possible  to  run  off-design  points  immediately  following  the  design  point  by  making 
"se  of  the  values  in  common,  or  to  begin  running  an  off-design  point  immediately  by  Inputting 
tiie  scaling  or  correction  factors.  The  first  method  is  usually  easier,  but  the  second  method  may 
be  desired  if  many  points  are  to  be  run  using  the  same  engine  parameters  with  no  changes 
except  for  power  setting,  Mach  number,  and  altitude. 

The  controls  which  fix  the  type  of  turbofan,  the  mode  of  operation,  the  method  of  calculating 
ram  recovery;  etc.,  are  explained  In  the  listing  of  subroutine  INPUT.  Subroutine  ZERO  deter¬ 
mines  what  values  in  common  will  be  zeroed  between  points.  None  of  the  design  values  or 
correction  factors  are  ever  zeroed. 

A  title  card  must  be  input  immediately  after  the  first  point  of  the  data  pack,  and  ITITLE  must 
be  set  equal  to  1  in  the  data  for  the  first  point.  This  is'  because  a  title  is  always  printed  for 
each  point  and  must,  therefore,  be  previously  defined.  The  input  format  for  the  title  If  12A6, 
and  the  resulting  72  spaces  are  centered  on  the  page  when  printed  out.  The  title  maybe  changed 
by  setting  ITITLE  =  1  and  inserting  a  new  title  card  after  the  Namelist  data  for  foe  point. 

When  an  afterburning  cr  duct-burning  point  is  run,  the  exhaust  nozzle  areas  are  allowed  to 
float  to  obtain  optimum  expansion.  This  means  that  there  can  be  no  balancing  at  the  point,  and  it 
is  necessary  to  prebalance  the  engine  cyole  in  a  nonaugmented  mode.  That  is,  an  identical 
point,  except  that  it  is  nonaugmented,  must  be  run  before  either  afterburning  or  ductburning. 
When  either  IAFTBN  or  IDBURN  is  greater  than  zero,  the  program  will  automatically  set 
IN  IT  =  1  and  use  the  balanced  values  from  the  preceding  point.  The  nozzle  areas  are  returned 
to  their  standard  design  values  after  completing  an  augmented  point.  Some  examples  of  this 
type  of  afterburning  are  given  in  the  sample  data  listing. 
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SECTION  IV 


SUBROUTINE  DESCRIPTIONS 


The  following  will  be  a  brief  description  of  what  each  subroutine  calculates  or  controls. 

See  Figure  1  for  a  subroutine  flow  chart 

SMOTE  Dummy  main  program  to  initiate  the  calculations  and  cause  the  input  of  the 
controlled  output  variables.  Because  of  the  looping  between  subroutines ,  control 
is  never  transferred  back  to  this  routine. 

ENGBAL  Main  subroutine.  Controls  all  engine  balancing  loops,  checks  tolerances  and 
number  of  loope,  and  loan*"  natrir. 

MATRIX  Solves  error  matrix. 

INPUT  Bead3  Namelist  data  and  title.  Trints  title. 

ZERO  Zeroes  common  and  certain  controls. 

COINLT  Determines  ram  recovery  xnd  performs  Inlet  calculations. 

ATM062  1952  ARDC  At? nosphere  Tables. 

RAM  Calculates  ram  recovery  defined  by  MIL-E-SOOSB  Specifications. 

GOES  Determines  Initial  values  of  independent  variables  (FCNF.  PCNC,  end  T4)  at 

each  point.  It  may  be  desired  to  change  these  equations  to  suit  a  particular 
engine.  The  closer  the  initial  values  are  to  the  final  values,  the  faster  the 
program  will  balance. 

CCFAN  Uses  Block  Data  to  perform  fan  calculations. 

CGCOMP  Uses  Block  Data  to  perform  compressor  calculations. 

COCOMB  Uses  Block  Data  to  perform  combustor  calculations.  May  use  either  T4  or 
WFB  as  the  main  parameter. 

CGHFTB  Uses  Block  Data  to  perform  high  pressure  turbine  calculations.  Calculates 
ERR(i)  and  ERR(2). 

COLPTB  Uses  Block  Data  to  perform  low  pressure  turbine  calculations.  Calculates  ERR(3) 
and  ERR*4). 

FRTOSD  Dummy  routine  to  transfer  values  from  common  FRONT  to  common  SIDE. 

CODUCT  Performs  duct  and  duct-burning  calculations.  May  use  either  T24  of 'jVFD  as  the 

main  parameter  for  duct-burning'.  Controls  the  duct  nozsle  and  calculates  ERE/5) 
if  in  separate-flow  mode. 

FASTBK  Dummy  routine  to  transfer  values  from  commons  FRONT  and  SIDE  to  common 
BACK. 
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COMIX 

COAFBN 

COMNOZ 

PERF 

OUTPUT 

COKOUT 

ERROR 

SYG 

TAPES 

THCOMP 

THTURB 

THERMO 

PROCOM 

SEARCH 

MAPBAC 

CONVRG 

CONDIV 

AFQUIR 


Performs  gas-mixing  calculations  if  in  mixed-flow  mode.  At  design  points  it 
calculates  areas  from  either  an  input  static  pressure  (P355)  or  an  input  Mach 
cumber  (AM  55)  if  PS55  =  0.  At  off-design  points  it  calculates  static  pressures 
and  Mach  numbers  from  the  design  areas.  Calculates  ERR(-5). 

Performs  after-burning  calculations.  May  use  either  T7  or  WFA  as  the  main 
parameter. 

Controls  the  main  nozzle  and  calculates  ERR(6). 

Calculates  performance  after  the  engine  ia  balanced. 

Prints  output  except  for  controlled  output.  Prints  the  main  commons  in  a  close 
format  after  each  point. 

Controls  and  prints  the  controlled  output  variables. 

Controls  all  printouts  if  an  error  occurs.  Prints  name  of  subroutine  where  error 
occurred  and  also  prints  the  values  of  all  variables  in  the  main  commons. 

Controls  printing  from  UNTT03.  Throughout  the  program  and  particularly  In 
ENGBAL.  certain  messages,  variables,  and  mairix  values  arv  written  on  UNIT08 
as  an  aid  In  determining  why  an  error  occurred  or  why  a  point  did  not  balance. 
These  values  ace  printed  out  if  subroutine  E  RROR  is  called  and  IDUMP  is  great¬ 
er  than  zero,  or  after  a  good  point  if  IDUMP  =  2. 

Defines  UNIT08,  which  is  just  a  “scratch”  disk  and  does  not  require  a  SSETUP 
card.  Normal  input  and  output  are  on  UNIT05  and  UNITOS,  respectively. 

Performs  isentropic  calculations  for  compressors. 

Performs  isentropic  calculations  for  turbines. 

Provides  thermodynamic  conditions  using  PROCOM. 

Calculates  thermodynamic  gas  properties  for  either  air  or  a  fuel-air  mixture, 
based  on  JP-4. 

General  table  look-up  and  interpolation  routine  to  obtain  data  from  the  Block 
Data  subroutines. 

Used  when  calculations  result  in  values  not  on  the  turbine  maps.  Changes  the 
map  value  and  an  independent  variable  \PCHF,  PCNC,  or  T4)  in  an  attempt  to 
rectify  the  situation. 

Performs  nozzle  calculations  for  a  convergent  nozzle. 

Performs  nozzle  calculations  for  a  convergent-divergent  nozzle. 

General  quadratic  interpolation  routine. 
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FANDAT  Block  Data  for  fan. 

CMPDAT  Block  Data  for  compressor. 

CMBDAT  Block  Data  for  combustor. 

HPTDAT  Block  Data  for  high  pressure  turbine. 

LPTDAT  Block  Data  for  low  pressure  turbine. 


SECTION  V 

PROGRAM  LISTING 

The  following  is  a  complete  listing  of  all  subroutines  required  to  run  SMOTE.  Immediately 
following  the  program  listing  is  a  listing  of  a  sample  input  data  pack. 
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4I8FTC  SHOTE  DECK, M94/2,XR7 
COMMON  /POINT/I DATPT 
I0ATPT=O 
CALL  CONOUT(l) 

CALL  ENGBAL 

STOP 

END 
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SlttFTC  ENGBAi,  0ECXyH94/2yXk? 
SUBROUTINE  ENGBAI 
COMMON  /  At*./ 

1WORO  , IDES  .JOES  ,KDES 
2ISASHX  , IOBURN , I AFTBN, I OCD 


,NUOE 
? 1 MCD 


}NUKMAP  yMAPSDG 


31TRYS  tLOOPERf NOMAP 
COMMON  /DESIGN/ 

1 PCNFGU , PCNCGU, T4GU 
2ZFOS  .PCNFOS.PRFDS 
3ZCDS  ,PCNCOS,PRCDS 
AT  AOS  rVFBDS  ,  DTCL'DS,  ETABDS,  H.43CD5 
5TFHPOS,CNK?OS»ETHPOS.TFHPCF,CNHPCF 
6TFLP0S,CNLP0S,ETLPnS,TFLPCF,CNLPCF 
TT  2-505  ,WFDOS  ,D7DUDS,ETAODS,kA23DS 


yCUMOl  ,DUMD2 
«  ETAfDSt  WAFOS 
.ETACDS.kACDS 


8T70S 

,«FADS  yOTAFDS 

, ETAAOS 

t MG6CDS 

yOPAFOS 

.DTAFCF 

yETAACF, 

9A55 

yA25  ,A6 

,A? 

,A8 

y  A9 

,A28 

yA29  , 

APS55 

COMMON 

,AM55  yCVDNOZ 
/  FRONT/ 

.CVHNOZ 

« A8SAV 

y A9SAV 

.A28SAV 

.A29SAV 

1T1 

yPl  yHl 

ySl 

*T2 

y?2 

yH2 

yS2  y 

2T21 

,P21  ,H21 

tS21 

«T3 

»P3 

,H3 

yS3  y 

3T4 

,P4  ,H4 

,S4 

yT5 

yP5 

yH5 

yS5  , 

4755 

yP55  ,H55 

yS55 

,elf 

,3LC 

yBLDU 

yBLOB  y 

5CNF 

,PRF  yETAF 

yWAFC 

ykAF 

ykA3 

,WS4 

,FAR4  » 

6CNC 

,PRC  yETAC 

ykACC 

,kAC 

yETAB 

yOPCOM 

yOUMF  » 

7CNHP 

,eta7hp,ohto:p 

yOHTC 

•  BLHP 

yWGS 

»PAR5 

.CS  y 

8CNLP 

,£7A7L?,0H7CLP 

,OHTF 

yfiLLP 

,kG55 

y FAR 55 

,HP£X7  y 

9AM 

yALTP  yElAR 

,ZF 

yPCNF 

,ZC 

yPCHC 

yWFB  » 

ATrFHP 

,TFFLP  yPCBLF 

yPCSLC 

yPCBLDU 

I.PCBLOB 

.PC3LHP 

-eceu? 

COMMON  /  SIDE/ 


XXP1 

fXkAr 

yXkAC 

»XBLF 

yX8LCU 

,XH3 

y  DU  MSI 

,  DU  MS  2 

t 

XXT21 

yXPZl 

•  XH21 

yXS21 

yT23 

,P23 

,H23 

,S23 

t 

3T24 

,P24 

,H24 

,S24 

*725 

,P25 

,H25 

,525 

t 

4723 

,P23 

yH28 

y  528 

*729 

,P29 

,H29 

tSZ9 

1 

5HA0 

ykfO 

»VcG24 

,FAR24 

,ETAO 

yDPDUC 

.BYPASS, DUMS3 

t 

6TS28 

,PS23 

yV28 

y  AM28 

*TS29 

»PS29 

,V29 

,AM29 

COMMON 

/  BACK/ 

XXT55 

,XP55 

yXH55 

,XS55 

,X725 

,XP25 

,XH25 

,X  S25 

t 

XXkFB 

yXMG55 

yXr AR55 

,XkFD 

,XkG24 

,XFAR24 

yXXPl 

.DUMB 

f 

i  3T6 

»P6 

»H6 

»S6 

,77 

?P7 

»H7 

,S7 

f 

i  478 

,P8 

yK8 

tS8 

,79 

*P9 

,H9 

»S9 

f 

5WG6 

.WFA 

.WG7 

y  FAR7 

,ETAA 

,DPAFT 

,Y55 

»V25 

t 

6?S6 

yV6 

.AM6 

yTS7 

*PS7 

,V7 

,AH7 

,AM25 

f 

7TS8 

,PS8 

yV3 

« AMS 

»TS9 

,PS9 

»V9 

y  AM9 

f 

8VA 

,FRD 

yVJD 

yFGNO 

,Vj« 

,FGKM 

yFGPO 

yFGPM 

J 

9FGM 

,FGP 

,kFT 

,  WGT 

,FART 

,fg 

,FN 

,SFC 

DIMENSION  VAP.  (6  )  yOEL  (61  y  EFRBI6I  »DELVAR{6) 

,£MAT (6,6 J  yVMAT (6) y 

INIT  ylOUHP  , 
IDSHOCylMSHOCy 
TOLALL  yERR (6) 

DELFG  ,DEL5N 
PRFCF  ,ETAFCF 
PRCCF  yETACCF 
OPCODS  yOTCOCF 
ETHPCF  yOHHPCF 
ETLPCF  yOHLPCF 
OPOUDS yDTOUCF 


IAMTP  , 
NOZFLT, 


OELSFCy 
WAFCF  , 
kACCF  y 
ETABCFy 
T2DS  . 
T21DS  y 
ETADCF, 


1AMATI6) 

DATA  AU0K0/6HENG5AL/ 
CALC  INPUT 

IF! INIT.EO.I J  GO  TO  50 
TFFK?=TFKPOS 

tfflp=tflpos 

50  LC0PER=0 

NUHHAP=0 
1  LOQP=0 

NOMAP =0 
IGG=2 

30  2  1=1 r6 
VMATIl }=0. 

AMAT { I 1=0. 

OELVARU  1=0. 


AFAPL— T?-£7=I2= 
Partn 


DO  2  L-lib 

2  cRATIl,L>=0. 

3  L‘J0P£R=lG0?6ft-H 
CAl  L  C-OfAN 
WO«D=AtiOKO 
If (LGOPEK.GT.ITRYS)  GO  70  1ft 
IP(NO«AP.GT.Oi  GO  TO  1 
KUKMAPeg 

VAR{ 1 )=2F*I00. 

VAR(2)=PCNF 

VARi3)=ZC*lC0. 

IF(H0DE.£0.0.0K.H)DE.S0.2)  VAR(4)=PCRC 
If (MODE .£0.1)  VAR(4)»T4/10. 

VAR  15  )*TFFhP 
VAR(6i*TffLP 
DO  4  1=1,6 

4  !F(A8$(£RR(1 ) )  .GT.TOLALL )  GO  70  5 
CALL  PERf 

CALL  ERROR 

5  1FIL00P.GT.0)  GO  TO  7 

NAP£DG=0 

MAPSET*0 
DU  6  1=1,6 
ERR8II J=ERRII ) 

6  DEL  1 1 )=0.00l«VAK( I J 
GO  TO  9 

7  IF(p.APEUG.EO.U)  GO  TO  70 
MAPE0G=0 

?*AP5ET*1 

VAR (LOOP  I =VAR  C  LOOP  J  +2 . *OEL (LOOP  J 
GO  TO  10 

70  IF(MAPSET.EO.O)  VAR (LOO?i=VAR(LOOP)+DEL (LOOP) 

IF (hAPSET.EO.l I  VAk (LOOP I = VAR (LOOP)— DEL ( LOUR  1 
*APSET=0 
DO  8  1=1,6 

8  6*»AT (  I  ,LOQP)= (£KR8{  I )— EP.K  (I  )  I/OEL  (LOOP) 

9  LOOP=LOOP+i 

IF  (LOOP. GT .6  )  C-0  TO  11 

VAR (LOOP ) =VAR ( LOOPI-DEL (LOOP! 

ID  ZF=VAR(1)/100, 

PCNF=VAK(2) 

ZC=VAX(3I/100. 

iF(*OOF.EO.O.OK.HUUE.E0.2)  PChC=VAR(4) 

IF (ROUE. EO. 1 )  T4=VAR(4)*10. 

TFFHP=VAR(5) 

TFFLP=VAR(6) 

IFIZF.LT. 0. )  ZF=0.05 
IF (ZC.LT.G. )  ZC=0.05 
GO  TO  ( 1 ,31, 1G0 

11  OU  12  1=1.6 

12  A*AT(I >=-ERRB(l) 

DO  14  1=1,6 
IZER0=0 

1)0  13  LOOP =1.6  , 

13  IF<eKAT(!,LOO?)-.EV.O.)  IZEKO=I ZcRU+1 
IF ( I ZERU.LT .6)  GO  TO  14 

WRITE (6,100! I 
IOOP£K=IT«YS+100 
GO  TO  18 

14  CONTINUE 

00  16  LU0P=!,6 


20 


AFAPL-TR-67-125 
Part  H 


I2F-K0=0 
UU  15  1  =  1,6 

15  IF  (6MAT{  I  ,LHi)P  >.EU,0.  )  I Z  ERO=  I  ZERO+1 
IF  (1 ZERO.LT .6 }  GO  TO  16 

KK 1 76(6,1011  LOOP 
LOOPER=I TRYS+100 
GO  TO  18 

16  CONTINUE 

CALL  hATR I X ( EHAT, VMAT , AMA  F 1 

DO  17  L00P=i,6 

DEL  VAR  (LOOP  )  =0.6<=VMAT  (LOOP) 

IF(A8S(0ELVAR(L00P) I .GT .0 .05$ VAR ( LOOP ) ) 

10ELVAR (LOOP ) =0.0 5- VAR (LOOP )~DEL VAR { LOOP } /A3S (DEL VAR  » LOUP ) ) 

1 7  VAR ( LOOP ) =VAR ( LOOP  J  +DEL VAR ( LOOP ) 

18  WRITE (8,102) LOOPER 
DO  19  1=1,6 

19  WRITE (8 , 103) AKAT ( I  I , ( EHAT ( I , L ) , L=1 , 6} , VMAT ( I) , DEL VAR ( I ) ,VAR(I ) 
ERRAVE=(ABSCERRB(1)  J+A8S(ERRB(2n+A8S(fcKRB(3J}  + 

1  ABS ( ERRS ( 4) l+ABS  C  ERR3( 5 ) J+A8S (ERRB (6) ) J/6. 

DELAV6= ( ABS ( OELVAk ( 1 ) I+ABS ( DELVAR (21) +ABS ( DEL VAR ( 3 ) )  + 

1  ABS ( DELVAR (4) 1 +ABS (DEL VAR { 5 ) J+ABS ( DELVAR (6 } J  j /6. 

WRITE (8,1 04 )ERRAVE,D£LAVE 
IG0  =  1 

IFtLOOPER.LT. ITRYS)  GO  TO  10 

CALL  ERRGR 

RETURN 

100  FORMAT (4H0R0W, 12, 16H  IS  ZERO  IN  EMAT) 

101  FORMAT (7H0C0LUMN,I2» 16H  IS  ZERO  IN  EMAT 1 

102  FORMAT (8KB  ERRB»28X23HERR0R  MATRIX  AFTER  LOUP , I4,29X4HVMAT, 
16X6HDEL VAR,7Xi4HVARI ABL  ESSSSSS  J 

103  FORMAT  (  1H0tF8.4,3X6F10.4t  10XF10,4,F11  .4 ,4XF1 1  .4,60556^35 i  - 

104  FORMAT ( 1H0,F8.4,32X14HAVERAGE  VALUES, 42XF11. 4, 6HS5SSSS) 

END 
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AFAPL-TR-67-125 
'Part  II 


5IBFTC  MATRIX  DECK, *94/ 2, XR 7 
SUBROUTINE  MATRIXC£,V,A) 

DIMENSION  E{6,6),V{6),fi{6)»PIV(7),T{6,7) 
DO  I  1=1,6 
TUf?)=Am 
DO  1  J=I,6 

1  T(I,JI=ECI,J) 

DO  7  1=1,* 

TEHP=0. 

DO  2  J=I ,6 

IFiTENP.GT.ABSCTIJ, U)}  GO  TO  2 

TEMP=ABS(T { J  »  J U 

IPIV=J 

2  CONTINUE 
IP1=I+1 

DO  3  J=iPl ,7 

3  PlV(Ji=TnPIV,JI/TUP!VfI  ; 

IFR0M=6 

IT0=6 

4  IF  UFKUH.FO.IPIV)  GO  TO  6 
RM=— T ( I FROM , I ) 

DO  5  J=1P1,7 

5  7(IT0,JI=T{ IFROM, J ) +RM*P I V(U ) 

110=110-1 

6  IFRGM=IFR0M-1 

IF ( IFROH.GE. I )  GO  TO  4 
DO  7  J=IP1,7 

7  T ( I , J )=P IV (J } 

DO  8  1=1,5 
J=7-l 

K=6— I 

DO  8  L=J ,6 

8  T (K,7 )=T C K ,  .  T(K,L)*T(L,7) 

DU  9  1=1,6 

9  V(I)=T(I,7) 

RETURN 

END 
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AFAPL-TR-67-125 
Part  II 


SIBFTC  INPUT  0ECK.M9 4/2, XR7 

SUBROUTINE  INPUT 
COMMON  /PUJNT/IUATPT 
COMMON  /  ALL/ 

1HORD  t IDES  .JOES  . KDES  .MODE  .1NIT  .1DUMP  » IAMTP 
2IGASMX.IDBURN.lAFTbN.IOCD  ,IMCD  . I D5H0C , I MShUC , N02  FL T 
31TRYS  .LOOPER.NOHAP  , NUMMAP, MAPEOG, TOLALL .ERR C6 J 
COMMON  /DESIGN/ 

1PCNFGU.PCNCGU.T4GU  . DUMD1  , DU«D2  .DELFG  .OELrN  .OfcLSFC 
2ZF0S  .PCNFDS.PRFOS  .ETAFDS. WAFOS  .PRFCF  .ETAFCF .WAFCF 
3ZC0S  .PCNCOS .PKCDS  , ETACDS, WACDS  .PRCCF  .ETACCF ,MAf CF 
4T4DS  .WF8DS  ,DTC'JDS,ETABDS,WA3C0S.0PC0DS,0TC0CF,FTABCF 
5TFHP0S. CNHPOS ,E  THPDS. TFHPCF.CNHPCF , EThPCF .DHHPCF .T2DS 
\  6TFLPDS.CNLPDS,£TLPDS.TFLPCr,CNLPCF,ETLPCF,0HLPCF,T210S 

7T24DS  .WFODS  ,0TDUDS,ETAD0S,ttA23DS»DP0UDS,CTDUCF»ETADCF 
1  8T7DS  .WF ADS  ,DTAFDS,ETAADS,MG(.CDS«CFAFDS,DTAFCF,ETAACF 


9A55 

.A25 

•  £6 

i  A7 

,A8 

,A9 

,A28 

»A29  , 

APS55 

COMMON 

,AH55  .CVDNOZ 
/  FRONT/ 

.CVMN02 

.A8SAV 

. A95AV 

.A28SAV 

. A29SAV 

IT1 

.PI 

.HI 

,S1 

.T2 

»P2 

»H2 

,S2 

2T21 

,P21 

,H2i 

.  S21 

♦  T3 

,P3 

,H3 

,S3 

3T4 

,P4 

.H4 

»S4 

»T5 

»  P5 

,H5 

»S5 

4T55 

.  PS5 

,H55 

.555 

,BLF 

.BLC 

.BLDO 

.BLOB  , 

5CNF 

,PRF 

,ETAF 

.HAFC 

*WAF 

,KA3 

»WG4 

,  FAR4  , 

6CNC 

.PRC 

,ETAC 

.  WACC 

•  WAC 

,ETAB 

.OPCOM 

.DUMP  , 

7CNHP 

.ETATHP 

.DHTCHP 

,OHTC 

.SLHP 

,WG5 

.FAR  5 

,CS  , 

8CNLP 

.ETATLP 

.DNTCL? 

.OHTF 

.BI.cP 

.W655 

.FAPS5 

.HPEXT  , 

9AM 

,altp 

,ZF 

,PCNF 

,ZC 

,PCNC 

,WFB  . 

ATFFH? 

COMMON 

.TFFLP  .PCSLF 
/  SIDE/ 

.PCBLC 

♦ PCBLDU. PCBLIB 

,PC.  -HP 

.PCBLLP 

XXP1 

•  XL  AF 

.XWfiC 

»XBLF 

, XBLOU 

,XK3 

.0UMS1 

,DU*S2  , 

XXT21 

,XP21 

.XH21 

,XS21 

»T23 

,PZ3 

,H23 

*SZ3 

3T24 

,P24 

-H24 

,524 

»TZ5 

.P25 

,H25 

»S25 

4T2B 

,P28 

,H28 

,S28 

,T29 

,P2  9 

»H29 

»S29 

5HAD 

,MFD 

t 

. FAR24 

.  ETAD 

.DPDUC 

.BYPASS 

,DU«S3 

6TS28 

COMMON 

,PS28  ,V2S 
/  BACK/ 

,AM28 

,TS29 

,PS29 

,  V29 

,AM29 

XXT55 

,XP55 

.XH55 

,XS55 

,XT25 

,XP25 

,XH25 

»XS25 

XXKF3 

.XWG55 

.XFAR55 

i.XWFD 

.XWG24 

. XFAR2A 

,XXP1 

.DUMB 

3T6 

.P6 

,H6 

»S6 

,T7 

»P7 

.H7 

.57 

4T8 

,P8 

.H8 

.SB 

»T9 

»P9 

,H9 

,S9 

5KG6 

.MFA 

•  WG7 

,FAR7 

.ETAA 

.DRAFT 

♦  V55 

,925 

6PS6 

»  V6 

.  AH6 

.  TS7 

»PS7 

>V7 

.AM7 

,AM25 

7TS3 

,PSB 

.  V3 

,AM8 

»TS9 

,P59 

.V9 

»AM9 

3VA 

,fro 

t  VJD 

.  FGHD 

,VJM 

,FG«M 

,F6P0 

,FS?M 

9FGM 

t?GP 

,WF? 

,WGT 

.FART 

»FG 

»FN 

,SFC 

NAMELIST  /DATAIN/ 

AITITLE .  IDES  .MODE  ,INIT  ,100'MP  .IAMTP  .IGASMX, 

BIDBUSN  lAFTBN.IDCD  ,  IMCD  .NC2FL7, ITRY5  .TOLALL. 

CZrliS  .PCNFDS.PRFOS  .ETAFDS, WAFDS  .PRFCF  .ETAFCF, WAFCF  , 
D2CDS  ,PCNC OS.PRCOS  .ETACDS. WACDS  .PRCCF  .ETACCF .WACCF  , 
ET40S  .HFBDSr  . DTC0D5 .ETABDS, WA3CDS.DPC0DS .DTCGCF .ETABCF, 
FTFHPDS . CNHPDS ,E  THPDS.TFHPCF, CNHPCF , ETHPCF .DHHPCF ,T20S  , 

GTFLPDS.CNLPDSjETLPDS .TFLPCF ,CNLPCF,£TLPCF»DHLPCF,T2iDS  , 


HT240S 

.KFODS 

.D70UDS.ETADDS 

.  WA23DS  ,  CROUDS  .DTDDCF 

.ETAOCF, 

IT7DS 

.WFADS 

.dtafos.etaaos 

,W&6Cfe,0PAF£)S,DTAFCF 

.ETAACF, 

JA55 

,  A25 

,A6  ,A? 

,A8  »A9  ,A20 

,  A29  , 

KT2 

»P2 

,74  ,Z? 

,PCNF  ,ZC  »PCNC 

» 

LHFB 

.TFFhP 

.TFFLP  ,AH 

,ALTP  .ETAS  .hPEXT 

» 

hpcblf 

.BLF 

.PCBLC  ,SLC 

.pcbldu.sldu  .pcblos 

.BLOB  . 

NPCBLHP 

,BLHP 

, PCBLLP, 3LLP 

» 
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liS^rsSgj 

WS§ 


Ir^s-;-^ 
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I  :C -■  -' 


i  *£&> 


|?Vi2* 


MM 


ISi 


liZZ-t 


vmssz 


» -<  :v£ 


h>; 


l-yi';.: 


|^*7v* 


S*  V  <%  ■ 


-.'-Sv^ 


i#  - 


■«.-! 
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AFAPL-TB-6'’-125 

P&Tltl 


OPS55  ,AN5B  ,724  «ETAD  ,WFD  , 

P77  ,£TAA  ,  UFA  , CVDNOZ , CVMNOZ , DELPS  *L>rLFN  ,OELSFC 

C  ***  1TITLE=1  HILL  READ  IN  TITLE 
C  ***  JOES  =1  FOR  CALCULATING  DESIGN  POINT 
C  **=  NODE  =0  FOR  CONSTANT  T4 

C  ***  NODE  =  1  FOR  CONSTANT  PCNC 

C  r**  MODE  =2  FOR  CONSTANT  WPS 

C  #**  IN  IT  =1  WILL  NOT  INITIALIZE  POINT 

C  «**  IDUMP  =1  WILL  DUMP  LOOPING  WRITE-OUTS  If  ERROR  OCCURS 

C  ***  IOUKp  =2  v:ILL  DUMP  LOOPING  WRITE-GUTS  AFTER  EVERY  POINT 

C  ***  IAMTP  =0  WILL  USE  INPUT  AH  AND  MIL  SPEC  £TAR 
C  I AH TP  =1  Hilt  USE  INPUT  AH  AND  INPUT  ETAH 
C  **«=  IAHTP  =2  WILL  USE  7  2  AS  71=T1  +  T2  AND  STANDARD  PI 
C  -**  IAHTP  =3  KILL  USE  P2  AND  STANDARD  TI 
€  ss*  J AMT P  =4  WILL  USE  T2  AND  P2 
C  !GASKX=-i  SEPARATE  FLOW,  INPUT  A6 
c  ***  IGA$HX=0  SEPARATE  FLOW,  A6=A55 

C  1GASMX=1  WILL  MIX  DUCT  AND  MAIN  STREAMS,  A6=A2S+AS5> 

C  -*#  IGASHX-2  HILL  HIX  OOC7  AND  MAIN  STREAMS  *  INPUT  A6 
C  I0BURN=1  FUR  DUCT  BURNING,  INPUT  724 
C  IDBURN-2  FOR  DUCT  BURNING,  INPUT  *FD 

C  iAFT5N=i  FOR  AFTERBURNING,  INPUT  T7 

C  *«■  I ArTBN=2  C0R  AFTERBURNING,  INPUT  WFA 

C  »**  IDCO  =1  DUCT  HOZZLr  WILL  3E  C-D 

C  ***  IH£D  =1  MAIN  NOZZLE  WILL  BE  C-D 

C  ***  kCZFlT=1  FOR  FLOATING  MAIN  NOZZLE 

C  ***  «UZFLT*2  FUR  FLOATING  DUCT  NOZZLE 

C  *«  NUZFLT =3  FOR  FLOATING  MAIN  AND  DUCT  N0Z1LFS 

C  «?s*  I TRVS  =N  NUMBER  OF  PASSES  THRO  ENGINE  BEFORE  QUITTING 

DIMENSION  TITLc(12) 

DATA  AWURD/6H  INPUT/ 

IOA7P7=iDA7PT*l 
PRINT  100,  I i*A TP  r 
CALL  ZERO 
WURD=AWCRD 
READ! 5  tDATAlN J 

IF 4  !  AF7BN.GT«.0*0R»iD3UR5i»&7  .01  lNiT=l 
IFtniTtF.EO.IJ  READ!  b ,  101 J  TiTLE 
’T III  fc-0 

WRITE  1 6, 1025  TITLE 

IFiMnUE.EO.O)  *RITE48»103)  IOEStAM, ALTP,T4  ,T24,77 
IF (MODE. £0.1)  WRITE f 5  » 104)  IDES, AH, ALTP,PCNC  ,T24,T7 
IFIMG0E.EQ.2)  WRITE{B,105)  IDES, AS,ALTP,WF 3  ,T24,T7 
CALL  CUINLT 
RETURN 

100  FORMAT  4 1 1H  DATA  POINT, 13 J 
id  FORMAT  < 12A6 f 

102  FORMAT { iHl ,30X1246}  . 

103  FORMAT «  IHG»7h  IDES- ,  1 3  - 1  0X7k  AH=,P7.3»6X7M  ALT«>-,P7.0, 

16X7H  T — * , FB . 2 ,  SX7H  724=  TF8. 2,5X70  ‘  T  7=,F8 .2  ,6HSS>pJSS ) 

104  FURHAT { 1H0,7H  IOES=,I3.10X7H  AM=, F7.3 ,6X7H  ALTP=,F7#0, 

16X7H  PCNC=,F8.3,SX7H  724=, FS. 2,&X?H  T7=,F8.2 »6H$SS*SS ) 

105  FORMAT  1 1M0»7H  10ES=  ,13, 10X7H  AM=,F7.3,6X?H  C.L7P=,F7.G, 

16X7H  WFB=,F8.4.bX7H  T24=,F8.2,5X7 H  T7=, F8 .2 ,6H>?ASSS > 

END 
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Part  II 


SloFTC  ZERO  0ECk,h9A/2,XH? 
SUBROUTINE  ZEkU 
COMMON  /  ALL/ 


49ES=C 

IMI7=0 

I03UR.N=0 

IfiFTBK=0 

iD3WK:=3 

IHSHOL=3 

T2C=T2 

P29=f >2 

7r.a*74 

DO  •  I *i ,63 

Zl(I»=0. 

00  ?  1=1, AS 
22(3 !=0- 
90  3  I=i,72 

Z3m=c» 

T2-*T20 
?2=P2C 
79-7*5 
CALL  SY6 1 1 » 
RETURN 


1W0R0 

,  IDcS 

» JOES 

,xOES 

.MODE 

,  IMT 

, I DUMP 

,!AM7P 

1 

i. 

2IGASMX 

3I7RYS 

COHHOn 

, !  D8U8N ,  1 AFT8N 
,LGO?ER, NOMAP 
/  FRONT/ 

,IOCO 

,f«UMMAP 

,  IMCD 
.HAPfcOG 

• *DSHOC 
, TDi ALL 

,IMSKOC 
,EKR (6 ) 

.NU2FLT 

f 

i 

.  > 

n  K 

'•  ( 

•» 

1T1 

,Pi 

,H1 

,Si 

,72 

»P2 

,rt2 

,S2 

t 

£ 

2T?1 

»P21 

,H21 

t  S21 

,12 

,P3 

,h3 

,53 

V  V 

3TA 

,P* 

,SA 

,75 

*P5 

,H5 

,55 

« 

A75S 

,P55 

,H55 

,555 

,  8LF 

tSLC 

•  BLDU 

{BLV8 

t 

.  ? 

5CNF 

,PRF 

iETAF 

,WAFC 

,WAt- 

»  K«3 

,5G* 

-PARA 

t 

v  Y 

6CNC 

,PnC 

,ETAC 

,  WACC 

,NAC 

,  fclAo 

.opcrci 

,50Hf 

f 

7CNHP 

,c7ATHF 

,DHTCHP 

»DHTC 

,BLHP 

,t»G5 

,FAR5 

,CS 

9 

: 

’’i 

BtUl? 

,E  T3TLP 

»9HTCLP 

t  0H7F 

,ULLP 

,  WG55 

,FAR35 

,H?fcXT 

9 

■ 

,  if 

9AM 

,AJ_T» 

,ETAR 

,2F 

,?CNF 

,zc 

,PCNC 

,  WFS 

9 

: 

3 

ATFFflP 

COMMON 

fTFFLP  ,PC2LF 
/  SIDE/ 

, PC3LC 

,PCS«  Dll 

,PC3L03 

,RCBLH? 

,PC8LLP 

f 

..." » 

XXPi 

,XH5F 

tXWAC 

,X3L? 

, X8LDU 

,  XH3 

,90MS1 

,C0MS2 

9 

t 

XXT2? 

,x?21 

,XK2I 

,XS2I 

.723 

,P23 

,H23 

,523 

9 

<r 

3T2A 

,P2A 

•  H  2  A 

♦  S2A 

,725 

,P25 

,M25 

,S25 

* 

4T2L 

,P28 

,H28 

,S2S 

,T29 

,P2  9 

,h2? 

•  529 

9 

* 

5WAD 

,3*F0 

,WG2* 

,?AR29 

#ETAO 

»DFDUC 

.CYPASS 

,G&hS3 

9 

* 

6TS28 

CGKMON 

,PS28  ,V28 

/  5ACX/ 

,AH28 

r  TS2V 

.PS2V 

»  V29 

,  AM29 

P 

XXT55 

,  XP55 

»/f»55 

,XS55 

,XT?5 

,XP25 

,XH25 

-XS25 

9 

: 

i 

XXWF8 

.XMG55 

.XFAR55 

fXHFO 

, XM62* 

,XFAR2A 

,XXPI 

,OOMB 

? 

376 

,P6 

,H6 

»S6 

,77 

,P7 

,H7 

,57 

9 

ATS 

,PM 

#H8 

,S6 

,79 

»P9 

,H9 

»S9 

9 

k 

'■  3_ 

5HG6 

,RFA 

»WG7 

,F*R7 

,E7AP 

-OPAFT 

,V55 

»  VZ3 

9 

.1 

6=S6 

,V6 

,AM6 

»  TS? 

,PS7  - 

,V  7 

»AN7 

,ab25 

f 

^  _  F- 

7TSM 

,PSS 

,va 

,AH8 

,TS9 

»P$9 

,V9 

-AH9 

SVA 

,F  P.9 

,  VJD 

,  FGHO 

,yjM 

,f  GMM 

,PGP9 

,FGPK 

# 

L- 

1; 

9FGM  *FGP  ,wFT  ,kST  ,FART 

DIMENSION  211635, Z2IASJ.Z3772I 
ECUIVALENCE  tZl ,T1 5 , IZ2, XP1 > , (Z3,X 
IOES=0 

755  5 

,FN 

,.CFC 

E- 

l 

t 

•  } 

11 

s-2 

t  s 

■  a 

.u 

j  >< 

VI 


-3 


,-rs 


r  ■ 

r  -  rr: 


3 


fefe% 
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1st  g&rp  r 


rasa 
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esr-  r  IS 

6-51 

tL~'  ‘-'. 


\  -t£S 


;  >> 
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S1BF7C  CQlWlT  0£C?;,K94/2,XR7 
SUBROUTINE  COJiiLT 
Ctf*KON  /  ALL/ 

i'eO.RO  ,JDES  ,JDtS  ,KOES  .KOOfc  ,ifclT  .IDlSHP  ,JAKTP 
2!&ASH3fjiDStmft»IAfYBH.  10C0  ,;KCD  ,  iOSHOC  ,lMSHUC,»iQ2FLT 
3ITRYS  ,LCO?SP.,f#a*5flP  ,NUHHA?  9fSAP£U6,T0LALL  ,£RR  C6  ) 

COMBOS  /DESIGN/ 

i.«lNr«J,PCf:esUfT4G0  -DUHOl  ,DU.R©2  . DtLFG  ,DELFN  .DEL3FC 
22FOS  jFCN'FOS.P&FOS  »ETAF£>S»«iAFOS  t?P.PCF  j£TAFCF,WAFCF 
32C05  sPCfiCOStPRCOS  ,£TACCS,*ACD5  ,P?-CCF  ,ETA(XF,KACCP 
4T40S  » is? 303  jOTCODS,  £7ASDS»R.43COS,OPCOOS ,OTCOCF,ETABCF 

5trKFCS»CNr«PDSjETKPOSf7FHFC?»CNHPCF,ETHPCFiOHHPCF,T2DS 
5TFLF0S»C*iLPDS,E7LPDS;TFLPCF»CMLPCF,ETLPCF,DHLPCF,T21DS 
7T24CS  t«?50S  ¥G70UDSjE7A00S»WA23DS»CPDU0S  »DTOOCFt  ETADCF 


6770S 

, a FADS 

,dTA?OS 

tETAADS 

,t£6CDS 

.OPAFDS 

,f»7AFCF,ETAACr 

9A55 

f  A25 

*A6 

*A7 

,A8 

*A9 

,A2B  ,A29 

APS  55 

,AH5S 

sCYBNGZ 

tCVKNQZ 

,  A&SAV 

,A9SAV 

, A28SAY »A29SA\ 

CORBOB 

/  FRONT/ 

1T1 

»P1 

i»‘s 

,Si 

»T2 

,P2 

,H2  »S2 

2T21 

♦  ?2i 

,H21 

,S21 

.73 

,P3 

r*2  *SA 

37* 

,?4 

»ka 

»S4 

?T5 

,PS 

,H5  ,S5 

4T55 

,P55 

,H55 

,S55 

,8LF 

,BLC 

,BLOU  .BLOB 

5CHF 

,PKF 

•ETAF 

,  XAFC 

,WAF 

,BA3 

,  HG'i  ,  FAR4 

6CNC 

,PRC 

♦fcTAC 

,WACC 

,WAC 

,  cTAB 

,OPCOM  .{JUMP 

7CNHP 

,ETATHP 

sDBTCHP 

,OHTC 

,3LHP 

,WG5 

,FAR5  ,CS 

8CMLP 

,ETATLP 

,OHTCLP 

,DHTF 

tBLL? 

,KG55 

,FAR55  ,HPEXT 

9AM 

,ALTF 

,STaR 

tZF 

,PCNF 

,ZC 

,PChC  ,MF8 

ATFFHP 

»tfflp 

•PCBLF 

tPCBLC 

,PC6L£XJ 

i.PCBLUS 

,PC3LHP»PCBLLP 

3ATA  AHORO/6HCQ1HLT/ 

WORO-cAVOHO 
AJ*773.26 
G=32. 174049 

ALT=*LTP*2.08555315+07/{2.0S55531E+07-ALTP) 

CALL  AT*nSUL7,Tl,XXl,XX2,XX3,Pi»CS»XX4»mKl 
IF5IAMTP.E0.2}  T1*T1+T2 
IFliAMTP.BE. II  CALL  RAM {AH, £7 AR I 
FAR«0.0 

CALL  PROCOMIFAS,Tl,C$,XX2,XX2,RI,PKil,Hl) 

S1«PHI1-R1*AL0&(PI) 

H2=Hi+{AM»CSi*»2/l2.*A,J*G» 

P2T=1. 

OO  2  1=1, 10 

CALL  THERKOeP2T5H2,T27,52T,AW, 0,0.0, II 
IF1A35IS27— SU.Lc.Ot0001=S1  ?  GO  TO  3 

P27=P1*EXPC I A*/l. 956375)  ={ {S2T-S1 )T( 1.986375/AH} *ALOGC P2T/P1 1 1 . 

CALL  ERROR 

RETURN 

IFIIAHTP.EO.3.0R.IAHTP.E0.4I  ETAR-P2/P2T 
P2=S7AR«P27 

IFUAMTP.BE.4l  CALL  THcRMlM  P2»H2,72«52  »XX5, 0,0.0,  1 1 
IFfi.'stTP  £0.4}  CALL  THERMO  I P2  ,  H2*  72*  S2*XX5,0  ,0.0,0) 
IFnHIT.EO.l)  SO  TO  6 
IF ( IOES.cO.il  GO  TO  4 

PCKF=«iESS{H0l>€>T4rTA0S.PCf«C,PCAfC0S,WFB.HFBDS,72,T2DS,PCNF0S) 

PCHFGO=FCrtF 

60  TO  5 

PCBFaFCHFOS 

pchfgu«pc»*f 

T205=T2 

2F*2F0S 

RETURN  * 
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SlHFTC  ATM.062  I1ECK,H94/2,XR7 

SUBROUTINE  ATMOS  (ZFT,  TM,  SIGMA,  RHO,  THETA »  DELTA,  CA,  AMU,  K) 
THIS  IS  A  SUBROUTINE  TO  COMPUTE  CERTAIN  ELEMENTS  UP  THE  1962 
U«S.  STANDARD  ATMOSPHERE  UP  TO  90  KILOMETERS. 

CALLING  SEOUENCE... 

CALL  ATHOS  {ZFT ,  TM,  SIGMA,  RHO,  THETA,  DELTA,  CA,  AMU,  XI 
Z FT  *  GEOMETRIC  ALTITUDE  (FEET) 

TM  =  MOLECULAR  SCALE  TEMPERATURE  (DEGREES  RANK1NE) 

SIGMA  =  RATIO  OF  DENSITY  TO  THAT  AT  SEA  LEVEL 

MHO  =  DENSITY  (LB-SEC**2-FT*»(-4)  OR  SLUGS-FT**3) 

THETA  =  RATIO  OF  TEMPERATURE  TO  THAT  AT  SEA  LEVEL 

DELTA  =  RATIO  OF  PRESSURE  TO  THAT  AT  SEA  LEVEL 

CA  *  SPEED  OF  SOUND  ( FT/SEC  I 

AKu  =  VISCOSITY  COEFFICIENT  (LB-SEC/FT «*2) 

X  a  1  NORMAL 

*  2  ALTITUDE  LESS  THAN  -5000  METERS  UR  GREATER  THAN  90  KH 
=  3  FLOATING  POINT  OVERFLOW 

ALL  DATA  AND  FUNDAMENTAL  CONSTANTS  ARE  IN  THE  METRIC  SYSTEM  AS 
THESE  QUANTITIES  ARE  DEFINED  AS  EXACT  IN  THIS  SYSTEM. 

THE  RAUIUS  OP  THc  EARTH  (REFT 59*  IS  ThE  VALUE  ASSOCIATED  WITH  THE 
1959  AROC  ATMOSPHERE  SO  THAT  PROGRAMS  CURRENTLY  USING  THE  LIBRARY 
ROUTINE  URL  NUT  KEUUIRE  ALTERATION  TO  USE  THIS  ROUTINE. 

DIMENSION  HB' 101 »Tm3( 10 I, DEL TAB ( 10) ,ALM( 101 
DATA (HBJ 1 ) ,  TMBllJ,  DEL  TAB ( I ) ,  ALK{ I) , 1 =1 , 10 ) / 

A  -3.0,  320.65,  1.753&3E  00,  -6.3, 

B  0.0,  283.15,  l.OUOOOE  00,  -6.3, 


C 

11.0, 

216.65, 

2.23361E-01, 

O.C, 

0 

20.0, 

216.63, 

5.40328E-02, 

1.0, 

E 

32.0, 

228.65, 

8.56663E-03, 

2.8, 

F 

47.0, 

27C.65, 

1.09455E-03, 

0.0, 

G 

52.0, 

270.65, 

5.622B9E— 04, 

-2.0. 

H 

dl.O, 

232.63, 

1.7971 8E-04, 

—4.0  , 

I 

79.0, 

1BU.65, 

1.024)  6-05, 

0.0, 

J 

38.743, 

180.65, 

1.6223  E— 06, 

0.0/ 

DATA  R6FT59/2.0B5S531E  07/,  GZ  /9. 

80665/ 

A 

AmZ 

✓28.9b44 

it  RSTAR  / 8. 

31432/ 

B  FTT{IKM/*.o4BE-04  /,  S  /110.4  /, 

C  AMUZ  /1.LD246-03  /,  CAZ  /1I16.45/, 

P  RHD2  /O. 0764 74  /,  GZENG  /32.174I/ 

C  CONVERT  GEOMETRIC  ALTITUOE  TO  GEOPOTEHTIAL  ALTITUDE 
HFT  =  IKEFT59/(REFT59*ZFTJ J»7FT 
C  CONVERT  hrT  AND  ZFT  TO  KIUJMETFKS 
Z  =  FTTUK«*ZFT 
H  =  FTTOXM*HFT 
K  =  1 

T«7  =  T«H(2) 

IF  (H.LT.— 3.0.0K.Z.GT .90.0)  GO  TO  16 
DU  10  M=l,10 
IF  (H-weiMlI  11,12,10 

10  CUNT  IM.IE 
GO  TO  16 

11  M  z  M— l 

12  DFLM  =  H-HBIMJ 

IF  IALMIMJ.FCJ.0.0)  GO  TO  13 
TMX  s  TM3(M)+ALrt(M)»0FLH 

Z  GRADIENT  IS  NON  ZERO,  PAGE  10,  EQUATION  1. 2.10-0) 

DELTA  *  0€LTAB(M)s;(TMB(M)/TPX)«*(GZ«AMZ/IRSTAR«ALR(M)))} 


V 
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13  TKK  = 
GRADIENT 

DELTA 

14  THETA 
SIGNA 
ALPHA 

CONVERSION 


15 

16 
17 


uu  III  14 

TrtBiMi 

IS  ZERO,  PAGE  10,  E0UA7I0N  I.2.10-C4) 

=  DELTAB  f  M) «£XP { -GZ*AHZ-DELH/ { RSTAR^ThB (M) J ) 
=  THK/ThZ 
*  D6LTA/THE7A 

SORT ( THETA--3) * { { TNZ+S ) / ( TMK+S 1 ) 

TO  ENGLISH  UNITS 
TM  =  1.8*THK 
RHO  =  RHOZ*SJGHA/GZENG 
CA  =  CAZ^SORTl THETA J 
AMU  *  AKUZ«ALPHA/GZENG 
CALL  5VERFUJ) 

GO  TO  (15,17),  J 
K  =  K+2 
GO  TO  17 
K  =  2 
RETURN 
END 


i 

} 


19 


AFAPL-TK-67-12S 
Part  n 


4I8FTC  RAH  D£CK»M94/2,XR7 
SUBROUTINE  RAM{AM,E7AR) 
IF(AM.GT.l.)  60  TO  2 
ETAR=1. 

1  RETURN 

2  IF(A*.GT.5.)  GO  TO  3 
ETAR=K-0.075*(  ( AM-I.  ) **1 .35  J 
GO  TO  1 

3  ETAR=800./UAH**4)+935.) 

GO  TO  1 

END 


SO 
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SlBFTC  GU&SS  OfcCK,hVA/2tXk7 

FUNCTION  GUESS{h,T,TO,P,PD,W,HD,D»DD,VO) 
IF(H.EO.O)  GUESS=VO*( { T/TD) **1 .60) * t (DD/l»**0.50) 
IF(M.EO.l)  GUESS=VD*UP/‘PD)**1.80)*(  (D0/D)**0.33) 
IFtH.E0.2I  GUESS=V0*{{W/KD)**0.33»*{(UD/D)**1.00) 
IF(H.E0.3J  GUESS=VO*(  { W/WD)**0. 00 )*({ P/PD) **0.50} 
IF (H.E0.4  J  GUESS=VD*((W/WD) **0.001*1 (P/PD)**G.50) 
IFIM.E0.5)  GU£SS-VD*( { T/TO 1**1. 10 )*( (DD/D)**0.60) 
IF IM.E0.6 1  GUESS=VD*( (P/PD) **1 .00) * ( I D/DD)**0.25) 
RETURN 
END 
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SI3FTC  COFAN  DECK,M9A/2,XR7 
SUBROUTINE  COFAN 
COMMON  /  ALL/ 

1K0RD  ,IOES  « JOES  ,XDES  ,NOOE  ,INIT  ,IDUMP  , IAMTP 
2IGASHX , IDBURN, I AFTBN, IOCO  ,IMCU  . IDShOC.IHSHUC.NUZFLT 
3JTRTS  tLOOPER, NOMAP  .NUMMAP.MAPEDG, TPLALL ,ERR C6 1 


1WORO 

21GASMX 

3ITRTS 

COMMON 

1PCNFGU 

22F0S 

3ZC0S 

AT4DS 

5TFHPOS 

6TFLPDS 


/DESIGN/ 

PCNCGU,T4GU  ,DUMD1  ,t»UH02  ,OELFG  ,UELFN  .DfcLSFC 
PCNFOSrPRFDS  ,ETAFDS,HAFOS  ,PRFCF  f ETAFCF »kafcf 
PCNCDS,PRCOS  ,£TACDS,VACDS  ,P RCCF  ,ETACCF,*ACCF 
WFBOS  tOTCUDSt ETASDSt WA3COS tOPCODS .OTCOCF  *ET  A8CF 
CNHPOS,ETHPOS ,TFHPCF .CNHPCF , ETHPCF .DHHPCF ,T2US 
CNLPDS.ETLPDS, TFLPCF,CNL?CF,ETLPCF  .DHLPCF ,T21PS 


7T2A0S 

,MFOOS 

,otduos,etaoos 

.VA230S 

,OPOUOS 

.DTOUCF.ETAOCF, 

8T70S 

,WFADS 

.otafds.etaaos 

, WG6C0S 

jDPAFDS 

,dtafcf,etaacf. 

9A55 

,A25 

• A6  ,A7 

,A8 

t  A9 

,A28  »&29  , 

APS5S 

,  AH55 

,CV0K07,CVMN0Z 

,  A8S.  V 

, A9SAV 

, A28SAV .A29SAV 

COH**CN 

I  /  FRONT/ 

1T1 

»P1 

,H1  ,S1 

,T2 

.P2 

,H2  *S2  , 

2T2I 

,P21 

,H21  »S21 

tT3 

,P3 

,H3  ,S3  , 

3T4 

,P4 

.HA  .S4 

,T5 

,P5 

,H5  ,S3  , 

4755 

,?5S 

,H55  ,555 

,BLF 

,BLC 

,BLDU  .BLOB  , 

5CNF 

,PRF 

,ETAF  , WAFC 

,  RAF 

,  WA3 

,WG4  ,FAR4  , 

6CNC 

,PRC 

,ETAC  ,«ACC 

,WAC 

,  ETAB 

,OPCOM  ,DUMF  , 

7CNHP 

.E7ATHP 

,DHTCHP,OHTC 

,BLHP 

,WG5 

,  r  AR5  ,CS  , 

8CNLP 

,ETATLP 

.DhTCLP.OHTF 

,BLLP 

,HG55 

,FAK55  ,HPEXT  , 

NAM 

,ALTP 

,ETAK  ,ZF 

,PCNF 

,ZC 

,PCNC  ,WF8  , 

ATFFHP 

,tfflp 

,PC3LF  ,PCBLC 

,  PC8LIX1 

, PC8L03 

,pcblhp,pcbllp 

COMMON  /  FAN/CNX(15),PRXC15,151fWACXC15,ISl,ETAXClS,15}, 
1NCN,NPTC15» 

DIMENSION  WLHC21 

DATA  AHOKDf WLH/6H  COFAN ,6H  1101  ,6H  (Mil  / 

WORD»AWORD 

THETA=SORTC 72/518.6661 
CNF=PChF/CIQ0.*THE7Al 
1FCZF.LT.0.1  ZF=0. 

IFIZF.GT.1.1  ZF=|* 

CNFS^CNF 

CALL  SEARCH (ZF, CNF, PRF.MAFCtETAF, 
1CHXC11,NCH,PRXC1,11,RACXC1,11,ETAXC1,11,NPTC1  1,15.15,100) 

IF ( CCNF-CNFS1 .GT.0.0005«CNF)  MAPEDG=1 

IF  C IGO.EO.l .OR.IGO.EO.2 1  HRITE(8,1000;  CNFS.WLHC IG01 

FORMAT  C 19H0*  «  *  CNF  OFF  hAP,F10.A,2XA6,HH>  «  •SSSSSS) 

KAF=HAFC*P2/ THETA 

IFCIDES.NE.il  GO  70  1 

PRFCF»CPRFDS-1 . l/CPRF— i. 1 

ETAFCF=ETAFDS/cTAF 

MAFCFsNAFOS/WAF 

WR I TE { 6 , 100 1 PRFCF , ETA  FCF , KAFCF , T20S 

FORMAT CllHOFAN  DESIGN, 13X6H  PRFCP=,E15.8,8H  ETAFCF=,E15.8, 

18H  HAFCF=,E15.6,8H  T2DS=,F15.8) 

PRF«PRFCF*CPRF-1.)+1. 

ETAF=ETAFCF*ETAF 

RAF»WAFCP«VAF' 

PCNF=1C0.*THETA«CNF 

DUM01*PCNF 

CALL  7HC0MP I PRF , ETAF, T2 , H2 , S 2 , P2 , T2 1 , H21 , S2 1 , P21 1 
1FCPCBLF.GT .0.1  BLF=PCBLF*WAF 
lFCJOES.EU.il  GO  TU  7 
JOES=l 

IFCINIT.EO.il  GO  TO  6 


}! 
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IF { I09S.E0.1 )  GU  TU  4 
IF (MODE. NE. 2)  GU  TO  2 

T4=GUESS{3tYl,Y2»PCNF,}>CNFUS,WFB,WFbr/S,Y7,YS,T40S  ) 
PCNC=GUESS(4tYl,Y2,PCNF,PCNFnS,WFe,WFB0S»Y7«Y8tPCNCUSl 
GO  TO  5 

IF (MODE.EO. 1 )  GO  TO  3 

PCNC=GU£SS(5,T4,T4DStY3,Y4t Y5, Y6, T21 , T21 OS ,PCNCDS ? 

GO  TO  5 

T4=GUESS(6tYl,Y2,PCNC,PCNCDS,Yf>,Y6,T21  tT210S,T40S) 

GO  TO  3 

PCNC=PCNCDS 

T4=T40S 

HFB=HFBOS 

T210S=T21 

ZC=2COS 

PCNCGU=PCNC 

T4GU=T4 

IN  I T=0 

CALL  COCO«P 

IF (NOMAP. EOc7)  PCNF=OUMDl 

RETURN 

END 
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*IBFFC  COCOMP  nECX,M9A/2.XR7 
SUBROUTINE  COCOMP 
COMMON  /  ALU 

2#ORD  *  IDES  .JOES  ,KOES  .MODE  ,IN1T  .IDUNP  .IAMTP 
2iGAS.*»X,IDBUkN,IAFTBN,!DCD  »IMCD  » IDSKJC . IMSKUC.NUZFLT 
3ITRYS  .LOOPER.KOKAP  .NUMHiP.MAPEUG.TOLALL ,ERR (fel 

common  /design/ 

I?CNFGU,PCNCGU»TAGU  .DUNOi  ,DUhD2  .DfcLFG  .UELFN  .OfcLSFC 
2XFDS  .PCNFOStPRFDS  .ETAFUS.WAFUS  .PRECF  .ETAFCFvWAFCF 
32COS  tPCNCOS »PRCDS  .ETACDS.WACDS  ,PRCCF  »ETACCF  ,HACCF 
4T4DS  .VF80S  ,DTCUDS»ETA8DS.WA3COS,UPCUI>StDTCOCP,E7ABCP 
STFHPDSiCNHPOS.ETHPDS. TFHPCF.CNHPCF.ETHPCF  .UHHPCF  ,T2I>S 
6TR.POS,C*LP0S,ETLPDS,TFLPCP,CNLPCFtfcTLPCF,DHLPCF.T2lO2> 
7T2*0S  tWFDDS  . D70UDS.  ETADDS»KA23DS,DPDUDS .OTDUCF  ^ETADCF 
877DS  tWFADS  .DTAFDS  »ETAAUS»W&6CDS. DPAFDS  .UTAFCF .fcTAACF 


9AS&  »A25  »A6  »A7  »A8  »A9  »A2B  »A29 

APS55  ,AH55  .CVDNOZ  .CVMNOZ  ,ARSAV  .ASSAY  .A2HSA V.A29SAV 
COMMON  /  FRONT/ 

1T1  ,P1  tHl  ,S1  »T2  »P2  ,H2  ,S2 

2721  ,?21  »H21  . S21  ,T3  .P3  .H3  .S3 

3T4  ,P*  ,H*  »S4  ,T5  ,P5  .H5  .S5 

AT5S  »PS5  »H55  »S55  .BLF  -BLC  .BLDU  .BLOB 

5CHF  »PP-F  .ETAF  rWAFC  ,WAF  .  WA3  ,WG4  .FAR4 

iCNC  .PRC  .ETAC  .WACC  .MAC  .ETAS  .DPCOM  .CUHF 

7CNHP  .ETATHP.DHTCHP.OHTC  ,BLHP  ,«G5  ,FAR5  ,CS 

8CNLP  .6TATLP.0HTCLP.0HTF  .BLLP  ,KG55  .FAR5S  .HPEXT 
9AM  .ALTP  .ETAK  «7F  .PTHF  .  tr  .orur  uld 


»PRF  .ETAF 
.PRC  .ETAC 
.ETATHP.DHTCF 
.6TATLP.0HTCI 
,altp  .etak 


.SI 

»T2 

»P2 

.H2  ,S2 

.  S21 

.T3 

»P3 

.H3  .S3 

,S4 

»T5 

»P5 

»H5  .55 

,S55 

,blf 

-BLC 

.BLDU  .BLOB 

,WFC 

,WAF 

.WA3 

,WG4  ,FAR4 

,  WACC 

.MAC 

.etab 

.DPCOM  ,CUMF 

,OHTC 

,BLHP 

,«G5 

.FAR5  ,CS 

.OHTF 

,BLLP 

*WG55 

.FAR55  .HPEXT 

,2F 

,PCHF 

.2C 

,PCNC  .WEB 

atffhp  .tfflp  .pcblf  , PCBLC  .PCBLOU.PCBLUB.PCBLHP.PCBLLP 
COHNON  /  CQMP/CNXf 15 J .PRX(15»15), WACXt 13.15) , ETAX ( IS, IS ) , 
1NCN.NPT ( 15  I 
DIMENSION  wl;h(2) 

OATA  AWQRO , WLH/6KCGC0RP , AH  fLO)  ,6H  (HI!  / 

W3RD*AWORO 

THfcTA»S0RTCT21/51B.668! 

CNC^PCNC/ < 100.«THETA ) 

IF  C2C.LT. 0. )  ZC=0. 

IFfZC.GT.l.l  ZC=1  • 

CNCS*CNC 

CALL  SEARCHtZC.CNC. PRC, HACC, ETAC, 

1CNX ( 1 } ,NCN .PRX (1.1) , WACXC 1.1) , ETAXC 1 ,1).NPT(1),15,15 , IGU) 
IFfHODE.EO.il  GO  TO  l 
IF(fCNC-CNCS).GT.0.0035*CNCI  HAPEDG=1 

1  IF  ( 1G0.E0.1  .OK.  IGtl.EO.2 1  WKITECB.iUCOI  CNCS.KLHl  IC-OI 
1000  FORMAT C19H0*  *  «  CNC  OFF  MAP, F10.4.2XA6, 1 1H»  •  •SAASAA) 

UAC*HACC<-P21/THETA 
lFfIOES.NE.il  GO  TO  2 
PRCCF=fPRCOS-l. 1/ JPRC-1 . ) 

ETACCF=£TACDS/ETAC 

MACCF=NACOS/HAC 

WRITE f6,100IPRCCF,ETACCF,WACCF,T2inS 
100  FORMAT UBHOCOMPRESSOR  DESIGN, 6X8H  PkCCE=,fcl‘5.8,BH  ETaCCF=,E15  B, 
18H  *ACCF=,E15.B»8H  T21CS=, F15. 8 1 

2  PRC*=PRCCF  »f  PRC-1. 1*1 . 

ETAC=ETACCF»ETAC 

NAC«=«ACCF»VAC 

CALL  ThCOM?  CPRC, ETAC, T2i ,H21 , S21 , P21 .73 ,H3 .S3 . P3 1 

IFfPCBLC.GT.O. I  BLC=PC8LC**AC 

KA3»*AC-BLC 

8LWJePCBLL*U*BLC 

8L0B«PCBL08»BLC 

BLHP*FC8LHPeBLC 


24 


A  FAPL-TR-57-125 
Part  n 


8LLP=PCBLLp-HLC 
IFjHOOb.NE.l  )  GO  TU  3 

IF(ABS{CNC-CWCS).Le.0.001«CMCS)  GO  TO  4 
WRITE  (8, 2000  JCK'CSiCNC 

2000  FORMAT { 10H0CNC  WAS=  »E15.8,1IW  A WO  ,E15.8* 

124H  CHECK  PCNC  iNPUTSSStJa? 

CALL  ERROR 

3  PCNC=100.*THETA*CNC 

4  CALL  C0C0M6 
RETURN 

END 


AFAFL-TR-8'.  -125 
F*rt  H 


-'-•Sllr 


>*  - "«.  .. 


-K\ 


J'l  - '  A 


il&f-TC  C0C0M8  OECX,K9*/2>XR7 
SUBROUTINE  CUCOMB 
COMMON  /  ALL/ 

1WORO  , JOES  •  JOES  ,*D£S  ,NOOE  ,1N1T  HOUR?  .IAMTp  , 
22GASMX  ,  IOBURN,  |  AETSM,  I QCL*  ,1HCD  ,  JOSHOC  ,  iMSBOC  .NOtELT, 
31TRYS  .LOOPt R.NUKAP  , NUKMAP, MAPEOG, TOLALL ,E»R CM 
COMMON  /DESIGN/ 

IPCNrGU.PCNCGU.TACU  .DUMOl  .OUHD2  ,DELE6  ,DELFN  ,DELSEC  • 
2ZFDE  .PCNFOa.PHFOS  .ETAEDS.MAFUS  .PKECE  ,£TAECF  .WuFCF  . 
3ZCO'  .PCnCDS.PXCDS  »£TACDS»yACD5  ,PRCCE  .ETACCE .WACCE  , 
ATAOS  jVFSO S  tCTCOOSjETASOS  ♦  V‘A3COS  jUPCOOS  »i)TCOCE  » 6TABCE, 
5TEHP0S  t  CHHPOS  .£ 7HPDS .  TFHPCf, CKHPCF  ,  ETMPLE  .DHHPCF  ,  T 20$  . 


1  i 

6TFLPOS»CNLP0S 

tETLPUS  tTf LPCF 

,C«LPCE,E7LPCE ,DHLPCf 

« T2 IOS  , 

l 

7T2ADS 

.WEDOS 

f DTOUDSt  E7A00S 

»  KA230S  tOPDODS .DTDUCE 

.ETADCf  » 

r 

t 

877DS 

»*<f  AOS 

•  OTAfOS  *£7AAyS 

.fc66CDS.DPAEQS,D7AECE 

.E7AACE, 

i  ; 

t 

\ 

* 

! 

J 

j 

9A55 

« A25 

,A6  ,A7 

,A8  ,A9  ,A2S 

,A29  , 

i  -  * 

AP-S55 

,AM5t> 

tCVuxQliCVMtGZ 

.ASSAY  , A95&V  .A2BSA Y 

, A29SAY 

5 

*  . 

COMMON 

/  FROM 

7/ 

L 

m 

.Pi 

,«l  ,S1 

s  T2  ,P2 

»S2  , 

£ 

f 

i 

2721 

«P21 

-,H2i  ,521 

,73  » P3  ,.«3 

.S3  » 

t 

\ 

37* 

.F  * 

,H*  ,5* 

»T5  , P5  ,H5 

>S5  , 

l 

*755 

,P55 

,H55  » S55 

,6LF  ,  SLC  ,  SLfXJ 

.8108  i 

:  : 

1 

\ 

5LNE 

,Pr.E 

,ETAf  ,MAEC 

,MAF  ,WA3  »MG* 

.FAR*  , 

t  . 

j 

ACNC 

•ar 

,-’*r  , uacc 

,MAC  ,£TAB  ,OPCO« 

.DOME  , 

\ 

i 

7CNHP 

. _ ' ATHP 

,L  ’  -P,DKTC 

-BLMP  ,Mf*5  »£ar5 

.CS 

t 

8CNLP 

.ETATLP 

,RH7(,LP,GH7r 

,BLLP  ,«555  .EAR55 

,H?£xf  , 

t 

9AM 

,ALTP 

»E  TAK  ,IF 

,PCNf  ,2C  ,PCNC 

,ME8  , 

ATFFHP 

( 7FE LP 

»PCBL?  t PC8LC 

, PC8L0U, PC3LL©  »PCBLHP 

.PC&LLP 

common  /  coMB/PsmsjjDfciTm.iSj.euu^.ipjfUPStUPTiii) 
DIMENSION  OlS|,DU«aOU5,I5l 
DATA  AWGH0/6Kti)CG*3/ 

WORO=AH£?RO 

Gi2i=0. 

0(31=0. 

P3PSl=iA.6V**P3 

WA3C=WA3*5'JETf73>/r3J>S.‘ 

IEHOE5.cO.lJ  VA3COi*MA3C 
DPCCJMs5EC00S« '  WA3C/W*  3CD3 ) 

IFCDPCOM.GT.l,  J  DECO-^l. 
p*ap*s; 1.~o°COh» 

IF  (7*.GT.3QOC«  )  T*=3000. 

IF fTA.GE.10C0. )  GO  TO  2 
T  .=1000. 

IF  INUOS.EO.n  fiAPfc05=l 
OTCO*7*-73 

IE  HOES. Vf-U  GO  TO  3 
OTC'JCE=OTCOOS/L’TCO 
U  i  COsOTCOC  F^tiTCU 
P3PS1N«P3PS! 

CALL  SEARCKI  — 1 .,P3P51N,07CQ» £  TAB. OUOMY , 

lPS!trifNPS>OcLm»i),c'AU,l).OUMBU«l.i»t«P7UI.lS,15,IGO) 

IF  I IG0.E0.7 )  CALL  ERROR 
iPIIOES.NE.1 J  GO  TO  5 
c  7  A  8CF  *E  TAfeOS/E 7*5 
ETA8*ETA8Cr*fc7A8 

HV»(H  i  l  (-.*59*31 7E- 19*7* }-.  203*1 1 AE-15  1*7**. 27836*36-11  }*T* 
i*. 2031501 E-07)oT*-. 2*53 1 ifcfc-03> *7*-. 9*332966-01 )*7*», 18*5557E*U5 
CALL  I Mf RKU (P* |HA  ,TA,XXi,XX2,0,0.0r0) 

EARA  =  (HA— M3  )  /  (NV*E7A5  J 
IF -FARA.LT. 0.  EAK*=0. 

RE8»=EAK*e»»A3 


|gS3 


AFAPL-TR-67-125 
Part  It 


b 

7 

ft 

9 


100 


If- (MJIDfc  .Ml- .21  t-'J  TO  K 
EHKW-twr-H-W^HXJ/WFB 


l)IR*S0KT<WFB/W»-8XI 
CALL  AcUUIR(0'l).T4 
GO  TO  (6,9.7  1.1  tit» 


FRHW.O. .20. .0.0001 »D1K ,T  AT » 1GU) 


T/,*»4T 


GO  TO  1 
CALL  ERROR 

WFB-WFBX  . . 

CALL  THFKMtJtP4,M4»T4*S4*XX2tl*F£H4*01 


U('A  ■UJ-  M  UA3 

IFflOES.EO.il  WRITE* 6, 100)  WA3C0S. FTABCP .UTCOCF 
FORMAT ( 17H0C0HBUS TOR  OFSI6N.7XBH  WA3CDS»»fcl 5.ft .ftH 
1BH  0TC0CF-.F15.B) 

CALL  COHPTB 


RETURN 

ENO 


FTABCF-.615.R 
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SiBFTC  COHPTB  0ECK,H94/2,XR7 
SUdkDUliKE  COHPTB 
COKttO X  /  ALL/ 

lKOKD  i  IDES  .JOES  ,£DES  ,HUOE  ,1NIT  »IOUHP  ,IAMTP  , 
2IGASHX.IDBURN, I AFT6N, IDCD  ,l*CO  ,iDSHGC,iHSHOC,HOZFLT, 
3ITRYS  ,LUGP£fc,NG*AP  ,  NU.RMAP,WAP£DG,TOLALL  .ERR  (6) 

CUHKOff  /DESIGN/ 

1PCNF6*J,PCNCGU,T4GU  .OUNDl  ,OUHD2  iDELFG  ,OELFN  ,Q£LS?C, 
2ZFDS  ,PCNFOS-?RFDS  .ETAFC-S.wAFDS  ,PRFCF  .ETAFCF.WAFCF  , 
3fC»S  .P'XCySfPRCOS  .ETACDS.yACOS  »PRCCf  ,ETACCF,WACCF  , 
4T4DS  i»»FeOS  *OTCGDS  tETABOS»MA3CDS»OPCOOS*DTCOCF*ETASCFi 
-  5TFhPOS,CKRPDS,£THPDS,TFHPCF,CNHPCF,ETHPCF,0KHPCF,T20S  , 
&TFtPOS»CfiLPOS  tcTLPOS  t TFLPCF  fCHLPCF ,  ETLPCF  .DHLPCF  »T2 IDS  , 
7T2V1S  ,KRK>S  ,OTSUOS»ETADOS,WA230S,DPOOOS,DTDUCF,6TAOCF. 
8T?OS  *KFAOS  t DTAFOS t ETAADS» WC6COS,DPAFl)S ,D7AFCF,cTAACF. 
9/.J5  ,A25  ,A6  ,A7  ,AS  ,A9  ,A23  ,A29  , 

A?355  ,AN55  ,CVDNOZ,CVHNQZ,A8SAV  ,A9SAV  ,A28SAV,A29$AV 

CONDON  /  FRONT 


Ifl 

,P1 

.HI 

,S1 

,T2 

,P2 

,H2 

,S2 

2T21 

,P21 

,H21 

,  S21 

,T3 

,P3 

»H3 

»S3 

3T* 

,P4 

*H4 

,S4 

,15 

,P5 

,H5 

,S5 

4T5S 

,P55 

,h55 

,  S55 

»3LF 

,BLC 

tBLOU 

,3LU3 

5CNF 

,PRF 

,etaf 

,WAFC 

fKAF 

,WA3 

,WG4 

,FAR4 

6CHC 

,PRC 

,ETac 

,  WACC 

,WAC 

,  ETAS 

,OPCON 

,DUHF 

7CNHP 

, ETATHP 

,  OH',  Cup 

,OHTC 

,3LHP 

,WG5 

,FAR5 

,CS 

8CNLP 

sETATLP 

.KiTCLP 

,OHTF 

,  BLLP 

,WG55 

,FAR55 

.HPEXT 

9AH 

,ALTP 

,cTAR 

,ZE 

»PCNr 

,ZC 

,PCNC 

,WFB 

ATFFHP 

,7Ff LP 

,PCBLF 

, PC8LC 

,?CBL0U,PC8LG8 

.PCBLHP 

.PC8LLP 

COMiON  /KTUS8/TFFX I 151  tZHXi  15,  131 ,UHTCX(  15,15  I ,E?ATX<  15,151, 
r  iMTFFSfMPTTFF 1151 

DATA  AWIJKO,WLO,WH1/6KCOHPTB,6H  UOJ  ,6H  (Hi  1  / 

WOPOsAKOKO 

IF? inE j.EU.O)  GO  TO  1 
CNHPCF=CNHPOS»SO“TITM/PCNC 

1  CNHP*ChHPCF  «PCNC/  SUR  T  f  1'4 1 
CNHPS=CNHf 
TPf«?S=»TFFHP 

CALL  SEARCH 1-1 . , TrFH? , CKHP ,OHTCHP , ETATHP , 

1TFFXI1 } ,NTf FS.CNX! 1,1 1 ,DH7CXt 1.11 ,E7ATX.C 1*11 .NPTTFF il)»15,13,IGOJ 
IFIIGO.EO.  1 .OK . IGO. Ey . 1 1 .Oft . 1  GO. EQ.21 )  WRlTE{8,I0001TFFHPS,Kn 
IF { 100.60.  2.GR.IG5.  EO.12.OR.1GO.EO.221  WRITE  Co »10001TFrKPS»W>ii 
IFJIGU. eo.10.ORcIGO.EO.il. OR. IGU. 60. 121  WRITE  (6,2000 1  CNHPC-.WLO 
IF(!GO.EO.2OcilR,lGO.EO.21.0«.IG0.E0.221  WRITE (S ,2000  1  CNh?S,*Kl 
lOOO  FORMAT Ci9H0e****7FFHP  OFF  KA?,Fl0.4,2XA6,ilh*«e*»SAS£Ail 
20U0  F0Rr-A7ll?H0**»**  ChHP  OFF  HAP,F10.4,2XA&,11H»*«*«**$*S.*1 

IFtlG0.NE.7I  GO  TO  3 

2  CALL  ERROR 
RETURN 

3  HAPGO«0 

IF  I ASSJ  TFFmPS-TFFmP 1. 1  E.0»O01sT?FHPS  1  GO  TO  <• 

RAft&O*! 

IF { ASSi CNHPS-CMi? I .GT .0.001 sCNHPS 1  HAPG0=3 
GO  TO  5 

<•  ;FtA6SfC«HPS-CNHPl.GT.0.001»C»iHP31  1*APSG=2 

5  IF (i*A?GU.&T.G)  CALL  HAPBACf 1 ,HAPGU,Tr FKPS ,TFFHP,CNMPS,CNHPsPCNC, 

174 .MODE ,OG*AP .NUSHAP 1 
IFfKTtNAP.GT.Ol  RETURN 
7 f  NCAL *WGAOSORT  {  TA 1  / 1 14. f>V6*PA  ! 

BTUEX7a0.7O6705«HPEXT 

OhTCC* ( HTUFXT+wAC® ( K3-H21 1 1/ ( WG4*T4) 

Ifn0ES.ty.01  CO  TO  6 
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TFHPCF=TFhPDS/TFHCAL 

f)HHPCF=OHTCC/OHTCHP 

ETHPCF=ETHPOS/£TATHP 

WRITE (6  r 102  JCNHPCF  ,  TFHPCF ,  ETHPCF ,DHHPCF 
102  FORHAT {20H0H»P«  TURBINE  DESIGN, 5X7HCNHPCF=,Elb. 8, BH  TFHPCF=,E15.B, 
lbH  ETHPCF=,E15.8,8H  DHHPCF= rEl 3.8 ) 

6  TFHCAL=TFHPCF«'TFHCAL 
DHTCHP=DHHPCF*DHTCHP 
ETATHPi-E  IHPCF^ETATHP 
DHTC=DHTCC*T4 

ERR ( 1 J  =  { TFHCAL-TFFHP ) /TFHCAL 
ERR ( 2 ) = {DHTCC-DHTCHP I/OHTCC 

call  THTURSI DHTC, ETATHP, FAR4, K4»S4 , P4, T5  »H5 ,S5 ,P5 ) 

IFUJLHP.LE.O.)  GO  TO  ? 
r *R5=WF8/ (WA3+8LHP) 

WG5=KG4+BLHP 

H3=C8LHP«H3+WG4*H5)/WGf> 

CALL  THERMO {P5th5,T5, S3, XX2,1, EARS, 1) 

GO  TO  8 

7  FAR5=FAK4 
HG5=WG4 

8  CALL  COLPTB 
RETURN 

END 


29 
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SIBFTt  COLPTB  OECK,N94/2,XR7 
SUBROUTINE  COLPTB 


COMMON 

iMORD 

2ISASHX 

BITRYS 

COMMON 

1PCNFGU 

2ZFDS 

3ZCOS 

4T40S 

5TFHPOS 

fcTFLFDS 

7T2ADS 

8T7DS 

9A55 

APS55 

COMMON 


/  ALL/ 

JOES  t JOES  ,KO£S  ,POOE  ,INIT  ,1DUNP  ,IANTP 
IOBURN.IAFTBNjIOCD  flhCU  , IDSKOC»IHSHUC,KOZfLT 
LOOPER  «JKAP  ,HU.MNAP,MAPEDG»T0LALLiERRC6J 
/DESIGN/ 

PCNCGU,T4GU  ,DUMD1  ,OUHD2  ,OElFG  tDELFH  ,DELSFC 
PCNFOS»PRFOS  ,ETAFDS,*AFDS  ,PRFCF  ,ETAFCF ,WAFCF 
PCNCOS.PRCUS  tETACOS.WACDS  ,PRCCF  ,ETACCF  .WACCF 
HFBOS  ,DTCODS, ETABDS, WA3CDS .OPCODS ,DTCOCF .ETASCF 
CKHPOS*ETHPi;S*TFHPCP,CHHPCFtETHPCF*DHHPCFIT20S 
CNLPDS»ETLPDS»TFLPCF,CNIPCF»STLPCF,DHLPCF  ,T210S 
HFCOS  ,DTDU0S»eTADDS»MA230S,0PDU0S  »0700£F  tETAOCF 
WFAOS  tOTAFOSt ETAAOSt >GdCOStOPAFOS  tOTAFCF  tETAACF 
A25  »A6  »A7  ,A8  ,A9  ,A28  ,A29 

AKS5  »CVDNUZ,CVMNOZ« A8SAV  ,A9SAV  ,A23SAV»4295AV 
/  FKttNT/ 


1T1 

,P1 

»H1 

,Si 

,T2 

,P2 

,H2 

,52 

2T21 

,P21 

.>11 

,521 

,T3 

,P3 

,H3 

,S3 

3U 

«H4 

?S^ 

,T5 

,?5 

,HS 

,S5 

AT  55 

»P55 

»H55 

,S55 

,8LP 

,  8LC 

,BLDU 

,BL05 

5CNF 

fPRF 

,ETAF 

,MAFC 

,  RAF 

,MA3 

,PG4 

,FAR4 

6CNC 

»PRC 

»ETAC 

,MACC 

,WAC 

,  ETAS 

tOPCOM 

,OUMF 

7CNHP 

,£TATHP 

, DHTCMP 

i  DHT  c 

,8LKP 

,WS5 

,FAR5 

,cs 

6CNLP 

,ETATLP,OHTCLP 

,OHTF 

,8LLP 

,  KG55 

,FARS5 

,HPEXT 

9AM 

,ALTP 

,ETAR 

,ZF 

,PCNF 

*2C 

,PCHC 

»*FB 

ATFFHP 

tTFFLP 

,PC8LF 

,PC3LC 

,PCBLOu 

,PCBLG8 

,  PC?.  LHP 

,PCBLLP 

1000 

2000 


C OHM ON  /LTURB/TFFX C IS ) ,CNX I 1 5 , 1 5 1 .OHTCX I 1 5 , 1 5 ) ,ETATX (IS , IS ) , 
1NTFFS.HPTTFFC151 

OAT*  AttORD,WLO,MHI/6HCOLPT3,6H  ILOI  ,6H  (HI)  / 

MOP.O-AWORO 

IF { IDES.EO.O)  GO  TO  i 
CN*.PCFxCNLPOS*SORT  < T5> /PCHF 
CHLP*CNLPCF*PCNF/SURT  C  75  ) 

CMLPS*CNLP 

TFFLPS*TFFLP 

CAUL  SEARCH l— 1,  »TFFLP,CNLP»DHTCLP*ETATLPs 

ITrFXIi )*>>TFPStCHXIl«l } »DHTCX( 1  flit  ETATXJ 1 >1 ) ,HPTTrf ill* 15,15,1601 
1F(IG0.E0.  1 .OR.IGO.EO.il .OR. 1 GO. EO. 21 1  WR]TE(B,1000)TFFLPS,WL9 
IFHGO.EO.  2.0R.IGO.E0.12.0R.IGC.E0.22)  KRITE(8,10Q0)TFFLPS»WHI 
ZFUGO.E&.IO.OR.IGO. EO.il. OR. IGO.EO.12I  KRITEI8, 20001  Cffl-PS.MLO 
IF  1 1GO.6Q.20.OR.1GO.EO.21.CR .ISC. E0.22 5  KKITE C3»2000)  CHIPS, XHI 
FORMAT ( 19H0****«TFFLP  OFF  MAP,F10v<>,2X£6,llK**«**SSSS»M 
FORMAT US«0*****  CHLP  OFF  MAP,F10*4,2XA6,11h»»»««S*$*5S; 
IFUG0.NE.71  GO  TO  3 
CALL  ERROR 
RETURN 
MAPGOsO 

IF(A3SCTFFLPS-TFFL?l.LE.0.C01*TrFLPSl  GO  TO  A 


MAFG0»1 

IFIABSI«LPS-£MLP1.GT.0,001*CNLPS)  MAPG0=3 
60  TO  5 

A  Ir (ABS1CMLPS— CHLPJ .GT«C»001*CMLPS)  HAPG0*2 

5  1FIHAPG0.ST.01  CALL  HAPBACl2,KAPG0,TFFLPS,TFFLP,CNLPS,Ctt.P,PCHF, 

IT*, MODE , NOMAP ,NUMMAP I 
IF INOMAP .ST.G1  RETURN 
TFLCAL*«G5*S0RTIT51/HA.696»F5) 

OHTCF-M AF*  C  H21-H2 » / { PGS*7  5 1 
IFCIPES.EO.O)  O  TO  6 
7FLPCF=TFLP0S/TFLCAL 
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DHL  PCF  =DHTCF/l)HTCL  P 

etlpcf=etlpds/etatlp 

HRITH(6,102iCNLPCF,7FLPCF,fcTLPCF,DHLPCF 
102  FURMAT (20H0L.P.  TURBINE  OESI GN, bX7HCNLPCF=,E15.8 »ftH  TFLPCF=, Elb.8, 
1HH  ETLPCF=,Ex5.8,8H  DHLPCF= , El  5.8 ) 

6  TFLCAL=TFLPCF*TFLCAL 
DHTCLP=OHLPCF*UHTCLP 
ETATLP=ETLPCF*ETATLP 
DHTF=DHTCF*T5 

ERR (3 )= { TFLCAL-TFFLP ) /TFLCAL 
ERR (4 )  =  CDHTCF-DHTCLP  J/OWTCF 

CALL  THTURBI DHTF  « ETATLP»  FAK5  »H3 » S5*  P5tT55  »H55  »S55»P55 ) 
IFIBLLP.LE.O.)  60  TO  7 
FAR55=WFB/ ( WA3+BLHP+BLLP > 

WG55=WG5+BLLP 

H55= ( BLLP*H3+HG5*K55 ) /MG55 

CALL  THERMO IP55,H55,T55*S55»XX2,1,FAR55,1) 

GO  TO  3 

7  FAR55=FAR5 
WG55=WG5 

8  CALL  FRTOSD 
RETURN 

END 
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*I8PTC  FRTOSD  DECK* M94/2*XR7 
SUBROUTINE  FRTOSD 
COMMON  /  FRONT/ 


1T1 

tPI 

»H1 

»S1 

2T21 

,P21 

»H2I 

*  521 

3T4 

*P4 

»H4 

tS4 

4T55 

»P55 

,Hd5 

,S55 

5CNF 

,PRF 

,ETAF 

»HAFC 

6CNC 

,?RC 

,ETAC 

,  sACC 

7CNHP 

.ETATHP 

tBHTCHP 

»OHTC 

8CNLP 

tETATLP 

,Ob'TCLP 

tDHTF 

9AM 

*ALTP 

»etar 

*ZF 

ATFFHP 

♦TFFLP 

»PC3LF 

*  PCBLC 

COMMON 

/  SIDE/ 

XXP1 

,XWAf 

*XWAC 

,xeiF 

XXT21 

t  XP2I 

»X>:21 

»XS2I 

3T24 

*P24 

*H2  4 

,S24 

4T28 

,P2B 

,H28 

*S28 

5MAD 

,Wr0 

*WG24 

*  FAR24 

6TS28 

»PS28 

*¥28 

»AHZ8 

X?1*P1 

XWAF=KAF 

XWAC=WAC 

X3LF=BLF 

XBLOU=BLDtl 

XH3«H3 

XT2I=T21 

XP21-P21 

XH21=H21 

XS2I=S21 

CALL  CODUCT 

RETURN 


*T2 

»P2 

*H2 

»S2 

»T3 

*P3 

,H3 

?S3 

»T5 

*P5 

»H3 

*Sb 

.BLF 

,blc 

,3  LOU 

iBLUb 

»  waf 

»WA3 

»*G4 

»FAR*» 

»  MAC 

♦  ETAB 

»OPCOM 

,DUMF 

♦  BLHP 

,WGS 

*FAR5 

*CS 

»BLLP 

*KGS>5 

*FAR55 

*HPEXT 

*PCNF 

*2C 

,PCNC 

»  WFB 

*  PC8LDU 

iPCBLOB 

»PC6«.HP 

♦PCBLL^ 

«  X3LDU 

,XH3 

*  DUMS 1 

»UUHS2 

*  T23 

»P23 

»H23 

»S23 

» T25 

,P25 

*H23 

*523 

,T29 

♦  ?29 

*H29 

»S29 

*  ETAD 

,OPOUC 

» BYPASS 

, DUMS 3 

»TS29 

,PS29 

*¥29 

*  AM29 
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SI8FTC  COOUCT  0ECK,H94/2*XR7 
SUBROUTINE  CODUCT 


COMMON 

IWORO 

2lG»:*L* 

3ITRYS 

COMMON 

1PCNFGU 

2ZFDS 

3ZCOS 

4T4DS 

5TFHPDS 

6TFLPDS 

7T24DS 

8T7DS 

9A55 

APS55 

COMMON 


/  ALL/ 

IOES  *  JDES  fKOES  ,HOOE  ,INIT  * l DUMP  *IAHTP  * 
103URN  *  I AFTBN  * I DC D  ,INCD  , IDSHOC* IMSHOC.NOZFLT, 
LOOPER.NOHAP  ,NUKHAP,MAPED6*T0LALL,ERR(6) 

/DESIGN/ 

PCNCGU,T4CU  .DUMD1  ,DUKD2  ,DELFG  .DELFN  tDELSFC, 
PCNFDSfPRFOS  *cTAFDS,WAFDS  ,PRFCF  *ETAFCF*KAFCF  * 
PCNCOSyPRCDS  , ETACDS* MACDS  ,PRCCF  »ETACCF*WACCF  » 
WF8DS  i DTCOOS , ETABDS * KA3CDS » DPCOOS * DTCOCF , ETABCF* 
CNHPDS»ETHPDS»TFHPCF,CNHPCF,ETK?CF,0HHPCF,T20S  » 
CNLPDS  »ETLPDS  »TFLPCF iCNLPCF*  ETLPCF  »DHLPCF  »T21DS  * 
MFOOS  , OTOUOS * ETADD5, WA23DS ,OPOUDS tOTDUCF  * ETADCF , 
WFACS  « OTAFOS » ETAAOS , WG6CDS » OPAFDS , DTAFCF , ETAACF , 


*A6 


*A7 


.A8 


*A9 


A25 

AM55  tCVONOZ  tCVMNQZ  * ASSAY  rA9SAV 
/  SIDE/ 


*A28  *A29  * 

*A26SAV»A295AV 


1P1 

*MAF 

*KAC 

*BLF 

*BLDU 

»H3 

jDUHSl 

•DUMS2 

2T21 

*P21 

»H21 

*S21 

»T23 

,?23 

*H23 

»S23 

3T24 

,P24 

,H24 

*S24 

»T25 

»P25 

.K25 

*S25 

4T2S 

,P28 

*H28 

*523 

»T29 

*P29 

»H29 

»S2S 

5HA0 

,KFD 

,NG24 

*FAR24 

*ETAD 

,dpduc 

*BYPASS 

»  DU MS  3 

6TS28 

*PS26 

,V28 

•  AH28 

»TS29 

*PS29 

*Y29 

*AM29 

DIMENSION  0(91 

OATA  AW0RD1 >  AW0R02/ 6KC00UCT *  6HDNQZZL/ 


WORD  -A WORD 1 
Q( 2)«Q* 

0(31-5. 

WAX-WAF-MAC-BLF 

WAD-WAX+5LDU 

P23-P21 

H23-(BL0U*H3+MAX*H21 J/N^D 

CALL  THERM0(P23-H23,T23»S23,XX2,1*0.0»1) 

BYPASS- (WAF-WAC) /MAC 

WA23C-HAD*SQRT ( T23 1 /P23 

IFUDES.EQ.i)  WA23DS-WA23C 

OPDUC '»DP0UDS*(MA25C/WA23DS  J 

IF (DPOUC.GT. I. )  OPDUC-i. 

P24«P23*(1.-0PDUC) 

IFIIGASMX.GT.OJ  I OS URN-0 
IF ( IOBURN.NE.O)  GO  TO  2 
T24-T23 
HFD-O. 

FAR24-0. 


GO  TO  7 

2  IF(IDBURN»EQ.2)  T24-T23+20G0. 

3  IF (T24.GT.4QC0. )  T24-4000. 

IF(T24.LT.T23)  T24-T23 

C  ***  IF  DESIRED*  ENTER  CALCULATIONS  FOR  ETAD  HERE 

HV-( ( ( ( { {-.45943l7E-l9*T24)-.20341l6E-15)*T24+.27fi3&43E-ll )*T24* 

1. 2051501E-07 ) *T24— .2453116E— 03 > *T24-.9433296£-01 ) +T24+. 1S45537E+05 
CALL  THERMO (P24*HA*T24* XXI *XX2*0*0.0*0) 

FAR24-IHA— H23)/(HY*ETAO) 

IFCFAR24.LT.0.)  FAR24-0. 

HFDX-FAR24+VA0 
IF( ICBULH .NE.2)  GO  TO  6 
ERRW- (HFO-HFOX ) /KFO 
0 IR-SORT (WFD/NFDX ) 

CALL  AFQUIR(Qll) »T24»ERRM»0.»?0* *0.0001 »DIR»T24T *IGO) 

GO  TO  (4*7*3)* IGO 


S3 


-I  O'  m 
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4  T24=T24T 

GO  TO  3 
CALL  ERROR 
MFO«HFDX 

CALL  THERMO IP24,H24»T24,S24,XX2,1»FAR24»0) 

HG24*HFD+HAD 

IFCIOES.EG.l)  HRITE(6»101)  KA23DS 
101  FORMAT  U2H0DUCT  DESIGN, 12X8H  WA23DS=,E15.8 ) 

C  ***  IF  DESIRED,  ENTER  OTHER  LOSSES  HERE 
T25*T24 
P25*P24 
H25=H24 
S25=S24 

1F(IGASMX*GT»0)  GO  TG  11 

N0RD=AW0RD2 

A28SAV=A2S 

A29SAV*A29 

NG2D*0 

I0N0Z=0 

IFJN0ZFLT.EQ.2.0R.N0ZFLT.6Q.3)  N0ZD=1 

IF  I IDES.EQ.1.0R»ID8URN,GT.0.0R.N0ZDoEQ. 1)  I0N02=1 

IFCIOCD.EQ.l)  SC  TO  8 

CALL  C0NVRG{~25,H25,P25,S25tFAR24,KG24,Pl,I0N0Z,A28,  P25R, 

1T28,H26,P28,S28,TS28 ,PS28 ,V28?AH28, ICON) 

GO  TO  19,9, 9,5), ICON 

8  CALL  C0NDIViT25,H25,P25,S25sFAR24,KG24,Pl,IDN0Z,A28tA29,P25R, 
lT28,H28,P28,S28,T29,H29,P29,S29,TS2e,TS29,PS28,PS29,V28,V29,AH28, 
2AH29,IC0N) 

IOSHOC=ICON 

GO  TO  (10, 10, 10, 5), ICON 

9  T29*T28 
H29»H28 
P29*P28 
S29=S28 
TS29*TS28 
PS29=?S28 
V29®V28 
AH29=AH28 
A29=A28 
I0SH0C=IC0N+3 

10  ERR  *5 )  =  f  P25R-P25) /P25R 

IF flDNOZ.EC.l )  WRITE!6,100)  A28,AK28,A29,AK29 
ICG  FORMAT {19H0DUC7  NOZZLE  DESIGN, 5X8H  A28=,E15.8,8H  AH23=,E15.8, 

18H  A29=,E15»8 ,8H  AH29=,E15«8> 

11  CALL  FAST8K 
RETORN 

END 
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SIBFTC  FASTBK  DECK,M94/2yXR7 
SUBROUTINE  FASTBK 
COMMON  /  FRONT/ 


1T1 

»PI 

»H1 

r  SI 

»T2 

»P  2 

tH2 

,S2 

2T21 

»P21 

»H21 

yS21 

*T3 

»P3 

tK3 

y  S3 

3T4 

»P4 

»H4 

yS4 

»T5 

»P5 

»H5 

yS5 

4T55 

,P55 

»H55 

»S55 

,BLF 

*  BLC 

tBLDU 

yBLOB 

5CNF 

»PRF 

»E7AF 

»HAFC 

,KAF 

»WA3 

»HG4 

yFAR4 

6CNC 

»PRC 

tETAC 

jHACC 

yWAC 

»ETAB 

,DPCOH 

yOUMF 

7CNHP 

, cTATHP , DHTCHP, DHTC 

yBLHP 

f  WG5 

tFAR5 

>CS 

8CNLP 

tETATLP  tDHTCLP iDHTF 

yBLLP 

»WG55 

»FAR55 

yHPEXT 

9AM 

,ALTP 

»ETAR 

yZF 

tPCNF 

yZC 

yPCNC 

»HFB 

ATFFHP  ,TFFLP  ,PCBLF 
COMMON  /  SIDE/ 

yPCBLC 

♦ PCBLDU , PCBLOB ,PCBLHP ,PCBLLP 

XXP1 

yXHAF 

t  XWAC 

tXBLF 

tXBLOU 

»XH3 

yDUHSl 

yOUKS2 

XXT2I 

yXP2i 

»XH22 

,XS2J 

»T23 

,P23 

t«23 

yS23 

3T24 

t?24 

,K24 

»S24 

»T25 

»P25 

»H25 

yS25 

4T28 

,P28 

t  H28 

*  S28 

,T29 

»P29 

»H29 

tS29 

5MAD 

yWFD 

y*G24 

» FAR24 

»ETAD 

yDPDUC 

»BYPASS 

T DU MS 3 

6T528  »PS26  ,V28 

COMMON  /  BACK/ 

*  AM28 

»TS29 

»  PS29 

yV29 

yAN29 

XXT55 

»XP55 

,XH55 

t  XS55 

,XT25 

,XP25 

»XH25 

yXS25 

XXWF8 

,XHG55 

yXFAR5! 

5,XWFD 

»XKG24 

,XFAR24,XXP1 

yfXJHB 

3T6 

j?6 

«M6 

*  S6 

»T7 

9*7 

»H7 

»S7 

4T8 

»P8 

»H& 

»S8 

•  T9 

»P9 

»H9 

yS9 

SHG6 

t  WFA 

»KG7 

»FAR7 

t  ETA  ' 

yDPAFT 

»V55 

*V25 

6PS6 

tV6 

»  AM6 

»TS7 

,PS7 

»v? 

»AM7 

,AK25 

7TS8 

»PS8 

rV8 

»  AK8 

tTS9 

»PS9 

»V9 

,AM9 

8VA 

»FRD 

»VJD 

yFGHD 

yVJM 

yFGHH 

tFGPD 

tFGPK 

9FGM 

XT  55  = 
XP55= 

»FGP 

=T55 

P55 

?HFT 

jKGT 

yFART 

,fg 

»FN 

ySFC 

XH55=H55 

XS55=S55 

XT25=T25 

XP25=P25 

XH25=H25 

XS25=S25 

XWFB=KFB 

XWG55=WG55 

XFAR55=FAR55 

XHFD=HF0 

XWG24=KG24 

XFAR24=FAR24 

XXPI=P1 

CALL  COHIX 

RETURN 
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1I8F7C  COM I 5  0=CX#K94/2»XR7 

SUBROUTINE  COKIX 
COMMON  /  ALL/ 

»IDES  * JOES  ,KOES  ,HDO£  ,INIT  .IOUMP  ,IANTF  » 
»IOBURN, IAFTBN, IDCD  ,1KCD  t  IOSHOC » 3KSW3C  »N'OZFLT  $ 
#LOOPfcR,NOHAP  ,NUB*AP,KAPEDG,TOLALL,ERRl&> 


1WOR0 

2IGASHX 

3ITRYS 

COMMON 

1FCNFGU 

2ZF0S 

32CDS 

4TADS 

STFHPPS 

6TFLP0S 

7T2+DS 

8T7DS 


/DESIGN/ 
PCNC60»T4GU 
PCKFDS,?RFOS 
PCNCDS.PRCOS 


jOUHOl  ,DU*02  tDELFG  ,OELFN  ,DELSFC 
,ETAFDS,HAFDS  ,PRFCF  ,ETAFCF»*AFCF 
,ETACOS,WACOS  *PRCCF  ,ETACCF,WACCF 
NF80S  .  OTCODS  ,  ETA3DS , HA3C0S , OPCOOS , DTCOCF ,  ETA5C  F 
CNHPDS  t  ETHPDS , TFHPC  F, CNHPC  F , ETKPCF , OHHPCF  « T20S 
CNLPDSt£TLPDSiTFLPCF,Cm.PCF,ETLPCF,DHLPCFfT21DS 
WFDOS  t  DTCUDS  « ETA  DOS , WA23DS  ,  DPOUDS  » DTOUCF , ETADCF 
WFA05  »DTAFDS»ETASDS»W&6CDS»DPAF05 »OTAFCF,ETAACF 


9A55 

?A25 

tA6 

»A7 

»A8 

»A9 

,A2  8 

»A29  , 

AP555  ,AM55  tCVONOZ 

COMMON  /  SACK/ 

»CVMN02 ,  ASSAY 

,A9SAV 

,A28SAV,A29SAV 

IT55 

,P55 

,K55 

,555 

•  T25 

»P23 

jH25 

*S25 

2WFB 

,*G55 

, FAR55 

,  HR) 

,WG24 

, FAR24 

,PI 

jDUHB 

3T6 

,P& 

,K6 

»S6 

tT7 

,P7 

,H7 

,ST 

4T8 

*P8 

,H8 

»S8 

tT9 

,P9 

,K9 

,S9 

5*56 

,KFA 

,*&7 

,FAR7 

jETAA 

,OPAFT 

»V55 

,V25 

6PS6 

»V6 

»  AKA 

,TS7 

,PS7 

,V7 

,AN7 

,AH25 

7T5B 

,PSS 

,  VS 

,AK8 

»TS9 

,PS9 

»V9 

,fc*9 

8VA 

»FRD 

*¥JD 

,FGND 

»VJH 

»FG*H 

,?&PD 

,FGPM 

9FGN 

DATA 

,FGP  ,WFT 

AW0RD/6H  COM IX/ 

,HGT 

,PART 

,FG 

,FN 

»SFC 

DIMENSION  0(1  £91 

WOXD*A*OAD 

AJ*773.26 

CAPSF»2II6,2170 

G*32.1740A9 

CALL  PR3COMIFAR5SIT55,BH,XX2,XX2,XXA,PH!5S,XX5> 
CALL  PR0C0HIFAR24»T25»XXI »XX2,XX3,XX+»PHI25tXX5 1 
IF f IDE3.EO.O)  GO  TO  6 

C  ***  CALCULATE  A55  AND  A?5  WITH  PS25«PS55 
IF(PS55.60.0.)  GO  TO  50 
T$55«T55*(PS55/P551**0.286 


DO  1  I»I»15 

CALL  PAOCOM{FAR55,TS55*CS55r AK55, CP55»ReX55,PHIS55,H555) 

?HIS*FHI55~REX55*AL0G(P55/PS55J 

DELPHI =PHIS~PHIS55 

IF (A5S (DELPHI l,LE«G«OOOi*PHIS 1  GO  TO  3 
TS55»TS55*EXPt4.0*DELPHi I 


2  CALL  ERROR 

RETURN 

50  TS55aO»875»T55 
00  SI  1*1, !5 

CALL  PROC0HIFAR55,TS55,CS55,AK5S,CP53,R£X55,PHIS55,HS55J 

VS5*AH55*CS55 

HSCAL*H55— V55**2/ I 2« *S»AJ } 

0ELHS*HSCAL-«SS5 

IF l A3S (CcLHS 5 .L  E. 0.0005*H3CAL )  GO  TO  52 

51  TS55*TC55+OELHS/CP55 
GO  TO  2 


52  P555*P55/EXP((PHIS5-?HISS5)/REX551 

3  IF(HS5.6T.HS55»  SO  TC  53 

HRlTE!S,10nP55tPS55,TS5,TS55,H55,H355 
iOl  FORMAT i 22H0S0RT  OF  H55-H5S5  NEG  ,6£15«8,6KiS*tit ) 

CALL  ERROR 


38 
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53  V5S=SCRT{2.»G<'A^{H55-KS55)J 
RK0=CAPSF*PS55/(A,i*F.EX55*TS55} 

A55=WS55/(RH0*V5S) 

AK55=V5S/CS55 
I*(IGASNX.GT.C'  SO  TO  54 
WRITE{6»iOMA55>AN55 

104  FORMAT <20H0TUR8IN’£  AREA  0ESISW,6X6H  A55*»E15.8»8H  AK55=,Ei5.8) 

50  TO  23 

54  PS2S=PS55 
7S25=Y2S*«PS25/P25>»*0.286 
DO  4  1=1, *.5 

CALL  PRO£03{FAR24,TS25,CS25,AtC25,CP2S,REX25,PKlS25»HS25) 

PHI 3*PHI 2S-R£X25*AL0G  (P25/PS25) 

DELPHI=PHI3“?HIS25 

IF{A8S(OEL?HI).LE.0.0001*PHIS)  60  TO  5 

4  TS25=T525*EXP(4.0*DELPNI) 

GO  TO  2 

5  IF(H25.GT.HS2S)  GO  TO  S5 
WRITE (8,102 »P25»PS25tT25,7S25,H25»HS25 

102  FORMAT ( 22H03CRT  OF  H25-HS25  NEG  ,6E15.6,6HSSS*SS ) 

CALL  ERROR 

55  V25=S0RT(2.*G*AJ»CH25-fii?55 ) 

RHP=CAPSF*PS25/(  AJ»R.£3C25*TS?5) 

A25=KG24/(RH0*V25) 

A.M25=y25/CS25 

WRITE(6,100IA55,AM55,A25,AH2S 
100  ?0RHA7 ( 25KQTURB I JiE/DUL'T  ARE4  DESI6N,7K  A55=,£15.3, 

18H  AH55=»E15~8,8K  A25=?E15.8,8H  AA25*,E15.8) 

GO  TO  20 

C  ***  CALCULATE  PS55  AND  PS25 

6  W3A=VGS5/A55 
C1=P55<50RT IS/ ( T55*A J 1 ) *CAPSF 
KC  05**0 
Q0(21*=0. 

00(33=0. 

AK55=0.S0 
TS55=0.375*TS5 

7  00  8  1=1,15 

CALL  PR0C0JiIFAR55,TS55,CS5S,AK55,CP5SrREX55,PHI-S55.HS55) 
V55=aS55*CS55 

HSCAL=H55-V55**2/  ( 2,  4-G*AJ } 

OELHS=HSCAL~t*5S5 

IF ( A3S(OELHS*.LF.0.0005*H5CAL)  GO  TO  9 
3  TS55=TS55+OELKS/CP55 

GO  TO  2 

9  WOAT=C 1*SQRT(AK55/RSX55 )*AH55/( I.+C AXS5-1. )*AM55**2/2. )*» 

1(( AK55+1. )/(2.*( AK55-1. 13) 

AKX-AH55 
I  GO  GO  =0 

10  OIP  =WOA/WOAT 
E*  =(WCA— WQATl/WQA 

CALL  AFOUIR (QQ(1 ) ,AHX,EW, C«, 30., 0.0005, DIR ,AHXT ,ICON ) 

GO  TO  {11,15,2),IC0.N 

11  IF ( AM XT. LE. 1.0)  GO  TO  13 

AH XT -0.7 
MC0K=RC0M+1 

IFCMCOH.LE.l)  GC  TO  13 
PCNF=0UM01 

WRITS (3,103)PCNF»AKX»P55»PS55,P25tPS2S 

103  FORMAT ( 22HOCQMIX  PO*F=,F7.4,4H  AH=,F8.6,5rt  P35=,r9.5, 

16H  PSS5=,F9.5,5H  ?25=,F9.5,8H  PS23=,F9.5,&HSS«SJI 
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AFAPL-TR-87-125 
Pari;  n 


PCNF-1.01*PCNF 

0UM01-PCNF 

12  130MAP-7 
RETURN 

13  ?F' IGQGO.CO.il  GO  TO  14 

AM55-AMXT 

GO  TO  7 

14  AM25-AMXT 
GO  TO  16 

15  IF ( IGOGQ.EQ. 1 )  GO  TO  19 

PS55-P53/EXP( IPHI55-PH1S55) /REX55) 

IF(IGASMX.LE.O)  GO  TO  28 
WOA»WG2«/A25 

Cl-R25*SeRTIG/<T25*AJI 1*CAPSF 
MCON-O 
00(21-0. 

00(31-0. 

AM25-0.25 

TS25-0.875*T25 

16  00  17  1-1,15 

CALL  PROCOM(FAR24,TS25,CS25,AK25,CP25,RlX25,PHIS25,HS251 

V25-AK25*CS25 

HSCAL-H25-V25**2/ ( 2 . *G*AJ 1 

DELHS-HSCAL-HS25 

IF ( A6S (OELHS 1 »LE .0 .0005*HSCAL1  GO  TO  IS 

17  TS25-TS25+0ELHS/CP25 
GO  TO  2 

18  WQAT-C1*SGRT( AK25/REX25)*AM25/ ( l.+( AK25-1. 1 *AH25**2/2 •  >** 
1( 'AX25+1. 1/(2. 91AK2 5-1. Ill 

AMX-AH25 
IGOGO-1 
GC  TO  10 

19  PS25-P25/EXP ( (PHI25-PHI525) /REX251 

20  WG6-WG24+HG55 

ERR (5 1" (PS25-PS551 /PS25 

WF6-MF0+MFB 

PAR6-WF6/ (WG6— WF61 

H6- ( WG24*M25+WG55*H55 1/MG6 

CALL  THERMO I 1. ,H6,T6 ,PHI6,AMX, 1 , FAR&,  11 

C1-PS55*A55*( l.+AK55*AM55**2 )+P525*A25*(  1  ••(■AK25-AM25  ->*2 1 

TS6-0.833*T6 

('0  25  1-1,15 

CALL  PROCOM (FAR6 »TS6*CS6, AK6fXP6 ,REXft»PHIS6*HS6 1 

C2-WG6*SQRT ( A J*REX6*T6/ ( AX6*G 1 1 

C3-C2/ (CAPSF*C1 1 

C4-CAK6-1. 1/2.— (C3*AK6 1**2 

C5U.-2.*AK6*C3**2 

C6-C5**2+4.*C4*C3**2 

IF (C6I21,22,23 

21  CALL  ERROR 
RETURN 

22  AM62—C5/  ( 2.*C<»1 
GO  TU  24 

23  AM62- ( SORT ( C6 1— C5 1 / ( 2 . *C4 ) 

24  1FCAM62.LE.0.1  GO  TO  21 

AM6-5QRT ( AM62 1 
Vfc-AM6*CS6 

HSCAL-H6— V6**2/ <2.*G< AJ) 

DELHS-HSCAL-HS6 

IF ( ABS (OELHS 1 »LE« 0.0005*H5CAL 1  GO  TO  26 

25  TS6-TS6+0ELHS/C.P6 
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AFAPL-TR-67-125 
Part  II 


! 


GO  TO  21 

26  IFUGASMX.EQ.2)  GO  TO  27 
A6=A25+A55 

27  C7=SQRT ( 1.+IAK6-1. }*AM62/2*  ) 
PS6=C2/(CAPSF*A6$AH6*C7J 
P6=PS6*EXP{ (PH16-PHIS6J/REX6) 

CALL  THERMO! P6,H6,T6tS6, XXI,  ly.6AR6»0) 
S6AV£= (WG24*S25+WG55*$55 ) /WG6 
I F { S6 « GE « S6A VE )  GO  TO  29 
S6=S6AVE 

P6=EXP ( AMX*(?HI6-S6) /I .986375  3 
GO  TO  29 

28  T6=T55 
P6=P55 
H6-H55 


So-555 

HG6=WGS5 

P56=PS55 

V6=V55 

AM6=AH55 

IF { IGaSMX.EQ.01  A6=A55 
29  CALL  C0AF8N 
RETURN 
END 


l 

5 

I 

! 

t 


i 

i 
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AFAPL-TR-67-125 

Partn 


SiSFTC  C0AP5K  DECK, *34/2, XR7 
SUSRCyTIHg  CDAF&N 
COMB*  /  iU/ 

1V0SD  »I0£$  (JOES  iKDcS  ,800c  .tH.'T  ,  JOUMP  .IAHTP  • 

ft*^*%I^o^oT*«FTf*,ICC0  ,IKC0  » 1 DSHQC  » IP^HOC  tKQZFLT  j 
3ITr»YS  fLOOPER.HCMAP  .mJKMAP.HAPEDGiTQLALL.ERRto J 


20 


22 


24 


COMMON  /DESIGN/ 
1PCNFGU .PCHCGU, T4GU 
2ZF0S 
3ZCQS 


.DUHDl  .OUK02  tOELFG  .D^LFN  .DELSFC 
tFCNFDS »PRFDS  .ETAFDS.WAFDS  .PRFCF  .ETAf  IF .WAFCF 
»fCNCOS,PRCDS  .ETACDS.WACDS  .PRCCF  .ETACCF.KACCF 
4T40S  tkrSGS  jOTCOOSj  STABD§ryA3C0S  tOPCODS  »£rTCOCF , ETABCF 
5TF8F0S  ,£3HJ>05,  ETrtPOS  ,7FHPCF  .Ct&tPCF ,  ETHPCF  .DHriPCF  .TEDS 

l^^^?S*H-H?StETJ‘?DS,Tfi'PCf:tC,iLPCP’£T-PCf:*DKl-pCF»T21DS 
7«240S  tfcrOOS  »DTOODS?ETADDSt KA230S. DPDUDS .DTCUCF .ETADC ~ 
8T7DS  .HFADS  «  DTAFDS,  ETAAQS. WG6C3S  .  DPAFDS  .D7AFCF  *  FT AACF 
»A25  ,A6  ,A7  .AS  |A9  .£28  ,429 

AFS5s  ,AM55  t C VONOZ , C  YWi02 , A SS«  V  .A9SAV  tA26SAV,A2s>SAV 
COHMGN  /  BACK/ 


1T5S 

2WFB 

3T6 

478 

5HG& 

6PS6 

7TS8 

SVA 

9F&M 


»H55 

.FAR5S 

*H6 

?  riS 

,KS7 

,AMS 

*Y8 

sVJD 

,HF7 


.255 

«y?o 

.58 

rSS 


i  ~A:- 


,7S7 

;A«S 

,?GS3 

,«s? 


tT25 

,ftG24 

»T7 

,?9 

,rTAA 

«PS7 

,T£9 

,VJ« 

,FAR7 


»P25 

.FAR24 

»P7 

,P9 

.DPAFT 

»Y7 

,PS9 

.FGMM 

,FG 


»H25 

.PI 

»H7 

»H9 

.V55 

,AM7 

.V9 

.FGPO 

.FH 


.  S25 

.DU  M3 

.57 

.59 

»V25 

,AM25 

.  AM9 

.FGPK 

,SFC 


C  m 


.P55 
,XG55 
»P6 
»?S 
jWFA 
,V& 

.PS8 
,FRD 
.FGP 

DIMENSION  Qi9l 
DATA  AWJRC/6HC0AFBH/ 
tfORD-AWORD 
Qt  21*0. 

0131*0. 

AJ*7?8.26 

CAPSF-2116.2I70 

0*32.174049 

MF6*UFB 

IFUGASMX.G7.01  MFj=WFt+!<?0 

WA6c-^6-MF6 

DRY  LOSS 

MP6C *KG6«SQRT (T61/P6 
lFHDES.EO.il  UGoC0S»KG6C 
0PAFT*0PAF0S*lKG6C/WG6e0S  > 

IF10PAFT.GT.1.1  DPAFT* 1. 

P7*P6«U.— DPAFT  1 
A7»A6 

FAR6*WF6/ScA6 

CALL  PR0C0M{FAR6»  f6»XXl,XX2»XX3«XX4»PHI6»XX6} 
WQA*MG6/A? 

Cl*P7«SQRTlG/( :6*AJ})*CAPSF 

AH7*AN6 

7S7=0.87S*T6 

00  22  1*1,15 

CALL  PROCOH1FAR6.TS7.CS7.AK7.CP7.REX7.PHIS7 .'IS?) 
V7*A«?*CS7 

HS CAL *H6-V7««2/ (2. *G*AJ) 

0ELHS*HSCAL-HS7 

IT (ASS (OELHS) .LE.O«0005*HSCAL)  GO  TO  24 

TS?*TS7+0ELHS/CP7 

GO  TO  8 

HQAT*Cl«SOftT(AK7/REX7}aAM?/Cl.*tAX7-l.)«AM7**2/2.)*» 

1CCAK7+1. J/I2.*(AX7-1. 11) 
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DIR*HQA/K£WT 

6WeCWQa-*QAT//«0A 

CALL  AFCOIR{0:'2},A.''..',£K.O»f30.ft*.CC05r0IRjAK7T,ia0J 
GO  TO  U*6,28,-3)»tSO 
26  AM7=AH7T 

IFCA«7.SE.i.9)  A«7*0.9 
GO  TO  20 

28  PS?»?7/£>:P(fPHI&-?KlS7)/RcX7! 

IFCiAFT8ft.GT.0J  GO  TO  4 
C  ***  HOH-aFTERBUSW  IftG 
3  T7=T6 

WFA=O.G 
FAR7*FAR6 
VG7=KG6 
GO  TO  13 

C  «**  AFTER  BURN  JffG 
K  1F{IAFT8N.£«j2)  77=7^2000. 

IF  (T7.LE.TS  l  SO  7C»  3 
RK065=CAPSF*fiS7/C;W*?.ex7*TS7) 

PS65=PS7 

V65*V7 

0(21=0. 

Q(3)=0. 

5  IF (T7.G7.4Q00. i  77=4300. 

C  ***  IF  DESIRED,  ENTER  £AtCL*LA7!CjfcS  FOR  STAA  HERE 

HV*(f  (CCC-.4594317E-.i«ST7J-.20341l6e~i5}*T7-i  .27S36*35-li>»T7 
l>.205150lE-07J«T7-.^*5TU6=“05f*7?“c9A33296E-e-l  i»T7^.i&«tS537E+G5 
CALL  TKERXO(P7,HAf7?*Xal<XX2:S,O.0»SI 
FAR7= (HA-H6  J / (HV*ETA£ } 

IFCFAR7.GT.0.J  GO  TO  6 

T7=»T6 

GO  TO  5 

6  KFAX=FAR7*HG6 
IFCIAFTBN.EO.I i  GO  70  9 
ERRW*(WFA-HFAX :/Hr A 
DIR=S03T  CHFA/Jt?AA  J 

CALL  AFOUIRCCI I J,T?j£fvRS< f *7.120.70. C001  *0IR*T?7r ISO? 

GO  TO  C7,10,8)r!GG 

7  7?«=77T 
GO  TO  5 

8  CALL  ERROR 

9  WPA=«FAX 

10  FAP.7*=(WF6*KrAJ/VA6 
WG7«=VG6*HfA 

C  «•  HCMENTUK  LOSS 

CALI  PROCOJtIFARTfT?  -XX l,XX2TXX3»RtX7,PHI7 »H7) 
RH07=CA?SF*P7/jAJ*R£X7«T7| 

V7-iiG7/CRK07=A75 
Q  ( 2  )  =0  . 

QJ3i=0« 

?S7»PS65— 0.01 

11  RHGT=WS7/IV7»A7| 

HS7eH7-Y7«*2/ 1 2.  *$=A-J ) 

CALL  rHERN0{1.0,HS7,TS7fPH;57.XX2,I,FAR7,l) 

IFCTS7.SS.301. )  GQ  73  110 

CALL  TH=R«0C1.0tHS7r*0C.,PHIS7fXS2f 1,FAR7,0) 
V7«=SGRTC2.*G*AJ»(H7-HS7J  ) 

GO  TO  11 

110  PS7*RH07*AJ*REX?*7S7/CAPSF 

FS  7A-PSS5+  CRH065*V65»*2-RHC7*V7**2 )  /  C  C-^CAPSF  I 
OJSeSORTC  A3S-:FS?/PS7AJ  ] 


<1 


..  •fit-.' 


A?A?L=TS-5?-IS5 
Part  II 


'  -■:* ' 


V  r 


l 

A 


EP*{PS7-PS7A)/PS7 

CALL  AFQUIR (QC1) ,V7,EP,0., 50., 0.0005, DIR, V7T.IG0) 
Y7*V7T 

IFIV7.LT.100.)  V7=i00. 

60  TO  (ll»12«6),2GO 

12  P7«PS7*£XP(CPHI7-PHIS7>/REX7) 

CALL  PROCOH  C  FAR? , TS7 , CS7 , XX2 , XX3 , XX4 , XX5 , XX6 ) 
AH7*V7/CS7 

13  CALL  THSRM0CP7,H7,T7,S7,XX2.1,FAR?,0J 
IF(IDES«EQ*1 )  KRITE(6»1C0>  K66CDS 

100  FORMAT (19H0AFTERBURNER  DESIGN, 5X8H  KG6CDS=,E15.8> 
CALL  COHNGZ 
RETURN 
END 
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A  S'  A  DT  -TR-S?  ^  12  5 
Part  if 


vs*f“^ 

Li3 


SlBFTC  C0MN9Z  DECK,M94/2,XR7 
SUBROUTINE  C0MN02 
COMMON  /  ALL/ 

1W0RD  ,IDES  » JOES  »KDES  ,MOOE  ,IN1T  »IDUHP  »1AMTP  » 
2IGASHX, 1DBURN, IAFTBN, IOCO  ,IKCO  t IOSHOC  > I HSHOC •NQZFLT  Y 
3ITRYS  |LOOPER».NOMA?  »NUKMAPf HAPEDGiTOLALL »ERR (6 } 

COMMON  /DESIGN/ 

1PCHFGU,PCNCGU,T4GU  »DUH01  »DUKD2  ,0ELFG  ,DELFN  ,DELSFC, 
2ZF0S  ,?CNFOSfPRFOS  ,ETAFDS,WAFDS  » PRFCF  .STAFCF ,WAFC r  . 
32CDS  , PCNCDS.PRCD5  ,£TACDS , MACDS  »PRCCF  jETACCF ,MACCF  , 
4T49S  »MFe9S  »DTC003t ETABDS»MA3CDS tDPCQDS iOTCOCF»ETABCF , 
5TFHPOS , CNHPDS i ETHPOS . TFHPCF , CNHPCF , ETHPCF ,OHKPCF , T2PS  r 
6TFLPDS,CNLPDS,ETLP0StTFLPCF,CKLPCF,ETLPGF,DHLPCF-*210S  , 
7T24DS  rHFCDS  ,DTDUDS,ETA0DS,HA23DS,0PDUDS,DT0L*CFcCTADCF, 
8770S  »HFADS  , OTAFDS  ,  ETAAOS , WG6CDS , OPAFOS  rOTA*CF , ETAACF , 
9A55  ,A25  ,A6  *A7  »A&  ,A9  ,A28  »A29  » 

APS55  ,AH55  ,CVONOZ.CVHNOZ»A8SAV  fA9SAV  ;A23SAV,A29$AV 
COMMON  /  BACK/ 


1T55 

»P55 

,H55 

tS55 

,T25 

,P25 

♦H25 

»S25 

9 

2HFB 

•VG55 

,FAR5S 

,HS2A 

,  FAR24 

,DUM3 

f 

3T6 

,P6 

,H6 

,$6 

«T7 

,P7 

»N7 

»S7 

» 

ATS 

,P8 

,H8 

•  SB 

,T9 

»P9 

,h9 

9 

5KG6 

,kf  A 

»WG7 

,FAR7 

•  ETAA 

-DPAFT 

»V55 

vV25 

1 

fcPSfc 

,V6 

,AM6 

,TS7 

•  PS7 

*Y7 

*AH7 

»AK25 

9 

rrss 

,PS8 

,V3 

,ANB 

,TS9 

»PS9 

?V9 

»AM9 

* 

6VA 

,  FRO 

,VJ0 

,FC*«n 

,VJK 

,FGMH 

,FGFO 

,F3PH 

9 

9FSH 

,FGP 

»VFT 

»  MGT 

,FART 

»F6 

*fm 

.s?c 

OATA  AH0RD/6HHN0ZZL/ 

WORD=AWGRO 
A8SAY=A8 
* 9SAV=A9 
>2>ZM=0 
'KN0Z=0 

IF (NOZFLT .ECel.0R,S02FLT»EQ»3)  KOZM=I 
IFIIOFS«EO»1.0R.IAFT8N.$T.O,OR.N02H.EO.U  IMNOZ*! 

IFdMCO.iO.U  GO  TO  1 

CALL  C9NVRG(T7,rt7*P?.S?,FAR7,KG7,Pl,l«NOZtA5*  P?R, 
lTe,H8,P8,S3.TS5,PSB,V8-#£H3f  ICON) 

GO  70  (3, 3,3, 2), ICON 

CALL  C0N0IY,T7,H7,P?,S7,FAR?,NG7,P1, IKN0Z,A8.A9.P7R» 
lT8,»t3,P3,SB,T9,H9,P9,S5,TS«,TS9,PS8,PS9,VB*¥SS,AMe,AK9,IC0H) 
1MSK«:=!CCN 
SO  TO  {4,4,4, 21, ICON 
CALL  ERROR 
T9=*8 
H9=KC 
?9=P8 
39*33 
TS9=TS8 
PS9*PC8 
79=98 
AK^-AH8 
A9=58 

IMSHOC  3 1 Cffii+3 
SMR{&}=(P7R-??i/P7R 

IFHHNGZ.eO.l?  HRITE16.100)  A8fAH3,49,AM9 

FORMAT (iAMCMOZZLS  DESIGN, J0X3H  A8=,£I5.8,3H  4M8‘,E15»3, 

18H  A9=.E15.3,e«  AK9=,E15.8» 

RETURN 

END 


g  ip, 

ill 


1.1 


m 

•sy 
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AFAPL-TR-87  - 12  5 
Part  D 


AISFTC  PERF  DECK » M94/2 » XR7 

SUBROUTINE  P£RF 


i  i 


1HORO 

.IDES 

» JOES  tKOES 

»NOOE 

tINIT 

*  IOUNP 

.  IAHTP 

2I6A5HX 

.IOSURK 

.IAFT3N.I0CD 

« IHCO 

. IDSKOC 

iIMSHOC 

.N02FLT 

3I7RYS 

.COOPER « NOKAP  rWJHHAP 

tHAPEDG 

tTOLALJ.  »F.RR  C6) 

CORN ON 

/DESIGN/ 

1PCNPGU 

ifCHCSU.TASU  ,DUR01 

.DUND2 

.DELFG 

-D5LFN 

.OELSFC 

2ZF0S 

.PCNFDS 

,PRFOS  pETAFOS 

tNAFDS 

tPRFCF 

tETAFCF 

.wafcf 

32CD3 

.PCNCCSi PRC55  jETACDS 

rWACDS 

,PRCCF 

»ETACCF 

pHACCF 

4T4DS 

t NEEDS 

tDTCODSoETAQOS 

; VA3CDS 

jOPCOOS 

,OTCOCF 

»ETABCF 

57FKP0S 

»CNHPDS»6THPDS,TrHPCF 

tCHHPCF 

»  ETHPCr 

.DKHPCF 

,T2DS 

6TFI.FOS 

tCNLF^S ,ETLPOS, TFLPCF 

tCNLPCf 

.ETLPCF.DHLPCF 

.T21DS 

7T240S 

tMFDOS 

.DTDUD3.ETAD0S 

.WA23DS 

jDPDUOS 

.DTOUCF 

.ETADCF 

8T70S 

.WFA05 

tOTAFDS, ETAADS 

,WG6CDS 

*DPAFGS 

*D7AFCF 

.ETAACF 

9A55 

.A25 

»A 6  ,A7 

»A8 

»A9 

,A2S 

.  A29 

APS55 

,AH55 

iCYuNOZtCVRNOZ 

.ASSAY 

*A9SAV 

,A2 SSAV 

.A29SAV 

COMMON 

/  FRONT/ 

1T1 

.Pi 

»H1  *S1 

.72 

,P2 

»H2 

tS2 

2T2X 

,?21 

,K21  ,S21 

*73 

»P3 

*K3 

.S3 

374 

♦  P4 

,H4  »S4 

»  T5 

»P5 

*H5 

,S5 

4755 

»P55 

*H55  ,SS5 

»Bi.r 

»BLC 

»BLOU 

,3L03 

5C«P  . 

,PRF 

tETAF  .WAFC 

,WAr 

»WA3 

»MG4 

.FAR4 

6CKC 

*PRC 

,ETAC  »WACC 

*WAC 

.ETAS 

tDPCOH 

i  DUMP 

7CNK? 

.etatmp 

»DHTCHP,CHTC 

»8LKP 

♦WG5 

♦  FAR5 

*cs 

8CK.P 

.etatlp 

>CHTCLPrDHTf 

.BLLP 

»*S5S 

• FAR55 

,  KP  E  XT 

9AN 

.ALTP 

»E7AR  »2F 

.  PCNF 

rZC 

tPCNC 

•  WFB 

A7FFHP 

»T??LP 

*PC9Lr  iPCBLC 

.FC5L0U 

,PC8L03 

•PCBS.H? 

.PCBLLP 

CGMSG8 

/  SIDE/ 

XXPI 

.XVAF 

tXVAC  ,ASLF 

sXCLtWi 

»XH3 

* CUNS1 

.0UHS2 

XX721 

,XP2I 

?XH21  iXS21 

*723 

,P2  3 

*H23 

.S23 

3T24 

»F24 

>H24  ,S24 

»T25 

,P25 

*H25 

.S25 

4T28 

t?28 

,K23  .523 

*729 

,F29 

»H29 

.S29 

SVAO 

»WFD 

*XG24  .FAR24 

»ET*0 

.DPDUC 

.BYPASS 

.DUNS2 

6TS28 

,PS28 

*V28  ,AN28 

*TS29 

*PS29 

»V2  9 

.AN29 

COfcKQK 

/  2AOC/ 

XX755 

»XP53 

sXH55  ,XS55 

.  XT25 

,XP25 

*XK25 

.XS25 

WWF9 

.XWG55 

tZFARSS.XyFO 

*XWG24 

,XFAR24 

.XXPI 

.OURB 

3T6 

»?S 

»H4  ,S& 

*77 

»p7 

»H? 

fS7 

478 

»PS 

*N3  S3 

»T9 

»P9 

,H9 

.S9 

5WG6 

»  UFA 

,V67  « FAR 7 

tETAA 

.DPAFT 

»V55 

»V25 

6PS4 

,V& 

.&H6  tTS » 

*PS7 

*V7 

,AN7 

,AM25 

7TS3 

*?S8 

.as  .Ate 

,TS9 

pPS9 

»V9 

.AK5 

8 l'i 

*FRD 

tVJO  *fg*q 

»VJR 

jFGHM 

»PGPD 

.FGPM 

°FSS 

.FO? 

*»FT  »W6T 

tFART 

»fc 

.FN 

.SFC 

DATA  AMORO/&H 

?£8F/ 

HORG-Avoao 

6*32*174049 

C£PSF«2i 14,2X70 

NfT**FO*SF£HWfA 

NAT««ftS-2L*ia 

yGr*w.«T*a?7 

far7*wf>/na? 

¥A»AH*CS 

FR&*¥A*VA?rS 

VKK0£=SV9 

t--6FN«Ca?SP?  f  ?S9-P  1 1 »  A9 
iF!IS3j«Xi.ST.0>  GO  TO  I 
YJ0«CV9*tQZ*Y29 
FM0eyj0i«524/G 
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A  FA  PL- TR -6? -125 
Part  II 


FGPC*CAPSF->{PS29-?1}*A29 
1  FGN=FGHH+FGHD 
FGP=FGPH+FGPD 
FG=FGH+FGP 
FN*FG-FRD 
SFCe3600**WFT/FN 
FG=DElFS*FG 
FN=DELFN*FN 
SFC=DELSFC*SFC 
CALL  OUTPUT 
CALL  ERROR 
RETURN 
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MMU 


AFAPL-TR-57-125 
Part  II 


1 


AFAPL-TR-67-125 
Part  II 


I  I 


WRITE(6,lO2)W0RD,AH*ALTP,T4,T24,6TAR 
CALL  CONOUT (2) 

WRITE (6 , 194) (W { IMSHOCf I ),I=1,4) ,FG»FN,SFC 
IF ( IGASHX.GT.O)  GO  TO  5 
W'lITE  (6  •  105  )  (W  { I  OSHOCt  I  1,1  =  1, 4) 

WRITE {6,1061 LOOPER 
W0R0=AW0RD2 

HR I T£ ( 6 , 107 1 WORD , 2F , PCNF , ZC, PCNC , T4 . MODE 
WRITE (6,109) 

WRITE (6, 109) (ANSI (I), I =1,80) 

WRITEI6 , 108 ) 

WRITE (6, 109) {ANS2II )» 1=1,80) 

WRITE  C6 ,108 ) 

WRITE  (6, 109)  (ANSMI)  ,1=1*48) 

WRITE(6,108) 

WRITE (6,109) (ANS4( I  1,1  =  1,72) 

IF(IOES.EO.l)  GO  TO  6 
A8=A8SAV 
A9=A9SAV 
A28=A28SAV 
A?°=A29SAV 

IF (  I0UMP.NE.2)  GO  TO  6 
WRITE (6 , 110) 

CALL  SYG ( 2 ) 

CALL  ENGBAL 
RETURN 

FORMAT (1HB,A6,14X7H  AH=,F7.5,6X7H  ALTP=,F7.0, 

16X7H  T4=,F8.2,25X7H  ETAR=,F7.4) 

FORMAT ( 1HB, A6 , 14X7H  AH=, F7.3,6X7H  ALTP=,F7.0, 

16X7H  T4=,F8c2, 5X7H  T7=,F8.2,5X7H  ETAR=,F7.4) 

FORMAT  (1HB,A6,14X7H  AB=,F7.3,6X7K  ALTP*=,F7.0, 

16X7H  T4=,F8*  2,5X?H  T24=,F8.2,5X7H  ETAR  =  ,F7.4i 

FORMAT (6H0MA IN  4Afc, 9X3HFG= , F9. 2, 18X3HFN=,F9.2 , 10X4HSFC= ,F8.5) 
FORMAT (6H  OUCT  ,4A6) 


FOP-MAT  (16H1CONVERGEO  AFTER,  1 4, 6H  LOOPS, /,1H81 
FORMAT C 1H  ,A6,9X,5E15.6, 14 J 
FORMAT ( 1H  ) 

FOP-MAT  (1H  ,8E15.6) 

FORMAT (1H1) 

END 


AFAPL-TR-37-125 

PartC 


3IBFTC  CON  OUT  DECX.K9A/2.XR7 
SUBROUTINE  CONOUTUCON) 

COMKON  /  ALL/ 

IWOXO  .IDES  j^DES  fXOES  ,«OCE  . INIT  ,IDGNP 

2IGAS*X,ID8U&N,UFT8N,IDC0  ,IMCO  , ICSHOC , IKSKOC 
3ITRYS  .LGOPER.NOftAP  .KUMHAP. NAPEDG.TOLALl  .ERA  (6) 
COMMON  /DESIGN/ 

1PCMFGU.PCHCGU.T4GU  ,DURD1  ,OUH)2  >DELFG  .OELFH 
2ZF0S  .PCNFDS.PRFOS  .ETAFOS.WAFDS  ,PRFCF  .ETAPCF 

3ZCDS  .PCWCDS.PRCDS  ,ETACDS,*ACDS  .PRCCF  .STACCF 

4T40S  tWrBDS  . OTCODS  . ETABBS  »  WA3CDS » DPCOOS * OTCOCF 
5TFKP0S.CNKPDS»ETHPDS»T?HPCF,CNHPCF,ETH?CF,0HHPCF 
6TFLP0S  rCNLPDS t  ETLPDS »  TFLPCF ,CNL  PCF , ETl PCF .DHL  PCF 
7T240S  ,*TDOS  rOTOUOS»ETAOOS,WA23DS.DPDUDS,DTDUCF 
6T70S  .WFAOS  .OTAFDS.cTAADS.lfCtCCS.GrAPnS.uUFCf 
9A53  ,A23  »  A6  ,A7  ,AS  ,A9  .A28 

APS55  ,A*5S  .CV0N3Z  .CVNNGZ .AESAV  .A9SAV  ,A23SAV 

COTTON  /  FRONT/ 


1PCMFGU.PCHCGU.T4GU 
2ZF0S  .PCNFDS.PRFOS 
3ZCDS  .PCHCDS, PACOS 


. INIT  .IDUNP  »1ANTP  , 
* ICSHOC , I MS HOC .NOZFLT , 


OEISFC 

KAFCF 

WACCF 

ETAECF 

T2DS 

T21D5 

ETADCF 

ETAACF 

A29 

A29SAV 

52 

53 
S5 

8LG3 

FAR4 

DUHF 

CS 

HPSXT 

,MFB 

PC8LLP 

DUHS2 

iS23 

S25 

,S29 

DUK33 

,AM29 


XXT55  , XP5S  ,XH55  ,XSS5  ,XT2S  ,XP25  ,XH25  ,XS25  . 

XXaFB  .XRG55  ,XFAR55,XHfO  .XMG2A  .XFAP.2A.XXP1  .DUMB  , 

3TF  »P6  ,H6  »S6  ,77  TP7  »H?  *S?  » 

ATS  ,P8  .H8  « S9  ,79  tP9  ,H9  ,S9  , 

5WG4  ,HFA  ,WG7  ,FAR7  ,ETAA  ,OPAFT  ,V55  ,V25  , 

6PS6  »V6  ,AM6  ,TS7  ,PS7  ,Y7  ,AN7  ,A*25  , 

7TS8  ,?S6  ,V8  ,AK8  ,TS9  ,PS9  ,V9  ,AN9  , 

8VA  .FAD  ,VJO  ,FSMO  ,¥JH  ,FGW  ,FGPO  ,FGPK  , 

9FW  »FGP  ,HFT  ,WGT  .PART  ,FG  ,FN  ,SFC 

DIMENSION  PARAMt 280 ) , KOROT (2S0 ) , 1 OUT ( 1031 ,AOUT (6 J t WOUT (6) 
EQUIVALENCE  {PARAH.PCHFGU} 

DATA  tKORDYU!,  1-1,981/ 

lfcHPCNFGU tfcHPCNCGUt 6HT4SU  .6H0UK01  .6H0UMD2  .6HDELFG  ,6*0ELFN  , 
2&H0EL SFC , 6HZFDS  .6HPCKFDS.6HPAF03  , 6HE TAFDS » 6HWAF0S  .6HPRFCF  , 
36HETAFCF *6KV<AFCF  .6HZCDS  ,6HPCNCDS,6HPRC0$  *6WETACDS,««ACDS  , 
4&HPRCCF  ,6HETACCF,6H*ACCF  .6HT4BS  .fcHVFBDS  »6H0TC00S.6HSTAB0S, 
5&WA3C0S  »  6H0PC0CS  »  6HDTC0CF  .-6KETA8CF  ,  6HTFHP0S  ,6KCNHPDS  »6HETHPDS  , 
66HTFHPCF ,6HCNH?CF,6HET«>CF,6H0HHPCF .6HT2DS  .6HTFLPDS  t6HCKLPDS* 
76reTLPOS*&MTFLPCF,6HCNLPCF,6«TLPCF,6HDHlPCF,6KT2IOS  .&HT24DS  , 
86HNF00S  ,6H0T0yDS»6HeTAD0S*6H»A23DS,6H0PDU0S,6K>T0UCF,6HETA0CF* 
96KT70S  ,6HiiFA0S  ,6«0TAFDS» 6HE7AAOS .6KVG6CDS *6«>PA>  DSt6K0TAFCF  » 
IAKcTAACF ,6HA55  ,6MA25  ,6NA6  ,6HA7  ,6fU8  »6NA9  , 

2BHA28  ,6HA29  ?6HPSS5  ,6HA»5S  ,6rtCVONOIr6«CvaNOZ,6HA8SAV  , 
3AHA9SAV  t6H52SSAV,NHA29SAV,6HTl  ,6H?l  r6HNI  t6HSl  » 

4AMT2  »6MP2  »6!!a2  ,6HS2  .6HT21  ,<>HPZl  ,6KH21  , 


ITI 

»P1 

.HI 

,S1 

.T2 

«P2 

»H2 

2T21 

,P21 

.H21 

>  S2i 

.73 

,P3 

»H3 

3TA 

»  PA 

.HA 

.SA 

»T5 

»P5 

,H5 

4T5S 

*P55 

,H55 

*S55 

.BLF 

»BLC 

•BLOW 

5CKF 

tPRF 

»ETAF 

•KAFC 

.RAF 

.rfA3 

.VGA 

6CNC 

.PRC 

.ETAC 

.  WACC 

,WAC 

»ETA8 

.DPCOM 

-:»»? 

.ETATWP 

.OHTCHP.DHTC 

»BLKP 

.WG5 

,FAR5 

8CNLP 

.ETATLP 

,£«TCIP 

,OHTF 

,SLLP 

,WG55 

.FAR55 

9AH 

.ALT® 

.ETAR 

♦ZF 

.PCNF 

»zc 

.PCHC 

ATFrHP 

COMMON 

.TFFLP  .PCBLF 
i  SIDE/ 

, PCBLC 

.PCBLOU.FCBLOB 

♦PCBLNP 

XXP1 

.  XVAF 

•  XMAC 

,XELF 

. XBLOU 

,XH3 

tOUHSI 

XXT21 

tXP21 

,XH2i 

,XS21 

»T23 

•  P23 

,H23 

3T2A 

,P24 

.M2A 

,S2A 

,T25 

.P25 

.H25 

4726 

.P28 

.H2S 

,S28 

.T29 

,P29 

,H29 

5KAD 

.WFO 

,HG2A 

.FAR24 

,ETAO 

.OPOUC 

.BYPASS 

&TS28 

C0W3N 

,PS28  ,V28 

/  BACK/ 

.AH28 

.TS29 

,PS29 

.V2v 

XXT55 

,  XP5S 

,XH55 

,XSS5 

,XT2S 

.XP25 

,XH2S 

XX«FB 

tXMGSS 

t  XFAR55 

,xafo 

.XMG2A 

.XFAP.2A.XXPI 

3TF 

»P6 

,H6 

»S6 

.77 

tP7 

»HY 

ATS 

»P8 

»H8 

« S9 

.79 

»P9 

»H9 

BUGS 

,WFA 

,WG7 

,FAR7 

,ETAA 

,  OP  A  FT 

.  V55 

6PS6 

»V6 

*  AN6 

»TS7 

,P$7 

»Y7 

,AH? 

7TS8 

»?S6 

.vs 

.AK8 

,TS9 

,PS9 

»V9 

8VA 

.FRO 

,VJO 

,FSMO 

»¥JN 

,FGW 

,FGPO 

9FW 

,FGP 

,HFT 

»W6T 

.FART 

»FG 

»FN 

»6HNI 

»<>HPZl 


•  9HS1 
.6KH21 


l|f/V V  ViidR 


-v 


AFAfL-TK-87-  125 
Part  B 


56HS21 

*&MT3 

,6HP3 

,6HH3 

,6HS3 

,6HTA 

,6HPA 

DATA  (WORDY ( I ) , I 

=99,169)/ 

6SHHA 

,6H$4 

,6HT5 

,6HP5 

,6hhs 

,  6HS5 

.6HT55 

TS8P55 

,6HHS5 

,6HS55 

«6H8LF 

,6N3LC 

,6HSL0U 

,6K3L0B 

86HCNF 

,t>HP5iF 

,6H£TAF 

,6HWAFC 

j6HWAF 

,6HVA3 

,6HWG* 

96HFAR4 

S6HCNC 

,6HPRC 

,cHETAC 

, 6KWACC 

,6HWAC 

,6HET*B 

liHorcaw 

,6HD*JHF 

-6KC.VHP 

,6HETATHP 

, 6H0HTCHP 

,6H0HTC 

,6HBLKP 

26HWG5 

,6NFAR5 

,6HCS 

, 6HCNLP 

,6HETATLP 

,6H0HTCLP 

vAHOHTF 

36H3LLP 

.6HWG55 

♦6HFAK55 

,4HHP£XT 

,6HAN 

,6HALTP 

,6HETAR 

A6H2F 

,6HPCNF 

,6HZC 

,6HPCNC 

,6hWFB 

,SHTFFHP 

,6HTFFLP 

56HPCBLF 

,5BPC6LC 

,6HPCSL0U 

»6HPC8LOB 

, 6HPCSLHP 

,&HPC3LL? 

,6HXP1 

66HXVAF 

»6HXWAC 

»6HX8LF 

,cHXBLOU 

,6HXH3 

,6H0UHS 1 

.6HOUNS2 

76HXT21 

,6H>'21 

,6hXH21 

.6KXS21 

,6HT23 

,6HP23 

,6KH23 

86HS23 

,6K124 

,6HP2A 

,6HH24 

,6HS24 

»6HT25 

,6KP25 

9&KH25 

,6KS25 

.6HT28 

,6HP28 

,6IW28 

•6HS28 

,6MT25 

OATA  (WORDY (I), I 

-190,280)/ 

I6HP29 

,6ttN  29 

,6HS29 

,6MWjO 

, 6HWFD 

,6KWG2A 

,6H?AR2A 

26HETA0 

,6HDP0Uf. 

,6N8YPrSS 

, 6HCUNS3 

,6HTS2S 

,6HPS28 

,6HV28 

36HAH28 

-6HTS29 

,6HPS2V 

,6HV29 

, 6HAK29 

.6HXT55 

,6HXP55 

46HXH55 

,6HXSSS 

,6HXT25 

, 6HXP25 

.6HXH25 

,6HXS25 

,6HXWFB 

56HXWG55 

,6HXFAS5S,6HXWFD 

,6HXWG2A 

,6HXFAR2A 

♦6HXXP1 

,6?4DU*S 

66HT6 

,6HP6 

»6HH6 

,6KS6 

,6HT7 

,6HP7 

,GHH7 

76HS7 

,6HT3 

,6H5>8 

,6MKg 

.6KS8 

,6HT9 

-6HP9 

8&HH9 

,6HS9 

,  fcHWGfc 

, 6HNFA 

,6KWS7 

,6HFAR7 

,6HETAA 

96H0P-FT 

.6HV55 

,6HV25 

,6HPS6 

,6HV& 

,8HAN6 

,6HTS7 

16HPS7 

,6KV7 

,6MAN7 

.6HAK25 

»6HTS3 

,6HPS8 

»&HV8 

26HAH8 

,6HTS9 

T6KP59 

,6HV9 

,6HAH9 

, &HYA 

?6KFR0 

36HVJ0 

toHFGHO 

,6HVJH 

.6HFGP0 

,6KFGPX 

.fcHFGK 

AfrrfPGP 

v6KWrT 

,6HWGT 

,6HFART 

,6HFG 

,6HFN 

,6HSFC 

OAlf  A  TKEENO,  BLANK,  1.1  HIT/BHTHEEN0,6H 
GO  TO  11, 12), ICON 
C  ***  INPUT  SECTION 

1  00  4  H*i,102 
NUK*N 

REAQ(S,I00)AINTCHANGE 
IF(AIN.EQ.THEEND)  SO  TO  5 
00  2  J*1,UKIT 

JJeJ 

3F(AlN.eQ.«0R0Y(  J)J  GO  TO  3 

2  CONTINUE 
WRITE ',6, 101  JAIN 
GO  TO  A- 

3  IQU?»NUKi«JJ 

IF  (CHANGE  .NE -BLANK  }  WOROY  (JJ  l-CHAMGE 

4  CONTINUE 
WRITE (6 ,102) 

5  NUH«NU»-1 
RETURN 

C  ***  OUTPUT  SECTION 

12  IFJKUH.E0.2)  CO  TO  16 
N=NUK 

J*4 

DO  15  I=l,NUHT6 
IF (H.GT .6 )  CO  TO  13 
J=N 

13  N«N~6 

OG  14  K=1*J 

LU+K-1 

M*IOUT(L) 

WQUT (K) =NGROY  CHI 
IA  AQUTIK)*PARAK{H) 


fc'^V  »■ 

gg?;. 


1 


'  -T;  '  Sc---'  .:. 


'.»£*•  ,-<  .«,*««  ,  -ji  ’ 

*C*  **&?  f-?*'  v 
tracts*  --  t.1 »  f 

.AXES  .*201  .!S?7  .IDS?-*?  ,|A??T? 

*  ,  ll.*£^Sf,;i=T8^,  Ilv*  .Jit-w  ,XD4J^)C,lKSHCK:»»CQ2FLr 

M  R'S  »l£X*P£R  .NOMAP  ,J£jn*A?,5>A?£0G,TCtAU  ,SR£I6) 
ftw-tl  /GfSIGM/ 

lPCNF£U,Pv.NCGU,T40*‘  ,CU«01  ,00*02  »D£LFG  -OElFH  »D£LSFC 
22rCS  ,P£NFC£,rAF05  , ETAFOS, WAF0S  .PAFCF  .ETAFCF  .WAFCF 
32COS  .PCNCCS.PRCDS  , ETACC5.WACDS  ,PRCCF  .E7ACCF ,WA££F 
AT405  ,WFBOS  » OTCODS » ETA9DS, WA3CCS ,DPCCDS .DTCOCF » ETABCF 
5. FHPD5»CSHPPS,£THP0S,7F«PCF  .CNHPCF ,£THPCF  .OHHPCF , ?ZDS 
6TFLP0S,GNtP0S#ETLP0S«TFi,PCF,C«.FCr  .ETLPCF  rDMLFCF  ,T2105 
7T24DS  .WFODS  .OTOUDS.ETAODS.WAZBDS.CPOiJOS.OTDIXP.ETADCF 
3T70S  »WFADS  .OTAFDS.ErAADSjHG&CCS.DPAFOS .DTAFCF ,ETAACF 
9A55  tAZ5  »A6  ,A7  »A8  ,A9  ,A23  rA29 

APS55  ,AH55  ,CYONOZ,CVHNOZ,A8SA¥  .ASSAY  .AZ8SAV.A29SAV 
COMMON  /  FRONT/ 


,P2  »«2  «S2 

#P3  ,H3  >53 

,P 5  ,H5  ,55 

•8LC  ,3LDU  ,BL08 

,HA3  ,WG4  ,FAR4 

,ETA8  .DPCQH  ,OOMF 

,HC.5  ,FAR5  ,CS 


e  t  BLOU  fBLOB 

5CNF  yPF.F  ,EiAF  »WAFC  ,WAf  ,WA3  ,WG4  ,FAR4 

6CNC  pPRC  ,ETAC  ,WACC  »WAC  ,ETAB  ,DPC0M  ,OOHF 

7CHHP  ,  ETATHP  ,  OHTCKP ,  OHTC  ,  BLHP  ,HC.5  ,FAR5  ,CS 

6CNXP  .ETATLP.OMTCLP.CHTF  ,BLLP  ,KG55  ,FAR55  ,H?EXT 


9AH  ,ALTP  ,ETAR  fZF  ,PCHF  ,2C  ,PC NC  ,WF8 

ATFFHP  tTFFLP  ,PC3LF  .PCBLC  ,PCBLDO,PCBL03,PCBLHP,PC8LLP 
COMMON  /SIDE/ 

XXPl  ,XHAF  ,XWAC  ,X8L?  rXBLOU  ,XH3  iDUMSI  .DUM5Z 

XST21  ,XP22  »XH2I  ,X$2i  ,T23  ,P23  ,H23  ,523 

3T24  ,P24  ,H24  ,524  ,T25  ,P25  ,H25  ,S25 

4123  ,P28  ,H28  ,S28  ,T25  ,P29  ,H29  ,529 

5HM  ,NFO  j«5?a  >FAR24  ,ETAD  ,0?2UC  ,EYPA5$?0UHS3 
6TS28  ,PS28  ,V28  ,AM25  ,T529  ,P529  ,V29  ,AH29 

COMMON  /  BACK/ 

XXT55  ,XP55  »Xli55  ,XS55  ,XT25  ,XP25  »XH25  ,XSk5 

XXHFB  .XVG55  .XFAR55.XHP0  .XWG24  ,XFAR24,XXP1  .DUMB 
376  ,P&  »H6  ,56  ,T7  ,P7  ,H7  ,5?  , 

4T8  ,P3  ,K8  ,58  ,T9  ,P9  ,H5  .59 


tPCNC 


XXPX  ,XMAF  ,XHAC 

<X?21  ,XP22  ,XH2I 

3T24  ,P24  ,H24 

4T2S  ,P28  *H28 


,XWAC_  ,X8Lr  ,XBLOU  ,XH3  ,0UMS1  .OUMSZ 

*XH2I  ,XS2i  ,T23  ,P23  ,H23  ,523 

,H24  ,524  ,T25  ,P25  ,H25  ,S25 

»H28  , S23  ,T25  ,P29  ,H29  ,529 

♦«G?A  .FAR24  ,ETA0  ,0?3UC  ,EYPA5$?0UHS3 
,V28  »  AH26  ,T529  ,P529  ,V29  ,AH29 


,FAR7  ,ETAA  ,DPAFT  ,V55 

,TS7  ,P57  ,V7  ,AM7 

,AK8  ,TS9  , PS9  ,V9 

,FGKO  ,VJM  ,FGMM  ,FGPD 
,HGT  rfiWT  ,?S  ,FN 


,DUH8 

,S? 

.59 

,V25 

,A*25 

,AM9 

tFGPH 

t$FC 


DIMENSION  TRASH1 {80} •TRASH2(80).TRASH3«48I .TRASH4C72) 

oS^vAIlSo%^S!^i/^>CNFGU,,^TftASH^,T1,,,7RASH3,XP1,t(TRASH4,XT55, 

WRITE (6. 100) WORD 
WORO*AMORO 

WRITE <6, 102}K0RDfZF,PCNF,ZC»PCNC»T4, MODE 
WRITE (6.103) 

WRITEI6.104) (TRASHlfl ) .1*1,90) 

WRITE (6,105) 

WRnE(6,104)  (TRASH2II )  ,1*1,80) 

WRITC.56,103) 

WaiTE(6,i04)!TRASH3tI),I«l,48) 

WRITE (6, 103? 

WRITE (6, 104) (TRA5K4(I ) ,I«1,72) 

HRITE(6,103) 

WRITE  C6,106)L00FER 


Fssss  * 


2 

100 

102 

103 

104 

105 

106 


5F(10URPeE9.O!  50  TO  ? 

WRITE (6, 105 J 
CALL  STG { 2 J 
CALL  ENG8AL 
RETURN 

FORHAT { 28H0AH  ERROR  HAS  BEEN  FOUND  IN  ,A6» 

FORMAT ( 1H0»A6,9X,5E15*6» 14) 

FORMAT (2H0  J 
FORMAT { lH0r6E15.6  J 
FORMAT ilHl) 

FORMAT! 25H0FAI LED  TO  CONVERGE  AFTER, 14, 6H  LOOPS) 


END 
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*I3FTC  SY6  0ECK,H99/2,XR7 
SUBROUTINE  SVG U CON J 
OIM6KSION  WORD! 132) 

DATA  0HE00L/6HS  / 

GO  TO  (1,2),! CON 
1  END  FILE  8 
REWIND  0 
RETURN 


C  TERMINATE  THE  FILE 

2  WRITE C8, 500) 

500  FORMAT { 12HSJ JSSSSSSFSJ ) 

END  FILE  8 

REHIND  S 

C  READ  RECORD 

5  READ (8, 501) (WORD! I ) » 1=1 » 132) 

501  FORMAT (132A1J 

C  CHECK  FOR  12  LEADING  DOLLAR  SIGNS 

00  10  1*1,12 

IFIWORDCI )— ONEDOLJlltlOf lx 

10  CONTINUE 
RETURN 

C  CHECK  FOR  6  TRAILING  DOLLAR  SIGNS 

11  DO  15  1=1,132 

1=1 

IF(WORO(!  )— 0NED0L)I5»12»15 

12  K=l+5 

DO  13  J=I,K 

IF (WORD (J)-ONEOOL) 15,13,15 
12  CONTINUE 

GO  TO  20 
15  CONTINUE 

WRITE (6,502) 

502  FORMAT C1H0, I 2HERR0R  IN  SYG) 

RETURN 

C  PRINT  LINE 

20  1*1-1 

WRITE (6, 501 ) (WORD* ,M*1,I ) 

GO  TO  5 
END 


S3 


•  ierrc  TKCOMP  DECS, S94/2,XR7 

SUBROUTINE  THCOMP (PR, ETA, T, H, S, P,TO,HO ,SG ,PO ) 
PO=P*PR 

TP=T*PR**0. 28572 
DO  1  1=1,25 

CALL  THERH0(P0,HP,TP,SP,X1,0,X2,0) 

OELS=SP— S 

IF(A85!DELS)»LE.0.00005*S)  GO  TO  2 

1  TP=TP/EXP{4.*DELS) 

CALL  ERROR 

2  HO=H+( (HP-HJ/ETA) 

CALL  THERMO! PO» HO, TO, SO, XI ,0,X2« 1 ) 

RETURN 

END 


AFAPUT?.-^'- 1 25 
PsrtH 


AP4PL- n?-67-l2§ 
Part  ?! 


SIBFTC  THERMO  DECK,H94/2,XR7 

SUBROUT I NE  THERMO i PX  , HX , TX , 5X , AMX»  L  ,  F  AR , K ) 
FX=0. 

IF(L.cQ.l)  FX=FAR 
IF(K-EQ.l)  GO  TO  1 

CALL  PROCOH<FX,TX,CS,AK,CP,R,PHI ,HX> 

GO  TO  3 

1  TX=4.*HX 

DO  2  1=1,15 

CALL  PROCOH(FXfTX,CS,AK,CP,R»PKI,H» 
DELH=HX-H 

IFCABSCDELHJ.LE.O.OOGOl^HX)  GO  TO  3 

2  TX=TX+4.*OELH 
WRITE  C8,100J 

100  FORMAT! 3 1H WO  CONVERGENCE  IN  THERKOSSS$$$ ) 

3  5X=PHI— R^ALOGtPX } 

AMX=1.986375/R 

RETURN 

END 


4*M?L-TF-f"-l2? 

Par*  U 


SJSPTC  PROCOM  DECK, M94/2»XR? 

SUBROUTINE  PR0C0M{FARXt7£XfC$EX,AKEX,CPEX,REX,PH! ,HEXi 
IFtFARX.LE. 0.067623)  GO  TO  1 
FARX=Q. 067623 
WRITE(S,101) 

1  IFITEX.GE.300.)  GO  TO  2 
T£X=300. 

WRITE (8,102) 

2  IF (TEX. LE. 4000.)  GO  TO  3 
TEX=4090. 

WRITE(8,103i 

3  IF(FARX.GE.O.Q)  GO  TO  4 
FARX=0.0 

WRITE (8,104) 

C  AIR  PATH 

4  CPA  ={(({{ I 1.0115540E— 25*TEX— 1.4526770E-Ci) *T£X 
l+7,6215767E-18}*TEX-i.5128259E-14)*7EX-6,?178376L-X2) 
2*7EX+6.5519436E-08)*T£X-5.i536879£-05)*TEX+2.502005lE-Sl 

HEA= {(((((( 1 . 264442 5 £-26*T£ X-2 * 0752 522  c-22 ) *T2X 
1+1. 270263OE-i8)^TgX-3. 025651 3E-I5)^TEX-;.6794594P-12)*T£X 
2+2. 1 839826E— 08 )  s=TEX-2.  5768440E-05  i*TE.X+2  .5G2005XE-0 1 }  *TEX 
3-1.7553886E+0D 

S£A=+2.5020051E~Q!3?AL0G  (TEX )  +  ( {  {  < 1.445076 7E-2 6*7 EX 

1- 2.4211288 E-22 i  *7gX+l . 5243153E-18  )*TEX-3.78*0648E-13 ) *TSX 

2- 2.23927'?0E-12)#TEX+3.2759743E-08S*TEX-5.1576879E-05!*7£X 
3+4.5432300E-02 

IF(FARX.LE.O.O)  GO  TO  5 
C  FUEL/AIR  PATH 

CPF  c{(5H  (7,2678710&-25*TEX~!.33356c8E-20J*TEX 
1+1.0212915E— IS  )3=T EX-4*  2051104E— 13)  *rTEX+9. 9686793c- 10) +TEX 
2— 1.3771901c— 063 4TEX+1.2258630E— 03) *T£X+? .381 6638E-02 
H£F- U  i  U  { { 9 «084838£E~26-T EX- 1,90509498-21 ) ^TEX 
l+1.7G2l525E-17)*TEX«8«4i02208S-14)*7EX+2.4921698£~IC}*T£X 
2-4 .59 06332 E-Q7  )<TEX+6.129315GE— 04  ?*TeX+7.3316638E— 02) 

3*TEX+3 .0581 530E+01 

S£F=+?.381663SE-02*AL0G(7EX)+{ ( ) { Hlr03S2670E-25»TEX 
*—2.2226 USE— 21 )*T£X+2.0425826E-17) -TEX-1 .05127?6E-13)*T£X 
2+3.322892SE-10)4T£X-6.3859£05E-07)C7EX+1.225863C£-035*TgX 
3+o .4 63398 E— 01 

5  CPEX= { CPA+F ARX-CPF I / { 1 .+ FARX ) 

H£X= (HEA+F ARX*HEF ) / ( 1 .+FAftX  > 

PH  I  =  ( SEA+FARX-S EF ) / (1 . +FARX ) 

AHa'=28.97-.946186*FARX 
R£X=1 • 986375/ AMW 
AKEX=CP£X/(CP£X— REX ) 

CSEX*SORT ( AK£X*REX*7EX*2503i . 37 ) 

RETURN' 

101  FORMAT t? HO, 63H2NPUT  FUEL-AIR  RATIO  ABOVE  LIMITS,  IT  HAS  BEEN  RESET 
273  0. 0676.23, 6HS6S$$S) 

102  FORMAT ( 1H0, 35 HP RO COM  INPUT  T7HP£RATURE  BELOW  3S0.,6HSSSSSS) 

103  FORMAT | 1H0,36HPR0C0M  INPUT  TEMPERATURE  ABOVE  4000. ,6H*S$$$$ } 

104  FORMAT (1H0,38HPR0C0H  INPUT  FUEL-AIR  RATIO  BELOW  ZER0,6HS$$$$S ) 

END 
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rtfirtfl  O  *•*  >'\  r>  r-  i*  z)  "v 


i-t  s 


teiti  C-r  CK  »  **•»  1 1  >  X9  ? 

-  i*;5  S£AS£H{  *■  » '»< B*C  tO..'  f  *■£•&*  *£  1  ******  s*5®1*  rHT  .  . 

y?**-£&5  jGii  A-H  t.ti&m  t^!  5  -C  i  !  NA**  *  »$»*  I  »£!Sr  -NAH  «N3 ft  ?  1  **-"*  ‘  *‘»  ^ 

7±i  WirFOS  j£W>P£K?t  IK? 

=**  iX  S'SO  5x  HUS?  S£  STf^tgti  to  to  ws 

»e«r  p=i*f?vT  i>SDF'wft? 5G^  BJ'7Kr-£3;  0-C  AND  l.C 
IF  KOT  ltfFyt,  P  KtiST  (SO  A  l  -1* 
s?.*  MCCOS*>QC  OK 
KC0BE*01  A  LO 
NC0?E=O2  A 
r-*C00==0?  cRP-OP- 
NCOt>c=U‘  9  LO 

?fCO0c=?O  3  HI 

KCGD£=a 
C-=0. 

0=0, 

c  ***  FJ?:0  A 
CO  1 
IH=I 

IftA.LT.ASdU  G3  10  Z 

1  CONTINUE 

iF(A«5T «.AX( IHi  )  NC0DE=2 
A=AXtIH) 

GO  TO  3 

2  Sf (IH*ST.i:  GO  TO  3 

«CD0E=1 

IH=2 

A=AXtlJ 

3  l£.  =  IH— 1 
lihh=«q;ihi 

tm=NO(IL> 

c  ***  find  e 

A-AX C 7 L II / ( AX ! IH I-AXI  XL) ) 

PP»P 

IFIP.GE-O- I  GO  TO  6 

8L=8XUL*i  )+PRH*(BX(IH,  l)-Bx(ILtl) I 

BH^BX  { lL*LIMi.)‘*PRK**(BX  C  IK»LIKH>— SX{  IL« LlHL)  ) 

IFtB.GE.3i)  GO  10  A 
N-C0DE=KC03E»10 
8=BL 
GO  TO  5 

A  IFI8.Lfc.vSHi  GO  TO  5 

NCOGE-sSCOOE+20 
&=3M 

5  PP«G. 5 

q;2)=c- 

0(31=0. 

6  BH=?P«(8X(IH,LIf?H)-BX(IH,lH+BXUH,IJ 
BL=PPc(BX  C I L»L1HL »— 8X t ILt 1 1 >+BX( IL*1 ) 

00  7  J=2»LI5W 

JH=J 

IF IBH.LT .8X(iH»J)J  GO  TO  8 

7  CONTtHUc 

8  JL=JH— 1 

00  9  K=2»i-«WL 
KH=K 

IF  (BL.LTaSX!  I  L»X?J  GO  TO  >,0 

9  CONTINUE 

10  KL=KH— 1 

Pfv=(BX(IH,rli)-SH»/(8XJ  IH^JKj-BXnHr  JL>> 
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|  1  I-1  UlJUUy  gll-HllM  MWm.  I 


AVAPL=TH-^5?-12-5 

Part  n 


CH*  CXf  IHf  JD— PR  ♦(CXUH,JH)-CX(IH,JL) } 
DH*  OXUH1JL1-PR  *CDXUH,JH)-DX(IH,JL)) 


PR*(BX{ILfKL)-BU/CBX(IL,KHi-BX(IL,KL)> 
CL*  CXCIL»KL)-PR  *(CX{IL»KH)-CXC IL,KL) ) 
OL*  OX{IL»KL  J—PR  <-COX(IL*KH)-DX(  IL t»L)  ) 


8T*3L+PRH*CBH-BL) 
CT*CL+PRM*(CH-CU 
OT  *DL+PRH*  C  DH—DL ) 


IF (P.GE.O. }  GO  TO  13 
01R*SQRT ( B/BT ) 

ERR«<B-BT)/B 

CALL  AFQUIR CQC 1 ) » PP , ERR ,0. ,25. ,0.001 , DIR ,PT * ICON} 
GO  TO  111, 13,12), ICON 

11  PP*PT 
IFCPP.LT. 0. )  PP*0« 

IFCPP.GT.l.)  PP=1. 

GO  TO  6 

12  NCO0E*7 

.3  B*BT 

C*CT 

D*DT 

RETURN 

END 


SO 


AFAPL-TR -67-125 
Part  II 


SIBFTC  HAPBAC  DECK , M94/2 ,XR7 

SUBROUTINE  MAPBAC(MAPfMAPGO#TFFS»TFFfCNS»CN,PCN,T#MODE» IGOtNUNI 
DATA  WH,WLfWT,WS/6H  H.P.  ,6H  L.P.  ,6H  TFF  .6HSPEED  / 

1  FORMAT ( 1HO ,  A6  » 1 2HTURB I NE  MAP  *A6t4HWAS«f£13.6»10H  AND  N0W«tE13»6» 

) 

2  FORMAT ( 1HO ; A6 1 A6f 22HMAS  ALSO  CHANGED  FROM  »E13*6«5H  TO  »E13.6V 
16H$$S$$$) 

IF(NUM.GT.O)  GO  TO  3 

NUMH=0 

NUML=0 

3  IGO=MAPGO+?*( MAP-1 ) 

GO  TO  { 100»200t 300f400 » 500»600 ) i IGO 
c  ***  HIGH  PRESSURE  TURBINE 
100  TFF=TFF+0« 1*( TFF-TFFS ) 

WRITE(8,1)WH,WT,TFFS,TFF 

RETURN 

200  CN=CN+0.05*(CN-CNS) 

IF(MODE.NE.l)  PCN=>PCN*(CN/CNS) 

IF(MODE.EQ.l)  T  =T  *(CNS/CN)**2 
WRITE (8*1 )WHfHSfCNS»CN 
IFINUMH.GT.2)  GO  TO  210 
NUM=1 

NUMH=NUMH+1 

RETURN 

210  DELCN=CN-CNS 

IFIDELCN.GE.O.)  RETURN 
TFF=TFF*< l.+DELCN/CN) 

WRITE (8» 2) WHfWT  »TFFS»TFF 
RETURN 

300  TFF=TFF+0.1*(TFF-TFFS) 

WRITE(8,1)WH,WT,TFFS,TFF 
GO  TO  200 

C  ***  LOW  PRESSURE  TURBINE 
400  TFF=TFF+0.1*(TFF-TFFS) 

WRITE  18 tllWLtWTtTFFS, TFF 
RETURN 

500  CN=CN+0.05*(CN-CNS) 

PCN=PCN*(CN/CNS) 

WRITE(8fl )WL»WS»CNSfCN 
IFINUML.GT.2)  GO  TO  510 
NUM=1 

NUML=NUML+1 

RETURN 

510  DELCN=CN-CNS 

1FIDELCN.GE.0. )  RETURN 
TFF=TFF*( l.+DELCN/CN) 

WRITE(8»2 ) WL»WT  »TFFS»TFF 
RETURN 

600  TFF=TFF+0.1*(TFF-TFFS> 

WRITE (8 » 1 ) WL  t WT*TFFS»TFF 

GO  TO  400 

END 
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non 


AFA  PL  -TR-67- 125 
Part  U 


SIBFTC  CONVRG  DECK,M94/2,XR7 

SUBROUTINE  CONVRG ITI ,HI , PI ,SI ,FAR *WG*PA , IOES ,AO,PR* 
1 TO , HO , PO , SO , TSO , P SO , VO , AMO * I CON ) 

ICON-1  SUBSONIC*  COMPARE  PI  WITH  PR 
ICON-2  SONIC*  COMPARE  PI  WITH  PR 

ICON-4  ERROR 

AJ-778.26 
CAPSF-2116.217 
G-32. 174049 

CALL  PROCOMtFAR *TI  *XX1»XX2«.'CX3«XX4*PHII  *XX6) 

C  ***  SONIC  CALCULATIONS 
J*0 

TSS«0.833*T1 

1  J-J+l 

CALL  PROCOM (FAR *TSS*CSS* AKS*CP«REXS»PHISS *HSS > 
HSCAL-Hl— C5S**2/ ( 2. *G*AJ ) 

DELHS-HSCAL— HSS 

IF { ABS(OELHS)— 0.0005*HSCAL)4*4*2 

2  TSS-TSS+DELHS/CP 
IF  C  J— 15  )  1 »  1 «  3 

3  ICON-4 
RETURN 

4  IF ( IOES) 12*12*5 

C  ***  I5ENTR0PIC  EXPANSION  CALCULATIONS 

5  J-U 

TS I -T 1*1  PA/PI  1**0.286 

6  J-J+l 

CALL  THERMO(PA*HSI ,7SI ,SSI *XX1 * 1,FAR,0) 

IFIABSJ SSI-SI  >-0.0001*SI  >8.8,7 

7  TSI«TSI/EXP(4.*(SSI-SI J> 

IF(J— 30)6*6*3 

8  VI S-SQRT ( 2. *G*AJ* (HI-HSI ) ) 

IF (VIS-CSS  >9*11  *  11 

C  ***  SUBSONIC  DESIGN.  CALCULATE  AO 

9  VO-VIS 
TSO-TSI 
PSO-PA 

CALL  PROCOM(FAR*TSOiCSO*XX2«XX3*REX*PHISO,HSO) 
RHO-CAPSF*PSO/ (AJ*REX*TSO) 

AO-WG/(RHO*VO) 

AMO-VO/CSO 

PR-PI 

ICON-1 

10  TO-TI 
HO -HI 
PO-PI 
SO-SI 
RETURN 

C  ***  SONIC  DESIGN,  CALCULATE  AO 

11  VO-CSS 
TSO-TSS 

PSO-PI*( TSO/TI >**(AKS/(AXS— 1. I > 
RHO-CAPSF*PSO/(AJ*REXS*TSO J 
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AFAPL-TR-67-123 
Part  II 


AO=HG/(RHO*VOJ 
AM0=1«0 
PR=PI 
ICON =2 
GO  TO  10 

k 

C  ***  NON-DESIGN,  CALCULATE  CRITICAL  CONDITIONS 

12  VO=CSS 
TSO=TSS 
PSO=PA 

RHO=CAPSF*PSO/{AJ*REXS*TSO) 

A0CRIT=WG/{RH0*V01 

AMO=1.0 

PR=PSO*( TI /TSO ) ** { AKS/ { AKS-i . ) ) 

I F  < AQ-A0CR1T 113,13,14 

C  ***  NON-DESIGN,  CRITICAL  AND  SUPERCRITICAL  CONDITIONS 

13  ?SO=PSO*AOCRIT/AG 
PR=PR*AOCRIT/AO 
ICON =2 

GO  TO  10 

C  m  NON-DESIGN,  SUBSONIC  CALCULATIONS 

14  PSO=PA 
JxO 

TSO*0.833*TSO 

15  J=J+1 

CALL  PROCOH{FAR»7SO,CSO»AKO»CP,REX,PHISO,HSOi 
RHO=CAPSF*PSO/tAJ*REX*TSO) 

VO=HG/(RHO*AOI 
HSCAL=HI-VQ**2/  i 2**G*A o  ) 

D5LHS-HSCAL-HSQ 

IF(A8S(DELHS)-0,0005*HSCAL) 17,17,16 

16  7S0*TS04DELKS/CP 
IF{J— 15515,15,3 

17  AMQ*VO/CSO 
PR5=PSG*(Tl/TSO)«CAKO/(AKO-1.5) 

IC0N=1 

GO  TO  10 


AFAPL-TR-67-125 

PjutH 


B1BFTC  COHOIV  0ECK.M94/2.XR7 

SUBROUTINE  CONOlVfTI.HI.PI.SI.FAR.WG.PAtlDEStAT.AOtPIR. 
iTT.HT.PT.ST.TO.HQ.PO.SO.TST.TSO.PST.PSO.VT.VO.AMT.AMO.ICQfn 
C  ICOM-1  SUBSONIC.  COMPARE  PIR  WITH  PI 

C  ICON-2  SONIC.  SHOCK  INSIDE  NOZZLE.  COMPARE  PIR  WITH  PI 

C  ICON-3  SONIC*  SHOCK  OUTSIDE  NOZZLE.  COMPARE  PIR  WITH  PI 

C  SCON-4  ERROR7 

DIMENSION  6(91 
Q(2)-0. 

Q(3)*C. 

A4-77A.26 
~  CAPSP-21 16*2170 
0-32*174049 

CALL  PR0C0MfFAR,TI.XXI.XX2,XX3.XX4,PKII,XX6> 

C  ***  SONIC  CALCULATIONS 
4*0 

TS5»0.S33*TI 

1  4-4*1 

CALL  PROCOMfFAR.TSS.CSS.AK.CP.R5KS.PHISS.HSS) 

HSCAL-HI— C$S**2/ ( 2.  «G*A4 ) 

OELHS-HSCAL-HSS 

lF«ABSfDELHSJ~0o0005»H5CAL)4,4,2 

2  TSS«TSS*DELHS/CP 
IFf 4*15)1. 1*3 

3  ICON-4 
RETURN 

4  IF (IDES) 11 .11.5 

C  ***  SONIC  DESIGN.  CALCULATE  AT 

5  VT-CSS 
TST-TSS 

PST-PI*! TST/T1 )•*' AK/f AK-X. H 
. RHG«CAPSF*PST/ f A4*REXS*TST ) 

AT-MG/fRHO-VT) 

AMT-1.0 

C  •**  IDEAL  EXPANSION  DESIGN.  CALCULATE  AO 

PSQ-PA 

4-0 

TS0-TI*IP?0/PI )»*.2B6 

6  4-441 

CALL  PROCOHCFAR.TSO.CSO.AK.CP.REX.PKISC.HSO) 

PHICAL-PH 1 1 -REX*ALOC  f PI/PSO ) 

GELPHI-PHIC&L— PHISO 

IFfABSf DELPHI J-Q*0001»PHICAL)8»8»7 

7  TSO«TSC*rXP(4.AOELPHI) 

IF<4*15)6,6»3 

8  VO-SORTf  2. -G»A4*f HI-MSO) ) 

AMO-VO/CSO 

AO«fAT/AMO)»f2.*(l.*IAK-l.)*AMO**2/2.l/fAK-l.))-*liAK*l.)/l2. 
-  If AK— 1. m 
PIR-PI 
ICON-3 
TO-TI 
HO -HI 
P0-P2 
SO-  SI 
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AFAPL-7R-87-125 
Part  n 


10  TT-TI 
HI -HI 
PT-PI 
ST-SI 
RETURN 

C  ***  ASSUME  SONIC  THROAT  A NO  ISENTROPIC  EXPANSION  TO  AO 

11  VT-CSS 
AKT-t.O 
TST-TSS 
RHO«MG/(AT*VT) 

P5T-RH0*AJ*REXS*TST/CAPSF 

P 1R-PST* i TI /TST 1 ** { AX/ ( AX-1 . I 1 
IF (PST-PA) 12,24,24 

12  TSG-G.95»TI 
NAN-0 

13  CALL  PR0CGK(FA*,TSQ,CSO,AX,CP,REX,PHISO,H5O} 

AKO-SQRT( 2.*( (TI/TS01-1. l/( AX-1. ) I 

A0CAL-IAT/AM0)*(2**(1.‘>(AX— l.l*AM0**2/2»>/ (AX*1.11»*(( AK*1.)/ 
1(2.«(AK-1.111 
£A« (AO- AO CALI /AO 
D!R«SQAT( AC/AOCAl ) 

CALL  AFQUJR(Q( I) ,TS9,EA,0«,  100., 0.0001  ,0IR,T50T,JC0N) 

CO  TO  (14,15, 31, JCON 

14  TSO-TSOT 

IF (TSO—TI 1140,13, 141 

140  TSC«2.*TI/(AX4l.) 

IF (TS0.GT.TSC1  CO  TO  142 

141  TS0»0.9S*T1 
GO  TO  13 

142  IF(Q (2 1.LT.3O.O.GR.AMO.LT.0 .95.0R.MAM.EQ.il  GO  TO  13 
TSQ-2.*TI/ (2.+0.9I *(AX— 1. ) 1 
NAM-1 
GO  TO  13 

15  PS0«PIR*(TS0/T11*4(AX/(AK-1.I1 
IF (PS0-PA) 17,16,24 

C  **«  CRITICAL  FLO*,  ISENTROPIC  EXPANSION  TO  PA 

16  V0-AM0*CS0 
ICON-1 
GO  TO  9 

C  ***  SUBSONIC  FLOW 

17  PSO-PA 
Q(21*0. 

Q(3l-0. 

-J»0 

TSO-0.833«TI 

18  J*J*I 

CALL  PRQC0N(FAR,T50,CS0,AK,CP,REX,PHI$0,HS0} 

RHO»CAPSF*PSO/(  AJAREXRTSO) 

V0-MC/(RHQ»AQ) 

HSCAL-HI— VQ*»2/(2.*G»AJ) 

OELHS-HSCAL-NSO 

IF (ABS(DELHS)—0.0005*HSCAL 120,20, 19 

19  TSO-TSO-KJELHS/CP 
IF(J-15)1S,1*,3 

20  ANO-VO/CSO 
PIR-PS0*(TI/TS0)**(AK/{AK-1.11 


«rflK»*TB»OTB3 

Partn 


TST-TSO 

21  CALL  PRa00M(FAP.,T$T,CST,AK,CP,REX,PKI'7.HST) 
PST«PIR*(7ST/TI)**( AK/(AK~1 . 1 ) 

RH0-P574CAPSF/  {  AJ*REX*T5T » 

VT*H 6/ C RHG*AT  J 
HSCAL»tiI-VT*»2/<2i,*G*AJ) 

EH- C  HSCAL-HST ) /KSCAL 
DIR-1.-MHSCAL-HSTI/(CP*TST) 

CALL  AFCU1R(«(1J, TST,£H,0,, 20., 0.0003, DIR, ISTT.JCON) 

CO  TO  (22,23,3) »JCOH 

22  TST-TSTT 
CO  TO  21 

23  AHT-VT/CST 
ICON-1 

CO  TO  9 

C  **♦  SUPERCRITICAL  FLOW,  ISEHTRQPIC  EXPANSION  TO  PA 

24  PSO-PA 
J-0 

TS0-TI*(PS0/PIR)»*,2«« 

CALL  PROCOK(FAR,TSO,CSO,AX,CP,REX,PH!SO,HSO) 

PKICAL-PHII— REX*ALQG(PJR/PSO) 
v  DELPHI -PH ICAL-PHI SO 

IF(A8SIDELPHIJ-0.0001*PH1CAL)27,27?26 

26  TSO-TSO-EXP (4.0-DELPHI J 
1F(J-1S)25»25,3 

27  V0-S0RT(2.*&*AJ*(HI~HS0)> 

AHO-VO/CSO 

AOIO-(ATM>«I#(2.»t  1.+IAK-1.  J*AH0**2/2.I/IAK+'i .  »**  l  (  AX*1.  J  / 
1(2.«(AK-1.}» 

ICON -3 
N-0 

IF{AO-AO!0)2S,9,29 

C  »«*  SUPERCRITICAL  FLOW,  ISEWTROPIC  EXPANSION  TO  AO 

2*  N-l 

29  TS0-O.t33*TI 
J-0 

30  J-J^l 

CALL  PROCON(FAR,TSO,CSO,AK,CP,REX,PH1S9,HSOI 
ANO-SQRT(2.*( I TI/TSOI-l. J/ ( AK-1. ) 1 

AOCAL«(AT/AHOI*(2.*C2.+(AK-l.)*AHO»»2/2.)/5AK*l.JI—CCAK+l.I/ 

l(2.*(Ar-l,U) 

oela-ag-aocal 

IF(ASS{DELA!-O.O0Ol*AO)3Z>32,31 

31  TSO-TSO*SORTIAOCAL/AO) 

IF  {>*50130,30,3 

32  IF(N 134,34*33 


C  *<**  ONDEREXPANOED,  SHOCK  OUTSIDE  NOZZLE 

33  PS0-PIR*(TS0/7I|**{AIC/{AX— l.II 
VO-ANO*CSO 

CO  70  9 

C  ««•  OVEREXPANDED,  FIND  SHOCK  POSITION 

34  PSX-PIR-J 7S0/TI )♦*( AK/t AK-1.) } 
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PSY=PSX*(2.*AK*AMO**2/UK+l» )-{AK-l. )/{ AK+1*) ) 
IF(PA-PSY)35,36,36 

C  ***  OVEREXPANDED »  SHOCK  OUTSIDE  NOZZLE 

35  ?SO=PSX 
VO=AMO*CSO 
GO  TO  9 

C  ***  OVEREXPANDED,  SHOCK  INSIDE  NOZZLE 

36  PSO*PA 
J=0 

TS0»0.833*TI 

37  J*J+1 

CALL  PR0C0H(FAR,7S0,CS0,AK,CP,REX,PHISD,HS0) 

RHO=CAPSF*PSO/CAJ*REX*TSOJ 

VO=HG/ ( RHC*AG } 

HSCAL=HI— yO**2/(2«*G*AJJ 
DELHS=HSCAL— HSO 

IF (ABS  COELKS)— 0.0Q05*KSCAL)39 ,39,38 

38  TSO=TSO+DELHS/CP 
IF(J-X5J37,37,3 

39  AHO=VO/CSO 
TO«TI 
HG=HI 

PO=PSO*C TO/TSOJ**! AK/  C  AK-J  •  I ) 
SO*PHII-REX*ALOGCPO) 

ICON*2 
GO  TO  10 
END 
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BISFTC  AFQUIR  0£CK,H94/2»XR7 

SUBROUTINE  AFGUIRIX,AIND, DEPEND, ANS,AJ,T0L,D1R,AN£V, ICON) 
DIMENSION  X(9) 

C  XII 1-NAME  OF  ARRAY  70  USE 
C  AINO-INOEPENOANT  VARIABLE 
C  DEPEND-  DEPENDANT  VARIABLE 
C  ANS-ANSWEft  UPON  WHICH  70  CONVERGE 
C  AJ-MAX  NUMBER  OF  TRYS 
C  TOL^EACENT  TOLERANCE  FOR  CONVERGENCE 
C  DIR-0IREC7IQN  AND  PERCENTAGE  FOR  FIRST  GUESS 
C  ANEW-CALCULATED  VALUE  OF  NEXT  TRY  AT  2N0EPEM0AKT  VARIABLE 
C  ICON-CONTROL  -1  GO  THRU  LOOP  AGAIN 
C  >2  YOU  HAVE  REACHED  THE  ANSWER 

C  *3  COUNTER  HAS  HIT  LIMITS 

C  XI2I-C0UNTER  STORAGE 
C  Xm-CHOOSES  METHOD  OF  CONVERGENCE 
C  AIM-THIRD  DEPEND  VAR 
C  XJSJ*7iURp  I MO  VAR 
£  XIA)— SECOND  OEPENO  VAX. 

C  X 171 -SECOND  SNO  VAR 
C  XICf-FIRST  DEPEND  VAR 
C  X|4f-F1*ST  I NO  VAR 

C  XI3)  MUST  BE  ZERO  UPON  FIRST  ENTRY  TO  ROUTINE 

Y-0. 

1F(AN$)1*2»1 
2  D£P-0EPEN0~AXS 

TOLANS-TOL-ANS 
GO  TO  3 

2  DEF-0EPENO 
T CLANS- TOL 

3  IFTA8$(D£P3— T0LANS15*5*4 

4  3F(XI2)-AJ)S,S*7 

5  ANEW* A I NO 
XI2I-0. 

ICON-2 

RETURN 

*  ANEW-Y 

XI2J-XI2JM. 

ICON-1 
RETURN 
7  ANEW-Y 

X 121-0. 

ICON-3 

RETURN 

«  IFIX(3)19«9*12 

C  ***  FIRST  GUESS  USING  DIR 

9  XI31-1. 

XIBl-OEP 
XI91-AIN0 
lEIAfNDf 10*11*10 

10  Y-0IR9AIN0 
GO  TO  « 

U  Y-C1R 

CO  TO  6 

22  IFIXI3J-1. 113*23*16 
C  *»*  LINEAR  GUESS 
13  XI31-2. 

XI63-OEP 

XI71-AIND 

JFIXfll-XIAJ 114*9,14 


» . 


AFAPI^TR-67-125 

PaztIJ 


14  IF (XC91-XC71J 25*9*15 

15  A*{X{91— X{7 ) }/(X (81— X(6) 1 
Y*X{91-A*X(81 

IFIABSC10.*XC9)l-A0S(Yns,9,6 
C  ***  QUADRATIC  GUESS 

16  X(41»D£P 
X(5)*AIND 


IF(X(7)—X(5) 118,17*18 

17  IF(X{61-X(41U3,9,13 

18  lFrXC6)-X<4m9,13,19 

19  IF(Xi91-X(5>)23, 20,23 

20  IF(X(81— X(4)121*22*21 

21  XC9)*XC7) 

X(81«X(61 
GO  TO  13 

22  Xt9:*X(7l 
XfS}-*X  16} 

X(3)*l. 

lFCXC9JJ10ai»10 

23  1F(X C81— X(4J J24c21»24 

24  F*{Xlfe}-Xf41J/CXi7J-X{51) 
A*lX{8)-X{41-F*{XI9}-X{5))l/(CX{91-Xf71)*tX(9)-X{511) 
8*F»A*{X(5}*X<7*| 

C*X 141+X 15 J»IA*X (71— FI 


IF J A) 242*240*242 

240  IF (B1241 *7*241 

241  V—C/S 
CO  TO  37 

242  1F<8I247j243,247 

243  IFtCI245,244,245 

244  Y*0. 

CO  TO  37 

245  C*— C/A 

IF(G 17,7,246 

246  Y*SQRT(£l 
YY*-SQRT(GJ 
CO  TO  270 

247  IF (C 1249*248.249 


248  Y*-B/A 


YY*0* 


CO  TO  270 

249  D»4.*A*C/B**2 

IF 11,-0113,25 ,26 

25  Y*~ B/f 2.*A1 
CO  TO  37 

26  E*SQAT(1.-D) 

27  Y«t-B/f2.*A>)*tl.«-£J 
YY=(-8/t2»*All»U,-E) 

270  J*4 

DEPHIN*ABS(Xf411 
DO  29  1 *6*8,2 

1FCDEPHIN— ABSIXCI 1 1129*29,28 

28  J*I 
D£PHIN*ABS(XII 1 1 

29  CONTINUE 
K*J*1 

IFf  CX  CK1-Y1*CXCX1-YVJ 132,32,30 

30  IFCABSIXCX!— Y7— ABStXlK)— YY1 137*37*31 

31  Y*YY 

CO  TO  37 

32  IF (J-6133, 34,34 
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33  JJ*J+2 
XK*K+2 
GO  TO  35 

34  JJ*J-2 
KK*K-2 

35  SLO?E=(X(KK>— XJX) J/CXtJJJ-XCJi) 
IFiSLOPE*X(J)*IXCK}-YJ >36,36,37 

36  Y*YY 

37  *(9>=XC7) 

X(8)*X(6J 

Xf7J*X(5) 
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SIBFTC  FANDAT  DECK,M94/2,XR7 
BLOCK  DATA 

COHHON  /  FAN/CN115)fPR(15tl5),WAC{15,15i,£7A(15tl5)tN^HPClJ>; 

DATA  N,NP/10,6, 3*7,5*10, 8,5*0/ 

DATA  CN/0.3t0.4,0.5,Q»6,0.7,0.8,0.9,1.0,U1,1.2,5*0.^ 


DATA  IPR( 

1,J>,KAC( 

1 , J),£~a{ 

1,J),J=1,  6)/ 

A  1.000 

,  73.0 

,0.7300 

,  1.030  ,  69.0 

,0.7350 

* 

B  1.070 

4  60.0 

,0.7400 

,  1.096  ,  50.0 

♦C.7300 

9 

C  1.112 

,  40.0 

j0.790G 

,  1-120  ,  26.0 

,0.6200 

/ 

DATA  tPRt 

2,J)<HAC{ 

2, J) ,ETA{ 

2, J )yJ3li  7)/ 

D  1.000 

y  86.0 

,0.7300 

,  1.050  ,  81.0 

,0.7500 

9 

£  1.100 

,  76.0 

,0.7650 

,  1.146  ,  70.0 

,0.7700 

9 

F  1.188 

,  62.9 

,0.7750 

,  1.230  ,  55.0 

,0-7500 

9 

G  1.250 

,  47.0 

,0.7150 

/ 

DATA  {PR< 

3«J)*WAC( 

3, J } ,£TA( 

3,J),J=1,  7)/ 

H  l.OvO 

y iDO.O 

,0.7250 

,  1.064  ,  96.8 

,0.7500 

9 

I  1.128 

.  93.0 

,0.7750 

,  1-200  ,  87.5 

,0.8000 

9 

J  1.290 

y  78.0 

,0.8150 

,  1-330  ,  72-0 

,0.8000 

9 

K  1.370 

,  64.0 

,0.7500 

/ 

DATA  {PR I 

4* J } aMAC( 

4,J),£TA( 

4 . J 1 , J=1 ,  7)/ 

L  1.000 

,115.0 

,0.7500 

,  1-092  ,112-8 

<0.7500 

9 

K  1.220 

,107.5 

,0.8000 

,  1.310  ,102.0 

,0.8250 

■9 

N  1.400 

,  54.0 

,0.8500 

—  1—474  ,  63-0 

,0.8000 

■9 

0  1.488 

,  80.0 

,0.7800 

/ 

DATA  (PR( 

5,J),HAC( 

5,J),£TAl 

5,J1,J=1,1C)/ 

P  1.000 

,132-0 

,0.7000 

,  1-160  ,131-0 

,0.7500 

-9 

0  1.296 

,128.5 

,0.8000 

,  1.370  ,126.2 

,0.8250 

t 

R  1.460 

,122.0 

,0.8500 

,  1-524  ,118-0 

.0-8700 

“f 

S  1.544 

,116.5 

,0.8730 

.  1-560  ,115.0 

,0-8700 

-9 

T  1.610 

,110.5 

,0.2500 

,  1-668  ,102-8 

,0.8000 

/ 

DATA  (PR{ 

6„!),HAC( 

6, J ) ,ETA( 

:  J-1,10)/ 

y  1.000 

,250.0 

,0.6650 

,  1.250  ,149.9 

,9.7500 

f 

V  1.400 

,148.0 

,0.8000 

,  1.500  ,145.6 

,0-8250 

-t 

K  1.570 

,143.0 

,0.8500 

,  1.638  ,140.0 

,0.8700 

-f 

X  1.680 

,137.0 

,0.8800 

,  1.716  ,134.5 

.0-8709 

Y  1.756 

,130.0 

,G.8500 

,  1.830  ,122-0 

,  3.7950 

/ 

DATA  (PR I 

;  7, J),WAC( 

7,J),£TA(  7,J J, J=1,1D»/ 

Z  1.000 

,170.0 

t 0.620C 

,  1.190  ,170.0 

.9.7B09 

7f 

A  1.380 

,170.0 

,0.7500 

,  1.548  ,1-68,0 

,0.8080 

■s 

B  1.650 

,166.0 

,0.8300 

,  1.724  ,163-5 

,0-8590 

1 

C  1-823 

,15B.6 

,0.8700 

,  1.904  ,152.8 

.0-3590 

J 

D  1.978 

,145.0 

<0.8090 

.  2.0 DC  ,1-42  .2 

.9o790G 

/ 

DATA  fPR( 

I  8,JJ,«AC( 

!  8,J),ETA{  B,dI,J=l--iC)/ 

£  1.000 

,190.0 

,0-5200 

,1.110  ,150.0 

,9-6200 

F  1.338 

,150.0 

, 0. f 000 

,  1.552  ,139.9 

^0.7590 

-% 

G  1.736 

,187.8 

,0.8000 

,  1.850  jJL85-0 

.9-8390 

9 

H  2.000 

,280.0 

^0*8500 

,  2-070  ,176.0 

,0-8300 

9 

I  2.120 

♦173.0 

,0,390G 

,  2-2 00  -#1-66.0 

,0*7600 

/ 

HAT A  CPR{ 

[  S»Ji,£TA{  ®,JI,J*1,2.0]/ 

J  1.000 

-,210.0 

,0.5500 

,  1.260  .210.0 

.^-6209 

ft  1.550 

,210-0 

.0  -7000 

-v  1-310  .210-0 

,0-7599 

L  2.000 

,213.0 

,0.7750 

.  2.129  ,205-5 

.0<73C9 

A 

K  2.250 

,3C3-0 

,0.7750 

-  2*334  ,205-0 

,  0.-75  &3 

*> 

H  2.420 

-200.0 

.Q.-fjtDO 

,  2.460  ,19r.6 

-S-7D0S 

J 

*>*?*  {PIWJ.O*Ji>a6C{10vJ^£TA1iO»J}^4=l,  SJ/ 

0  1.000 

,225*0 

,  ,225—0 

■5 

P  l-TCV 

,225.0 

vD.7030 

,  2.090  ,225.0 

,9-7309 

* 

0  2.t£rt> 

♦225.0 

-i  5o350  ,225.0 

.0.7250 

■r 

2  2.450 

.-224.5 

*0.70 SO 

-r  2.650  ,221.0 

.0.6200 

f 
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SIBFTC  CMPDAT  0ECK,M94/2,XR7 
BLOCK  OATA 

COMMON  /  COHP/CN ( 15) ,PR ( 15, 15) ,MAC( 15,15) , ETA (15, 15 )  ,N,NP( 15) 
OATA  NyNP/10, 2*6, 2*8 .4*10, 2*8, 5*0/ 


DATA  CN/C 

•5  y0.6 f 0.7 

(0.8y0,89»0.92y0.95yl 

.0,1.1 

y 1 • 15 ,5*0.  / 

OATA  (PRC 

1 y J ) yMAC ( 

ly J) »ETA( 

1 , J ) ,  J*1  • 

6)/ 

A  1.000 

y  21.4 

y0.5700 

1.599  , 

21.1 

y 0.6000 

• 

B  2.018 

y  20.8 

y 0.6290 

2.335  * 

20.5 

•0.6300 

9 

C  3.024 

y  19.6 

,0.6500 

3.533  , 

18.7 

,0.6200 

OATA  (PR( 

2yJ)yWAC( 

2, JJyETAC 

2, J ) , J*1 y 

6)/ 

D  l.OCJO 

y  24.9 

,0.5700 

1.689  y 

24.9 

,0.6200 

9 

E  2.497 

y  24.7 

,0.6700 

3.096  y 

24.3 

,0.7000 

9 

f  3.636 

y  23.8 

,0.7200 

4.234  , 

23.1 

,0.7000 

/ 

OATA  (PR( 

3yJ)yWAC( 

3y J) y ETAC 

3y J ) , J“1 , 

8)/ 

G  1,000 

y  29.4 

,0.5650 

1.599  y 

29.4 

•0.6200 

9 

H  2.198 

y  29.4 

,0.6600 

2.719  y 

29.2 

,0.7000 

t 

1  3.515 

y  28.9 

,0.7500 

4.174  , 

28.5 

,0.7650 

9 

J  4.796 

y  27.9 

,0.7500 

4.946  y 

27.6 

,0.7430 

/ 

OATA  (PRC 

4,  J / «MAC( 

4y J) y ETA( 

4y  J) y J*1 y 

8)/ 

K  1.000 

y  35.6 

y 0.5600 

1.719  y 

35.6 

•0.6200 

9 

L  2.737 

,  35.6 

,0.7000 

3.575  , 

35.6 

,0.7500 

9 

N  4.353 

y  35.3 

y 0.78 00 

5.012  y 

35.0 

,0.8000 

9 

N  5.671 

y  34.3 

•0.7800 

5.970  y 

33.8 

,0.7650 

/ 

OATA  <PR( 

5  y J 1 •  WAC  ( 

5yJ)yETA( 

5,J),J-1,10)/ 

0  1.000 

y  42.7 

•0.5520 

2.084  , 

42.7 

,0.6200 

9 

P  3.132 

,  42.7 

,0.7000 

4.054  y 

42.7 

,0.7500 

9 

0  4.922 

*  42,5 

,0.7800 

5.431  y 

42.4 

,0.8100 

9 

R  6.000 

»  42.0 

.0.8300 

6.407  , 

41.8 

,0.8100 

9 

S  6.988 

,  41.2 

•0.7800 

7.138  y 

40.9 

,0.7780 

/ 

OATA  (PR( 

6yJ)»MAC( 

6yJ)yETA( 

6,J),J«1,10)/ 

T  1.000 

,  45.4 

,0.5500 

2.323  y 

45.4 

•0.6200 

9 

l)  3.353 

,  45.4 

,0.7000 

4.293  • 

45.4 

,0.7500 

9 

V  5.240 

.  45.2 

,0.7800 

5,731  y 

45.0 

,0.8100 

• 

M  6.389 

»  44.8 

•0.8300 

6.808  y 

44.5 

,3.8100 

9 

X  7.377 

y  43.9 

•0.7800 

7.587  f 

43.7 

,0.7760 

/ 

OATA  (PR ( 

7,j),WAC( 

7yJ)yETA( 

7, J 1 , J*1 ,  10  )  / 

Z  1.000 

«  48.1 

,0.5400 

2.617  y 

48.: 

,0.6200 

9 

A  3,683 

y  48.1 

,0.7000 

4.653  y 

48.1 

,0.7500 

9 

B  5.611 

y  48.1 

,0.7800 

6.222  • 

48.0 

,0.8100 

9 

C  6.868 

,  47.7 

•0.8200 

7.168  y 

47.6 

,0.8100 

9 

0  7.796 

y  47.3 

•0.7800 

8.133  y 

47.2 

,0.7700 

/ 

OATA  (PR( 

8yJ>yWAC( 

8, J), ETAC 

8,J)f J-1,10)/ 

E  1 .000 

y  51.6 

,0.5200 

1.491  y 

51.  .6 

,0.5500 

9 

F  3*03 6 

y  51.6 

,0.6200 

4.192  y 

51,6 

,0.7000 

f 

G  5.192 

y  5-1.6 

,0.7500 

6.263  y 

51.6 

,0.7800 

9 

H  6.689 

y  51.5 

,0.7950 

7.347  y 

51.3 

,0.8109 

9 

I  8.078 

y  51.1 

,0.7800 

8.725  y 

50.7 

,0.7500 

/ 

OATA  (PR( 

9*J)yWAC( 

9,J) yETAC 

9,J),J*1, 

8)/ 

J  1.000 

y  58.7 

y 0.4600 

3.395  y 

58.7 

,0.5800 

9 

K  5.731 

y  58.7 

,0.7000 

7.507  y 

58.7 

,0.7500 

9 

L  8.186 

y  58.7 

,0.7550 

0.784  y 

59.6 

,0.7500 

9 

M  9.683 

y  58.5 

,0.7000 

10.042  • 

58.4 

,0.6750 

/ 

OATA  ( PR ( 10,2)  ,WAC( 10, J  )  »ETA(  J 

Oy  J)  y  J**l  y 

B>/ 

N  1.000 

y  61.4 

,0.4500  , 

2.976  y 

6^.4 

,0.5500 

* 

0  4.916 

.  61.4 

• 0.6200  i 

6.838  y 

61.4 

,0.7000 

9 

P  7.587 

y  61.4 

,0.7250  , 

8.485  , 

61.4 

,0.7250 

9 

0  9.365 

y  61.4 

,0.7000  i 

10.581  y 

61.4 

,0.6200 

/ 

END 
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$IBFTC  CMBDAT  DECK,M94/2,X&7 

COMMON  /”cOM8/PSI  115  J  »DELT ( 15*15)  * ETAI 15 >1 5}  *N»Kf'  C 15  > 

DATA  N»NP/8, 8*11,7*0/  , 

DATA  PSI/10.*20.,30.*4Q.,50.*60.*i 

DATA  DEL 1/15*800., 15*900., 15*975., i  065. **5*11 7.,* *15  1260.* 
A15*1325 . , 15*1450 .,15*1550., 15*1685 .  - 5*1800. *60*0  ./ 

a820,0.9020,0.920C  ,0.9360,0.9500,0.9600,0.9680,T*0** 

IliliSliiSipsiiig; 

10  °00C,0  9195, 0.9376, 0.9543*0.9690, 0.9800,0.9850,0.9850,7*0., 

END 


J 


i 

] 

/ 

; 


i 

i 

t 


j 

* 

3 

I 

« 

-f 

1 

i 


/ 
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SIBFTC  HPTOAT  0ECK,M94/2,XR7 
BLOCK  DATA 

COMMON  /  HTUR8/TFF(15)  »CN(  15, 15) ,OH( 15,15) ,ETA (15,  15 )  ,N»NP(  151 


DATA  N ,NP/10»S*15» 12 ,5*0/ 


DATA  TFF/ 

16.86 

,18.27 

,20.17 

,20.66 

,20.90 

121.08 

.21.25 

,21.43 

,21.64 

•21.92 

,5*0.0 

DATA  (CN ( 

1,J),DH( 

1, J),ETA( 

1,J),J>1 

,15)/ 

AO. 191 

,0.0029 

,0.6150 

,0.344 

,0.0051 

,0.7000 

80.526 

,0.0076 

,0.7780 

,0.727 

,0.0097 

,0.8000 

C0.957 

,0.0120 

,0.8000 

,1.167 

,0.0137 

,0.7875 

01.340 

,0.0148 

,0.7692 

,1.569 

,0.0157 

,0.7340 

El. 761 

,0.0161 

,0.7000 

,1.971 

,0.0158 

,0.7550 

F2.193 

,0.0150 

,0.6000 

’,2.454 

,0.0130 

,0.5250 

G2.641 

,0.0108 

,0.4720 

•2.842 

,0.0074 

,0.4000 

H3.005 

,0.0031 

,0.3000 

/ 

DATA  (CN ( 

2, J) ,DH( 

2,J) ,ETA{ 

2, J ) , J*1 

,15)/ 

,0.7000 

10.191 

,0.0034 

,0.6000 

,0.402 

,0.0072 

JO. 593 

,0.0102 

,0.8000 

,0.727 

,0.0122 

,0.8200 

K0.861 

,0.0140 

,0.8270 

,1.000 

,0.0158 

,0.8300 

LI. 129 

,0.0173 

,0.8275 

,1.301 

,0.0191 

,0.8210 

Ml. 474 

,0.0205 

,0.8162 

,1.741 

,0.0223 

,0.8000 

N2.009 

,0.0234 

,0.7610 

•2.316 

,0.0235 

,0.7000 

02.751 

,0.0217 

,0.6000 

,3.158 

,0.0169 

,0.5000 

P3.445 

,0.0115 

,0.4150 

/ 

OATA  (CN ( 

3,J),0H( 

3,J),ETA( 

3, J ) ,  J*1 

,15)/ 

00.191 

,0.0041 

,0.5700 

,0.445 

,0.0090 

,0.7000 

R0.670 

,0.0130 

,0.8000 

,0.890 

,0.0166 

,0.8400 

S1.091 

,0.0194 

,0.8448 

,1.263 

,0.0216 

,0.8420 

T1.493 

,0.0241 

,0.8400 

,1-722 

,0.0263 

,0.8315 

U2.009 

,0.0234 

,0.8170 

•2.258 

,0.0298 

,0.8000 

V2.603 

,0.0310 

,0.7495 

,2.861 

•0.0311 

,0.7000 

M3. 100 

,0.0306 

,0.6578 

,3.330 

,0.0292 

,0.6000 

X3.445 

,0.028.1 

,0.5800 

/ 

* 

DATA  (CN( 

4,J),0H( 

4, J I ,ETA ( 

4,J), J-l 

,15)/ 

YO. 191 

,0.0047 

,0.5580 

,0.260 

,0.0061 

,0.6000 

Z0.488 

,0.0108 

,0.7000 

,0.708 

,0.0148 

,0.8000 

AC. 933 

,0.0184 

,0.8400 

,1.167 

,0.0220 

,0.8500 

31.416 

,0.0252 

,0.8500 

,1.593 

,0.0274 

,0.8480 

Cl. 837 

,0.0302 

,0.8440 

,2.019 

,0.0320 

,0.8400 

02.325 

,0.0349 

,0.8270 

,2.564 

,0.0371 

,0.8170 

E2.890 

,0.0397 

,0.8000 

,3.205 

,0.0  75 

,0.7710 

F3.445 

,0.0445 

,0.7500 

/ 

OATA  (CN( 

5, J) »DH( 

5, J)  ,ETA( 

5, J) » J*1 

,15)/ 

GO. 191 

,0.0050 

,0.5500 

,0.306 

,0.0079 

,0.6000 

HO. 536 

,0.0130 

,0.7000 

,0.765 

,0.0173 

,0.8000 

10.99:, 

,0.0212 

,0.8400 

,1.301 

,0.0259 

,0.8600 

J1.512 

,0.0289 

,0.8599 

,1.799 

,0.0324 

,0.8566 

K2.086 

,0.0360 

,0.8519 

,2.383 

,0.0400 

,0.8460 

L2.698 

,0.0446 

,0.8426 

,2.928 

,0.0486 

,0.8400 

M3. 138 

,0.0536 

,0.8400 

,3.215 

,0.0576 

,0.6437 

N3.225 

,0.0595 

,0.8475 

/ 

OATA  (CN( 

6,  J ) ,0H( 

6,J),ETA( 

6,J),J*1 

,15)/ 

00.191 

,0.0061 

,0.5250 

,0.364 

,0.0108 

,0.6000 

P0.632 

,0.0173 

,0.7000 

,0.880 

,0.0227 

,0  >8000 

01.115 

,0.0270 

,0.8400 

,1.311 

,0.0306 

>0.8600 

R1.531 

,0.0346 

,0.8721 

,1.722 

,0.0379 

,0.8800 

SI. 952 

,0.0425 

,0.8840 

,2.167 

,0.0*72 

,0.8870 

T2.316 

,0.0508 

,0.8875 

,2.471 

,0.0551 

,0.8876 

U2.545 

,0.0576 

,0.8868 

,2.588 

,0.0601 

,0.8838 

V2.607 

,0.0628 

,0.8798 

/ 

OATA  (CN 

7,  J ) ,0H ( 

7, J ) »ETA{ 

.7,  J) ,  J*1 

,15)/ 

* 

/ 
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HO. 191 

,0.0072 

,0.5005 

,0.440 

,0.0148 

,0.6000 

X0.698 

,0.0212 

,0.7000 

,0.976 

,0.0277 

,0.8000 

Y1.225 

,0.0335 

,0.8400 

,1.411 

,0.0374 

,0.8600 

Zl-541 

,0.0403 

,0.8700 

,1.651 

,0.0428 

,0.6800 

A1.780 

#0.0459 

,0.8855 

,1.899 

,0.0490 

,0.8900 

B1.995 

,0.0518 

,0.8910 

,2.040 

,0.0540 

,0.8900 

C2.086 

,0.0562 

,0.8880 

,2.124 

,0.0594 

,0.8826 

DZ.1**3  ,0-0630 

OATA  (CN(  J ) , OH ( 

,0.8695 
8,J2 ,ETA( 

/ 

8,J),J=1, 

151/ 

E0.191 

,0.0079 

,0.4995 

,0.493 

,0.0176 

,0.6000 

F0.746 

,0.0245 

,0.7000 

,0.899 

,0.0284 

,0.7580 

G1.043 

,0.0320 

,0.8000 

,1.167 

,0.0353 

,0.8200 

HI  .306 

,0.0389 

,0.8400 

,1.397 

,0.0414 

,0.8500 

11.493 

,0.0439 

,0.8600 

,1.627 

,0.0475 

,0.8710 

J1 .708 

,0.0504 

,0.8750 

i 1.780 

,0.0536 

,0.8750 

K1.837 

,0.0576 

,0.8690 

,1.852 

•0.0598 

,0.8600 

LI .856 
OATA  (CNI 

♦0,0624 

9,J),9HC 

,0.8495 

9,J),£TA( 

/ 

9,J»,J=1, 

151/ 

H0.191 

,0.0084 

,0.4855 

,0.344 

,0.0143 

,0.5320 

H0.545 

,0.0209 

,0.6000 

-0.689 

.0.0256 

,0.6500 

00.323 

,0.0297 

,0.7000 

,0.938 

,0.0331 

,0.7380 

PI. 033 

,0.0360 

,0.7690 

,1.148 

,0.0398 

,0.8000 

Q1.244 

,0.0432 

,0.8100 

,1.340 

,0.0472 

,0.8210 

R1.397 

,0.0500 

•0.8255 

,1.435 

•0.0522 

•0.8270 

•0.8150. 

Si. 474 

,0.0551 

,0.8230 

•1.493 

•0.0576 

T1.497 

,0.0601 

,0.8000 

/ 

DATA  CCNUe,JJ,DHC10,JJ,ETA{10,JI,J*l, 

121/ 

U0.191 

,0.0119 

•0.4210 

,0.287 

,0.0162 

,0.4695 

VO. 388 

,0.0205 

,0.5000 

,0.478 

•0.0241 

^0.5300 

H0.574 

,0.0283 

,0.5620 

»0i651 

,0.0317 

,0.5875 

X0.703 

,0.0342 

,0.6000 

,0.751 

,0.0371 

,0.6110 

Y0.785 

,0.0396 

,0.6170 

,0.823 

•0.0428 

,0.6240 

Z0.842 

END 

,0.0454 

,0.6195 

,0.847 

,0.0477 

•0.6050 
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AIBFTC  LPTDAT  DECK,M94/2,XK7 
BLOCK  DATA 


COMMON  /  LTURB/TFF(15),CN(15 
DATA  N,NP/il,9*15,12»9»4*0/ 
DATA  TFF/ 

*15I,0H(15 

•15) *ETA(15 

*15),N,MPC 

131.50 

*36.60 

,41.60 

,40.05 

*50.22 

*51.89 

253.41 

*54.60 

*55.69 

,56.89 

*58.10 

*4*0.0 

OATA  (CN( 

1*J)*DH( 

1*J>)*ETA( 

1*J)*J« 

1*15)/ 

AO. 365 

*0.0029 

*0.7120 

,0.529 

*0.0042 

*0.7300 

60.730 

*0.0056 

*0.7472 

,0.967 

*0.0070 

*0.7300 

Cl. 204 

*0.0081 

,0.7140 

,1.405 

*0.0089 

*0.7000 

01.606 

*0.0095 

*0.6850 

,1.752 

,0.0098 

*0.6730 

E2.C07 

*0.0100 

*0.6452 

*2.203 

*0.0098 

*0.6200 

F2.445 

*0.0092 

*0.6000 

*2.591 

*0.0085 

,0.5750 

62.792 

H3.2B5 

*0.0070 

*0.0001 

,0.5310 

,0.3850 

*2.920 

/ 

*0.0056 

,0.5000 

OATA  tCN( 

2*J),0H< 

2,J!,ETA< 

2* Jl, J* 

1* 

151/ 

10.365 

•0.0C42 

*0.8000 

,0.547 

*0.0063 

*0.8100 

JO. 785 

*0.0086 

,0.8200 

*1.058 

*0.0111 

*0.8300 

K1.277 

*0.0128 

,0.8300 

*1.432 

,0.0139 

*0.8290 

LI. 679 

*0.0153 

*0.8100 

*1.879 

*0.0162 

*0.8000 

M2. 044 

*0.0167 

,0.7850 

,2.208 

,0.0171 

*0.7600 

N2.354 

*0.0173 

,0.7450 

*2.600 

,0.0170 

,0.7000 

02 r 647 
P3.285 

*0.0162 

*0.0123 

,0.6800 

,0.5900 

,3.029 

/ 

*0.0151 

•0.6450 

OATA  (CNf 

3*J)*0M( 

3* Jl ,ETA( 

3,J),J- 

1* 

15)/ 

00.365 

*0.0050 

*0.8000 

*0.586 

*0.0081 

,0.8300 

R0.658 

*0.0114 

*0.8600 

,1.067 

*0.0139 

*0.8630 

SI. 241 

*0.0159 

,0.8670 

*1.423 

,0.0178 

,0.8700 

T1.606 

*0.0195 

,0.8720 

,1.825 

*0.0215 

•0.8720 

U2.007 

*0.0229 

*0.8700 

,2.226 

*0.0243 

,0.8670 

V2.409 

,0.0252 

,0.8600 

,2.573 

*0.0259 

,0.8500 

M2. 755 
X3.285 

,0.0265 

,0.0263 

,0.8300 

,0.7600 

,3.029 

/ 

,0.0268 

*0.8000 

OATA  (CN( 

4,J),0H( 

4»J| *ETA( 

4*  Jl* J*l* 

15)/ 

70.365 

*0.0053 

,0.7995 

,C.420 

*0.0061 

*0.8000 

Z0.675 

,0.0098 

,0.8400 

*0.876 

*0.0125 

*0.8600 

A1.095 

*0.0153 

,0.8680 

,1  277 

*0.0176 

•0.8730 

81.496 

,0.0201 

,0.8800 

♦1.733 

*0.0226 

*0.8830 

C2.025 

,0.0254 

,0.8835 

,2.190 

*0.0266 

*0.8830 

02.354 

,0.0279 

,0.8800 

,2.582 

*0.0293 

,0.8740 

E2.847 

F3.285 

,0.0305 

,0.0315 

*9.8500 

,0.8200 

,3.102 

/ 

•0.0312 

*0.8350 

OATA  (CMC 

5, J) *DH( 

5*J)*ETA( 

5, J) ,J* 

1, 

15)/ 

60.365 

,0.0057 

,0.7750 

*0.502 

*0.0078 

*0.8000 

HO. 730 

,0.0113 

,0.8480 

,0.967 

*0.0140 

*0.8600 

11.277 

,0.0191 

,0.8750 

*1.55) 

*0.0226 

*0.8900 

J 1.71 5 

,0.0248 

,0.8912 

,1.952 

,0.0276 

*0.8940 

K2.080 

,0.0290 

,0.8955 

*2.244 

,0.0307 

*0.8970 

L2.412 

*0.0323 

*0.8961 

,2.646 

,0.0343 

*0.8900 

K2.920 

M3. 285 

,0.0362 

*0.0382 

*0.8790 

*0.8600 

,3.157 

/ 

,0.0376 

,0.8671 

OATA  (CMC 

6, J ) ,0H  C 

6»J)  *ETA( 

6*  J)  *  J" 

1, 

15)/ 

00.36‘J 

,0.0061 

,0.7600 

,0.611 

*0.0103 

*0.8000 

P0.830 

,0.0139 

,0.8450 

,1.095 

*0.0181 

*0.8600 

01.277 

,0.0208 

*0.8730 

,1.515 

*0.0243 

*0.6900 

R1.733 

,0.0273 

,0.8950 

*1.934 

*0.0300 

*0.9000 

52.194 

*0.0334 

,0.9005 

,2.409 

*0.0362 

*0.9010 

T2.628 

,0.0390 

,0.9004 

,2.837 

*0.0415 

,0.9000 

U3.H1 

V3.285 

,0.0457 

,0.0510 

,0.8900 

,0.8735 

,3.230 

/ 

•0.0485 

,0.8800 
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DATA  (CMC 

?,J),OH( 

7,«J),ETA( 

7, J) » J*l»15)/ 

,0.8000 

s, 

HO. 365 

,0.0070 

,0.7310 

,0.766  ,0.0142 

* 

*  <•  .5 

XI. 004 

,0.0184 

,0.8300 

,1.255  ,0.0226 

#0.8600 

-t 

Y1.460 

,0.0259 

,0.8750 

,1.664  ,0.0290 

,0.8900 

* 

- 

Z1.825 

,0.0315 

,0.8930 

,2.007  ,0.0346 

,0.8975 

■t 

A2.171 

,0.0376 

,0.8999 

,2.280  ,0.0400 

$0.9000 

t 

r 

B2.372 

,0.0424 

,0.8980 

,2.482  ,0.0454 

,0.8937 

* 

C2.536 

02.609 

,0.0474 

,0.0527 

,0.8900 

,0.8710 

,2.591  ,0.0502 

/ 

,0.8799 

9 

DATA  tCNI 

8,J),DH( 

8,J) ,£TA( 

8»J),J*1»15)/ 

,0.7450 

K  *  £ 

E0.365 

,0.0077 

,0.7100 

,0.602  ,0.0125 

,»■ 

F0.785 

,0.0163 

,0.7680 

,0.967  ,0.0198 

,0.8000 

t 

1  5; 

G1.223 

,0.0245 

,0.8380 

, 1.442  ,0.0284 

,0.8600 

9 

HI. 624 

,0.0321 

,0.8712 

,1.825  ,0.0362 

#0.8780  - 

9 

:  s-3 

11.934 

,0.0387 

,0.8800 

,2.062  ,0.9418 

,0.8775 

9 

4L  -  i 

J2.135 

,0.0438 

,0.8760 

,2.190  ,0.0456 

,9.8722 

;~S 

— 

K2.244 

,0.0479 

,0.8660 

,2.281  ,0.0502 

,0.8600 

9 

L2.285 

,0.0514 

,0.8480 

/ 

DATA  ICHl 

9,J1,0H( 

9, J),ETA( 

9, J> »J=1,15}/ 

#0.7000 

- 

r  -  *< 

M0.365 

,0.0086 

,0.6780 

,0.547  ,0.0128 

•9 

’  _  !  « 

NO. 657 

,0.0153 

,0.7125 

,0.821  ,0.0190 

,0.7350' 

9 

i 

01.004 

,0.0226 

,0.7690 

,1.159  ,0.0256 

,0.8000 

■9 

:  ?i 

P1.223 

,0.0270 

,0.8060 

,1.369  ,0.0301 

,0.8225 

Q1.515 

,0.0334 

,0.8395 

,1.606  ,0.0357 

,0.8450 

-  "5 

R1.733 

,0.0390 

,0.8470 

tU8Z5  ,0.0421 

,0.8445 

'V  *-~ 

SI .379 

,0.0446 

,0.8330 

,1.898  ,0.0463 

,0.8235 

9  _ 

T1.907 

,0.0485 

,0.8080 

/ 

-* 

OATA  (CNI10,J J,OHf 10»J) ,ETAf 10»JI» J*l*12)/ 

■  '  . 

lfO.365 

,0.0098 

,0.6380 

,0.474  ,0.0125 

,0.6550  j: 

■9 

C 

V0.639 

,0.0167 

,0.6700 

,0.748  ,0.0195 

,0.6850 

9  ..  : 

*  * 

HO. 85b 

,0.0223 

,0.7000 

,0.967  ,0.0251 

,0.7110 

.#.  . 

v  -Cr 

£  *3 

XI. 095 

,0.0290 

,0.7180 

,1.191  ,0.0321 

,0.7180 

-.ti. 

Y1.259 

,0.0348 

,0.7170 

,1.296  ,0.0368 

#0.7140 

t  -  - 

Z1.325 

,0.0399 

,0.7000 

,1.329  ,0.0401 

,0.6890 

/ 

DATA  CCN!11, J),0HC11,J),£TA{11, J J,J*1,  91/ 

- 

% 

AO .365 

,0.0111 

,0.6000 

,0.438  ,0.0137 

,0.6000 

-9.  \  - 

80.547 

,0.0169 

,0.6120 

,0.639  ,0.0197 

,0.6170 

9-  - 

** "  1 

C0.730 

,0.0226 

,0.6210 

,0.821  ,0.0254 

,0.6258 

*  - 

00.876 

,0.0275 

,0.6259 

,0.931  ,0.0298 

,0.6230 

*  ■  V-:, 

E0.963 

END 

,0.032 2 

,0.6009 

/ 

--  • 

AFAPL-TR-67-125 

PtztS 


SO  AT  A 

•  fcmf 

ZHF 

IF 

9*F 

WAFC  „ 

MAF 

rate 

CMC 

1C 

FAC 

UACC 

MAC 

T2 

P2 

T2I  - 

P21 

T3 

P3  r. 

PCRLF 

H-iCF'v.::  -  , 

-  PCBLC 
BLC 

PCBLOB 
BLOB 
FCBLHF 
BLHP 
,  PCBLLP 
BLLP 
T* 

P* 

MA3 

MFB  - 
MCA 
"FARA 
ETA* 

DPCQH 
TFFMP 
CMMP 
DHTCHP 
OKTC 
T5 
P5. 

TFFLP 
CMLP 
OKTCLF 
OHTF  " 

T55 
P55 


PCBLOU 


AFAPL-TR-67-125 

Paitn 


ETAf 

ETAC 

ETA7HP 

ETATLP 

AM55 

AN25 

T6 

P6 

PS6 

AM  6 

V6 

UG6 

T7 

UFA 

VG7 

FAR7 

67AA 

DRAFT 

PS8 

AMS 

V8 

PS9 

AM9 

¥9 

PS2S 

AM2S 

Y2S 

PS29 

AM29 

¥29 

SYPAS5 

HPEXT 

WFT 

MGT 

VA 

FRO 

CVMNOZ 

VJM  \ 

CVDNOZ 

VJD 

FGM 

FGP 

THEEMO 

SOAIA IN  ITITLE*1 » 1PES*1,N00E*0,IDUNP*1 ,  1AMTP-0 « 

IGASKX-l,IOC0«Of  IMCD*0,N0ZFLT*0  »XTRYS*2i0,T0LALL*0 .005. 

ZFOS*O.833333,PC«FDS*lOOA,PAF0S*2.OO,ET£FDS*O.85O,lfAFDS*18O., 

ICDS*O.ei4598,PCNCOS*lOO.iPACDS*6.OO.ETACDS»O.83O,«ACDS*75.0, 

T4O5*24OO.,DTCOOS-126O.,ETAB0S*O.985,DPCOOS*©.C5r 

TFKP0S«21.25,CNHP0S*2.040,ETHP0S*P.S90. 

TFLPDS*53.41fCNLPDS*2.280,ETLP0S*0.9CC, 

DPOOOS-0,04  »0PAFDS«0.04t 

AM»0.0,AL7P*0.0*WP£XT*0.0.PC8LF*0.0f PC8LC*0.05»  _  . 

PC8LDU *0.20, PCS LOS *0.0»PC8LHP *©.80» PC8 LLP*0.0» 

?S55*1 .88  » CVONOZ*0.985, CVMN0Z*0.985 .  0€LFC*1 .0  ,0B.FN*1 ,0 ,DELSFC*1 .08 

SMOTE  DESIGN  POINT 

SOATAIN  T4=1750. ,AM*OoO»ALTP *00000., ITITLE-1S 

SLS  ISLE 

SDATAIN  T4-2400. • AM*0 .0, ALTP-OOOOO. , IT ITLE-18 

SLS  MILITARY 

SDATAIN  T4*2400. , AM*0 .0, ALTP *00000. , ITITLE*! , 

IAFTBN*1,T7*3200.»ETAA*0.91S 


AFAPL-TR-B7-125 
Part  II 


SLS  TAKE-OFF 

SOATAIN  T4*2400.*AK«1.2*ALTP*00500..ITITLE«1S 

SET-UP  LOW  ALTITUDE  DASH 
SOATAIN  T4*2400. , AH«1.2 •  ALTl’«00500* ,  IT  ITLE*1 , 
2AFTBN*i,T?*3200.*ETAA*0.91S 

LOU  ALTITUDE  DASH 

SDAtAIN  T4«2IOO»,AH*0.8,ALTP«36100.*imLS=IS 

SUBSONIC  CRUISE 

SDATAIN  T4»240Q. « AN*1 .2 »ALTP*50000. t ITITLE*1$ 

SUPERSONIC  AT  MILITARY  POWER 
S DATA IN  T4»2400. , AM*1.6,ALTP*50000. , ITITLE*1S 

SET-UP  SUPBISONIC  WITH  AFTERBURNER 
SDATAIN  T4*240G. * AM* 1.6 , ALTP*50000. , ITITLE-1 t 
IAFTBN»1,T7*1700.,ETAA«0.94S 

SUPERSONIC  PARTIAL  A/B 

SDATAIff  T4=2450«*AM»1< 4,AJ-TP*50000.  ,  ITXTLE*1 , 

I AFTBM»1 »T7*320G*»£TAA*0*9iS 

SUPERSONIC  FULL  A/B 


Security  Clarification 


DOCUMENT  CONTROL  DATA  -  R&D 

(Security  clmaeiticmtien  o/  title,  body  ot  abstract  and  Indexing  annotation  m tat  bo  entered  tehon  the  overall  report  la  do  a  milled) 


1  ORIGIN  A  TIN  G  ACTIVITY  (Corporate  author) 


Itc.  rs^lat  ikcurity  c  LAIIIRICATION 


Air  Force  Aero  Propulsion  Laboratory 
Wright- Patterson  Air  Force  Base,  Ohio  45433 


Unclassified 


3  REPORT  TITLE 

SIMULATION  OF  TURBOFAN  ENGINE 

Part  n.  User's  Manual  and  Computer  Program  Listing 


4  DESCRIPTIVE  MOTES  (Typo  ot  top art  mat  Inchitlro  d oft) 


5  AUTHOR^  ff-aal  ruMia.  lire!  name.  Initial) 

McKinney,  John  S. 

Captain,  USAF 


6  REPORT  DAYS 


November  1967 


4a.  CONTRACT  OR  ORANT  MO. 

b  PROJECT  MO.  3066 


c  Task  No.  306603 


•a.  oaiatNATOMu  report  mumbbrW 


AFAFL-TR-67-125,  Part  H 


»*.  gTHBR  REPORT  NOW  (Airy  a*  «r 


i  o  a  v *  il  ability/limit atiom  notice*  This  document  1b  subject  to  special  eaqr> 
transmittal  to  foreign  governments  or  foreign  nationals  may  be  made  only  with  prior 
of  the  Air  Force  Aero  Propulsion  Laboratory,  Wrlgfat-Patterson  AFB,  Ohio, 


I  1.  supplementary  notes 


13  ABSTRACT 


This  report  describes  a  digital  computer  program  titled  SMOTE  (Simulation  of  Tuibofan 
Engine).  SMOTE  is  a  computer  program  for  balancing-cycle  turbofan  pmgi«PM»  ^yt<iA  of 
turning  both  design  and  off-design  points.  The  program  Is  written  In  Fortran  TV  language 
nd  was  designed  for  use  on  an  IBM  7090  Digital  Computer,  although  it  lias  also  been  ms  on 
m  IBM  System  360.  Performance  maps  (Block  Data  format)  of  foe  .major  AAgiiiei  components 
ire  required.  Information  for  setting  up  the  Block  Data  and  Input  data  is  given  in  foe  report. 
Mso  included  in  the  report  is  a  complete  program  listing  with  a  description  of  each  sub¬ 
routine  and  a  sample  data  pack.  (  ) 

K 

(Distribution  of  foe  abstract  is  unlimited.) 


DD  /JK.  1473 


UNCLASSIFIED 
Security  Claaaification 


UNCLASSIFIED 


Security  CjurtficjUgi 


KEY  WOKOi 


Turbine  Engine  Simulation 
Tuxtotae  Engine  Computer  Programs 
Turbine  Engine  Cycle  Analysts 


LINK  Jl 


HOLE  I  WT 


-INK  I 


•VOLE  |  WT 


»  I _ ■  . . i- 

INSTRUCTIONS 


1.  ORIGINATING  ACTIVITY:  £m«  tbs  turn  and  address 
of  (Im  contractor,  subcontractor,  giant**,  Depart msnt  of  Da- 
feaaa  activity  or  othar  organization  'corporate  author)  issuing 
tha  report. 

2a.  REPORT  8ECUHTY  CLASSIFICATION:  Enter  tba  orar- 
all  security  classification  of  the  report.  Indicate  whether 
“Restricted  Data"  is  included.  Marking  is  to  be  in  accord¬ 
ance  with  appropriate  aacurity  regulations, 

2b.  GROUP:  Automatic  dowry adins  is  a  pacified  in  DoD  Di¬ 
rective  5200.10  and  Anted  Forces  Industrial  Manual.  Enter 
the  etoup  number.  Also,  whrn  applicable,  rhow  that  optional 
■makings  have  been  used  for  Croup  3  end  Group  4  an  author¬ 
ized. 

3.  REPORT  TITLE:  Enter  the  complete  report  title  in  all 
capital  letters.  Titlee  in  nil  cases  nhould  be  unclassified. 

If  a  meaningful  title  cannot  be  aalectsd  without  classifica¬ 
tion,  ahow  title  classification  In  all  capitals  in  paranthanis 
immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES:  II  appropriate,  enter  the  type  or 
report.  e.(. ,  interim,  progress,  aumma/y,  annual,  or  final. 

Give  the  inclusive  dates  when  a  specific  reporting  period  is 
covsred. 


5.  AUTHOR'S):  Enter  tbs  rams'*)  of  autboi<s)  ns  shown  no 
or  in  the  report.  Enfci  leet  name,  first  name,  middle  initial. 
If  military,  show  rank  and  branch  of  nervier.  The  name  of 
the  principal  author  is  an  absolute  minimum  requirement. 

6.  REPORT  DATE:  Enter  the  data  o'  the  report  as  day, 
month,  year,  or  month,  yean  If  mare  than  one  data  appears 
on  the  report,  use  dale  of  publication. 

7 a.  TOTAL  NUMBER  OF  PAAE&  The  total  page  count 
should  follow  normal  pagination  procedures,  be.,  enter  the 
number  uf  pages  containing  information, 

76.  NU“»ER  OF  REFERENCED  Enter  the  total  number  of 
references  cited  in  the  report. 

8  a.  CONTRACT  OH  GJtANT  NUMBER:  If  appropriate,  enter 
the  applicable  number  of  the  contract  or  grant  und«  which 
the  report  waa  written. 

8b,  8c,  b  td.  PROJECT  NUMBER:  Eater  the  appropriate 
military  department  identification,  such  *s  project  number, 
subproject  nuaber,  system  numbers,  task  number,,  etc. 

9 a.  ORIGINATOR'S  REPORT  NUMBER(S):  Enter  the  offi¬ 
cial  repot  t  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.  This  number  moat 


be  unique  to  this  report. 

So.  OTHER  REPORT  NUMBER'S;:  If  the  report  ha 
assigned  any  other  report  numbers  ( either  by  the  originator 
oi  by  the  sponsor!,  aict>  enter  this  number's). 


1C,  AVADUABn-lTY/UMlTATlON  NOTH.  <Dfc  fitter  any  lim¬ 
itations  on  'tether  dissesunation  uf  the  newt,  .  ■'hur  than  those 


imposed  by  aacurity  claaeWeatlcn,  using  standard  statements 
each  as: 

(lj  "Qualified  requesters  may  Obtain  copies  of  this 
report -&OM  DDO"  t 

(2)  "Foreign  announcement  and  dissemination  of  this 
report  by  DDC  is  not  authorized." 

(9)  "U.  8.  Uovorament  agencies  may  obtain  copies  of 

this  report  .dJroDlly  from  DDC.  Other  qualified  DDC 
users  shall  requrat  thmugh 

_  _ _  ’•* 

(4)  "U.  3  military  agencies  may  obtain  copies  of  this 

report  directly  tom  DDC.  Other  qualified  users 
shall  requeat  through 

»• 

(T)  "All  distribution  of  this  report  is  controlled.  Qual¬ 
ified  DDC  users  shall  request  through 

»* 

e  -  • 

If  the  i sport  has  bean  furnished  to  the  Office  of  Technical 
Bar  -ices,  Dsn  art  pent  of  Commerce,  for  sale  to  the  public,  indi¬ 
cate  this  iact  and  enter  the  price,  if  knows. 

1L  SUPPLEMENTARY  NOTES:  Use  for  additional  explana¬ 
tory  notes. 

12.  8P3N8G  iriO  MILITARY  ACTIVITY:  Enter  the  name  of 
the  dap  art  mental  project  office  or  laboratory  sponsoring  (pay- 
inf  tor)  tbs  rose  arch  and  development.  Include  address. 

13.  ABSTRACT:  Enter  on  abstract  giving  a  brief  and  factual 
comma -y  of  tba  document  indicative  of  the  report,  even  though 
It  may  also  appear  elsewhere  In  the  body  of  the  technical  re- 
port.  If  additional  space  is  required,  a  continuation  sheet  uhall 
be  attached. 

it  is  highly  desirable  that  the  abstract  of  classified  reports 
be  uaclacsifisd.  Each  paragraph  of  the  abstract  shall  end  with 
an  indication  cf  the  military  security  classification  of  the  in¬ 
formation  in  the  paragraph,  represented  as  fTS).  !S),  (C).  or  (V). 

There  le  no  limitation  on  the  length  of  the  abstract.  How¬ 
ever,  the  suggested  length  is  from  130  to  22S  words. 

14.  KEY  WOODS:  Key  words  are  technically  meaningful  terms 
or  short  phrases  that  characterise  a  report  and  may  be  used  as 
tadwx  entries  for  cataloging  the  report-  Key  words  must  be 
selected  so  that  no  aacurity  classification  is  required.  Identi¬ 
fiers,  such  as  equipment  model  designation,  trade  name,  military 
project  coda  name,  geographic  location,  may  be  us  ad  as  key 
words  but  will  be  followed  by  an  indication  of  technical  con¬ 
test.  The  assignment  of  links,  rules,  and  weights  is  optional. 


UNCLASSIFIED 


Security  Classification 


